Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Qm Math Orbital And Spin

1 · Prerequisites

2 · Summary

The orbital-and-spin layer of the non-relativistic quantum mechanics mathematical prerequisites, all of it mathematics: the finite-dimensional spin-s matrices of su(2), their symmetric-power and tensor representations and block decompositions, the Tsirelson bound for two binary sharp models with its singlet attainment, the Legendre and associated Legendre equations with the normalized spherical functions, their completeness in L2(S2,σ), the diagonal sphere Laplacian and weight operator with explicit domains, the harmonic-polynomial decomposition of polynomial spaces and of the isotropic oscillator levels, the strongly continuous unitary rotation action on L2(S2,σ), its self-adjoint orbital generators with the angular-momentum commutators and rotational domain identification, and the dipole angular selection rules. Every item carries domain: mathematics and depends only on mathematical suppliers, published supersets or items of qm-math-oscillator; nothing here asserts a physical interpretation of spin, measurement, transition rates or empirical evidence, and the companion page qm-math-orbital-and-spin-examples stays mathematics as well.

The spin layer begins with Finite spin matrices: for s∈12N the standard matrices Jzem=mem and J±em=(s∓m)(s±m+1) em±1 on C2s+1 satisfy J+∗=J−; the matrices Jz, Jx=(J++J−)/2 and Jy=(J+−J−)/(2i) are self-adjoint, obey the su(2) commutators and have Casimir J2=s(s+1)I. For s=12, Ji=12σi for i=x,y,z, with the multiplication law σiσj=δijI+i∑kεijkσk. Symmetric-power SU2 representation identifies the symmetric power of the spin-12 module with the spin-n2 module inside the n-fold tensor power; Finite spin-block decomposition splits any self-adjoint su(2) triple on a finite-dimensional space into spin blocks with the standard ladders and reads the multiplicities off dim⁡ker⁡(Jz−MI); and Tensor spin composition composes two spins, Vs⊗Vt=⨁j=∣s−t∣s+tWj with each Wj unitarily equivalent to the spin-j model and an explicit ladder construction. These finite-dimensional results feed Tsirelson bound: for bounded self-adjoint involutions the Bell operator C=A0⊗(B0+B1)+A1⊗(B0−B1) satisfies C2=4I−[A0,A1]⊗[B0,B1] and ∣⟨Cψ,ψ′⟩∣≤22∥ψ∥∥ψ′∥, the bound is attained by the singlet with the standard settings, and the explicitly specified common-outcome ±1 model obeys the CHSH bound 2 instead.

The angular layer is analytic. Legendre Rodrigues equation develops the Rodrigues polynomials Pl and the Legendre equation and orthogonality; Sphere integrals in angle coordinates identifies the polar surface measure with the chart measure sin⁡θ dθ dφ and computes σ(S2)=4π; Associated Legendre function normalization establishes the associated functions wlm, their equation and the normalized spherical functions Ylm=blmwlm(cos⁡θ)eimφ/2π for 0≤m≤l, with Yl,−m=(−1)mYlm‾ for 1≤m≤l; Weighted angular completeness proves completeness of (wlm)l in the weighted space; and Sphere ONB assembles the complete orthonormal family (Ylm) of L2(S2,σ). Sphere Laplacian diagonal domain defines the self-adjoint diagonal operators AYlm=l(l+1)Ylm and LzYlm=ℏmYlm on explicit square-summability domains with the finite span as a common core, matching the sphere differential expressions away from the poles and reducing the blocks Vl. Harmonic polynomial dimension connects the sphere to polynomials: Pl=Hl⊕r2Pl−2, dim⁡Hl=2l+1, the Ylm are restrictions of explicit harmonic homogeneous polynomials, and the isotropic oscillator level M contains exactly the orbital degrees M,M−2,⋯≥0, each once. Strongly continuous rotation action constructs the unitary action U(R)F=F∘R−1 with U(RS)=U(R)U(S), its preservation of the harmonic blocks with intertwining against the restriction isomorphisms, and the norm-continuity of every orbit map with the truncation estimate. Orbital self-adjoint generators extracts the self-adjoint generators Li of the coordinate-axis rotation groups, proves Li=−iℏ(x×∇)i on the blocks, the commutators [Li,Lj]=iℏ∑kεijkLk on the core, the identification L3=Lz, and the unitary identification of the three domains by explicit quarter turns. Finally Dipole angular selection rules expands Qqhlm into harmonic homogeneous pieces of degrees l+1 and l−1 and proves that the angular matrix element ⟨Yl′m′,QqYlm⟩ can be nonzero only for l′=l±1, m′=m+q: necessary conditions only, with no transition-probability claim. All items are draft mathematics under construction; the physical orbital and spin pages consume these results as stated mathematical premises and add their own postulates and interpretation.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

PropositionStatement: AI-adaptedProof: AI-adaptedOpen item page →

Finite spin matrices

Statement

Let s∈12N={0,12,1,32,… }, put N:=2s∈N and let Vs:=CN+1 carry the standard inner product ⟨z,w⟩=∑k=0Nzkwk‾, linear in the first argument (Real and complex inner product spaces, with the inner product linear in the first argument). Write (em)−s≤m≤s for the standard basis, indexed so that em is the vector with entry 1 in position k=s−m and entry 0 elsewhere; thus ⟨em,en⟩=1 if m=n and ⟨em,en⟩=0 otherwise. Define C-linear maps Jz,J+,J−∈End⁡C(Vs) by their values on this basis (Linear map between vector spaces over the same field):

Jzem:=m em,J+em:=(s−m)(s+m+1) em+1,J−em:=(s+m)(s−m+1) em−1,

where es+1:=0 and e−s−1:=0, and put

Jx:=12(J++J−),Jy:=12i(J+−J−),J2:=Jx2+Jy2+Jz2.

Then:

(a) Jz, Jx and Jy are self-adjoint endomorphisms of Vs (Self-adjoint and normal endomorphisms of a finite-dimensional real or complex inner product space), and J+∗=J−, J−∗=J+ (The adjoint T∗:W→V is characterised by ⟨Tv,w⟩W=⟨v,T∗w⟩V);

(b) on Vs the commutators are

[Jz,J+]=J+,[Jz,J−]=−J−,[J+,J−]=2Jz,[Jx,Jy]=iJz,[Jy,Jz]=iJx,[Jz,Jx]=iJy;

(c) J2=J−J++Jz2+Jz=J+J−+Jz2−Jz=s(s+1)I, and [J2,Jz]=[J2,J±]=0;

(d) the matrix elements of the raising operator are ∥J+em∥2=(s−m)(s+m+1) and ∥J−em∥2=(s+m)(s−m+1), and ⟨J+em,em+1⟩=(s−m)(s+m+1), which is nonzero exactly when m<s;

(e) for s=12, identify V12=C2 with the basis ordered as (e12,e−12) and let

σx:=(0110),σy:=(0−ii0),σz:=(100−1)

be the Pauli matrices. Then Ji=12σi for i=x,y,z, σx2=σy2=σz2=I, σxσy=iσz, σyσz=iσx, σzσx=iσy, and consequently σiσj+σjσi=2δijI and σiσj=δijI+i∑kεijkσk for all i,j, where δij=1 if i=j and 0 otherwise, and εijk=0 when two of i,j,k coincide, and otherwise εijk is the sign of the permutation (x,y,z)↦(i,j,k); the sum ranges over k∈{x,y,z}. The same convention with indices 1,2,3 uses the order-preserving relabelling (1,2,3)=(x,y,z).

All statements are finite-dimensional linear algebra: no physical interpretation of Jx,Jy,Jz as angular momentum, and no measurement postulate, is asserted here.

Facts & Assumptions

Given: s=N2 with N∈N; the space Vs=CN+1 with its standard inner product; the standard basis (em)−s≤m≤s; the maps Jz,J±,Jx,Jy,J2 displayed in the Statement.

[F1]

The standard basis is orthonormal: ⟨em,en⟩=δmn in the sense of the Statement, and every z∈Vs has the unique expansion z=∑m⟨z,em⟩em with coefficients ⟨z,em⟩; a linear map is determined by, and may be prescribed by, its values on a basis (Basis of a vector space: a linearly independent spanning subset; and ordered basis: an injective finite list whose image is a basis, Linear combination of a finite list, and the span span⁡(S) as the smallest linear subspace containing S, Linear map between vector spaces over the same field).

[F2]

For T∈End⁡(Vs) the adjoint T∗ is the unique endomorphism with ⟨Tz,w⟩=⟨z,T∗w⟩ for all z,w; T is self-adjoint when T=T∗, equivalently when ⟨Tz,w⟩=⟨z,Tw⟩ for all z,w (The adjoint T∗:W→V is characterised by ⟨Tv,w⟩W=⟨v,T∗w⟩V, Self-adjoint and normal endomorphisms of a finite-dimensional real or complex inner product space).

[F3]

Complex conjugation is additive, satisfies zw‾=z‾ w‾, z‾‾=z and zz‾=∣z∣2≥0, and zz‾=0 only for z=0; also i2=−1 and i‾=−i (The complex numbers as R[x]/(x2+1), with the real embedding and imaginary unit i, Conjugation is an involutive real-field automorphism, zz‾=∣z∣2, and modulus is definite, multiplicative, and subadditive).

[F4]

For real t≥0 the square root symbol denotes the unique nonnegative r with r2=t; the numbers (s−m)(s+m+1) and (s+m)(s−m+1) are nonnegative integers for every m∈{−s,−s+1,…,s}, and they vanish exactly at the endpoints m=s, respectively m=−s, where an endpoint vector is killed (The natural numbers N (von Neumann)).

[F5]

Matrix multiplication of the displayed 2×2 matrices is entrywise-finite multiplication; the identities σx2=σy2=σz2=I, σxσy=iσz, σyσz=iσx, σzσx=iσy are checked by direct multiplication using [F3].

Proof

1.1F1F2F3F4algebra

[Well-definedness; (a)] By [F4] the coefficients are nonnegative reals and the endpoint conventions are compatible, so prescribing the values Jzem,J±em on the basis defines unique linear maps Jz,J+,J−∈End⁡(Vs) by [F1]; Jx,Jy,J2 are then defined by the displayed formulas and are linear. For all basis vectors, ⟨Jzem,en⟩=mδmn=⟨em,Jzen⟩, and ⟨J+em,en⟩=(s−m)(s+m+1) δn,m+1, while ⟨em,J−en⟩=(s+n)(s−n+1) δm,n−1, which at n=m+1 equals (s+m+1)(s−m)=(s−m)(s+m+1), and both sides vanish otherwise; hence ⟨J+em,en⟩=⟨em,J−en⟩ and ⟨J−em,en⟩=⟨em,J+en⟩ for all m,n, so J+∗=J− and J−∗=J+ by [F1], [F2]. Taking adjoints in the definitions of Jx,Jy and using (J++J−)∗=J−+J+=J++J− and (12i(J+−J−))∗=12i‾(J+−J−)∗=1−2i(J−−J+)=12i(J+−J−) gives Jx∗=Jx and Jy∗=Jy. This is (a).

1.2F1F4algebra

[(b), ladder relations] For every basis vector, Jz(J+em)=Jz(s−m)(s+m+1)em+1=(m+1)J+em=J+(Jzem)+J+em, so acting on the basis and using linearity, [Jz,J+]=J+; similarly Jz(J−em)=(m−1)J−em=J−(Jzem)−J−em gives [Jz,J−]=−J−. Moreover, using the endpoint conventions, J+J−em=(s+m)(s−m+1)em and J−J+em=(s−m)(s+m+1)em, so [J+,J−]em=[(s+m)(s−m+1)−(s−m)(s+m+1)]em=2m em=2Jzem, and by linearity [J+,J−]=2Jz.

1.3F1F4algebra

[(d)] For −s≤m≤s, orthonormality [F1] and the definitions give ∥J+em∥2=(s−m)(s+m+1) and ∥J−em∥2=(s+m)(s−m+1), and ⟨J+em,em+1⟩=(s−m)(s+m+1); the coefficient vanishes exactly when m=s by [F4].

2.1step 1.2F3algebra

[(b), remaining commutators] From Jx=12(J++J−), Jy=12i(J+−J−) and the relations [J+,J−]=2Jz, [J−,J+]=−2Jz of step 1.2, expand the products bilinearly: [Jx,Jy]=14i([J+,J+]−[J+,J−]+[J−,J+]−[J−,J−])=14i(0−2Jz−2Jz−0)=−4Jz4i=iJz, because 1/i=−i by [F3]. Likewise [Jz,Jx]=12([Jz,J+]+[Jz,J−])=12(J+−J−)=iJy, and [Jz,Jy]=12i([Jz,J+]−[Jz,J−])=12i(J++J−)=−iJx, which is [Jy,Jz]=iJx. This completes (b).

2.2step 1.2F1algebra

[(c)] The definitions give Jx2+Jy2=12(J+J−+J−J+); combining this with [J+,J−]=2Jz from step 1.2 yields J2=J−J++Jz2+Jz=J+J−+Jz2−Jz. Direct expansion on each basis vector gives (J−J++Jz2+Jz)em=[(s−m)(s+m+1)+m2+m]em=s(s+1)em and (J+J−+Jz2−Jz)em=[(s+m)(s−m+1)+m2−m]em=s(s+1)em, using (s∓m)(s±m+1)=s2+s−m2∓m; hence J2=s(s+1)I by [F1]. Since J2 is a scalar multiple of the identity, it commutes with Jx,Jy,Jz and with J±=Jx±iJy, giving (c).

3.1F3F5algebra∎

[(e), Pauli case] Let s=12, so N=1 and the basis is (e12,e−12). From the definitions, Jze±12=±12e±12, J+e−12=e12, J+e12=0, J−e12=e−12, J−e−12=0; reading these actions as matrices in the ordered basis (e12,e−12) gives Jz=12σz, J+=12(σx+iσy) and J−=12(σx−iσy), whence Jx=12(J++J−)=12σx and Jy=12i(J+−J−)=12σy. The displayed matrices satisfy σi∗=σi by inspection. Direct multiplication using [F3] and [F5] gives σx2=σy2=σz2=I, σxσy=iσz, σyσz=iσx, σzσx=iσy; taking adjoints of the three off-diagonal products and using that the σi are Hermitian with i‾=−i gives σjσi=−σiσj for i≠j, so σiσj+σjσi=2δijI and σiσj=δijI+i∑kεijkσk for all ordered pairs i,j∈{x,y,z}. For i=j, every εiik is zero, and the product formula reduces to σi2=I; for i≠j, only the remaining coordinate k has nonzero εijk, with sign positive in the three cyclic orders and negative in their reversals, exactly matching the products just computed. This is (e).

The displayed relations are those of the standard (2s+1)-dimensional spin representation; the proposition alone asserts only the finite matrix identities written in its Statement.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Symmetric-power SU2 representation

Statement

Fix an integer n≥0, let V:=C2 carry the standard inner product linear in the first argument, with orthonormal basis (e+,e−) where e+:=(1,0) and e−:=(0,1) (Real and complex inner product spaces, with the inner product linear in the first argument), and let V⊗n be the n-fold algebraic tensor power with its finite string basis (The elementary tensors of two bases form the product basis of the tensor product) made orthonormal. Use Sn=Sym⁡({0,…,n−1}) as in The finite symmetric group Sn, one-line notation, and cycle notation. For n≥1 define the bijection bn:{0,…,n−1}→{1,…,n} by bn(a)=a+1, and for π∈Sn put π^:=bn∘π∘bn−1. Define Uπ to be the place action of π^ in Commuting symmetric-group and linear actions on a tensor power: explicitly, for x1,…,xn∈V, Uπ(x1⊗⋯⊗xn)=xπ−1(0)+1⊗⋯⊗xπ−1(n−1)+1, extended linearly. Use the diagonal action W⊗n of W∈M2(C) and the operators Δ(T):=∑j=1n1⊗(j−1)⊗T⊗1⊗(n−j). Put

P+:=1n!∑π∈SnUπ,Sn:=ran⁡P+⊆V⊗n,SU(2):={W∈M2(C):W∗W=I2, det⁡W=1},

where W∗:=W‾T is the conjugate transpose of W and I2 the identity matrix (Rectangular matrix multiplication and the identity matrix In, including zero-sized shapes, For n≥1, the determinant over a commutative ring by the Leibniz formula, and ∣det⁡A∣ for a real matrix). For W∈SU(2) define Un(W):=W⊗n∣Sn. For k∈{0,…,n} let sk be the sum of the elementary tensors eε1⊗⋯⊗eεn over all 2n strings (ε1,…,εn)∈{+,−}n containing exactly k entries −, and put vk:=(nk)−1/2sk.

For n=0 use V⊗0=C, the empty string s0=v0=1, 0!=(00)=1, the trivial group S0, Uπ=P+=IC, W⊗0=IC for every W, and Δ(T)=0. These conventions define S0=C and the trivial representation U0, with Jx=Jy=Jz=0.

Then:

(a) P+ is the orthogonal projection onto the closed subspace Sn, Uπ and W⊗n commute for every π∈Sn and every W, and W⊗n(Sn)⊆Sn; the assignment W↦Un(W) is well defined on all of M2(C), and Un(W′)Un(W)=Un(W′W) holds for all W,W′;

(b) SU(2) is a subgroup of the multiplicative monoid of M2(C) with identity I2; for W∈SU(2) the operator Un(W) is unitary with inverse Un(W−1)=Un(W)∗, so Un:SU(2)→U(Sn) is a group homomorphism (Strongly continuous unitary representations, invariant linear subspaces and intertwiners, Linear isometries, and orthogonal or unitary operators on finite-dimensional inner product spaces);

(c) the vectors v0,…,vn form an orthonormal basis of Sn; in particular Sn is (n+1)-dimensional;

(d) the operators Ji:=12Δ(σi)∣Sn, i=x,y,z, where σx,σy,σz are the Pauli matrices of Finite spin matrices, are self-adjoint and satisfy, with J±:=Jx±iJy, J2:=Jx2+Jy2+Jz2,

Jzvk=(n2−k)vk,J+vk=k(n−k+1) vk−1,J−vk=(k+1)(n−k) vk+1,

with v−1:=0 and vn+1:=0; moreover [Jz,J±]=±J±, [J+,J−]=2Jz, [Ji,Jj]=i∑kεijkJk and J2=n2(n2+1)ISn; the map vk↦en/2−k is a unitary equivalence of (Jx,Jy,Jz) with the spin-n2 matrices of Finite spin matrices;

(e) for each i the curve Ai(t):=cos⁡(t/2)I2−isin⁡(t/2)σi (t∈R) lies in SU(2), satisfies Ai(s)Ai(t)=Ai(s+t) and Ai′(0)=−i2σi, and t↦Vi(t):=Un(Ai(t)) is a strongly continuous one-parameter unitary group on Sn (Strongly continuous one-parameter unitary group) whose infinitesimal generator Gi (Infinitesimal generator of a unitary group) is Gi=−iJi on all of Sn; equivalently 1t(Vi(t)v−v)→−iJiv for every v∈Sn;

(f) Un is strongly continuous on SU(2): if Wm→W0 in M2(C) with all Wm∈SU(2) and W0∈SU(2), then ∥Un(Wm)−Un(W0)∥→0 in the operator norm on Sn.

This is finite-dimensional linear algebra: the theorem constructs a unitary representation of SU(2) and its generators. It does not assert that SU(2) or any spin observable describes a physical particle, and it makes no measurement or statistics claim.

Facts & Assumptions

Given: n≥0; the space V=C2 with its orthonormal basis (e+,e−); the tensor power V⊗n; the place permutations Uπ and diagonal actions W⊗n; the projection P+; the subspace Sn; the strings sk and vectors vk; the group SU(2); the Pauli matrices σx,σy,σz; the operators Δ(T) and Ji.

[F1]

For n≥1, iterating The elementary tensors of two bases form the product basis of the tensor product gives the finite basis of all 2n strings eε1⊗⋯⊗eεn; for n=0 the basis is 1∈C. Define the inner product by ⟨∑cεeε,∑dεeε⟩=∑cεdε‾. Expansion on elementary tensors gives the product inner product. The finite tensor product is thus the standard finite-dimensional coordinate inner product space; no completion or countable choice is used. The universal multilinear property is Finite iterated tensor products represent multilinear maps independently of parenthesization.

[F2]

For n≥1, conjugation by bn is a bijection from the zero-indexed Sn of The finite symmetric group Sn, one-line notation, and cycle notation to the group on tensor positions {1,…,n} of Commuting symmetric-group and linear actions on a tensor power; it preserves composition because π^σ^=bnπσbn−1=πσ^, and π^−1(j)=π−1(j−1)+1 for 1≤j≤n. Thus the stated Uπ are well-defined linear maps by [F1], permute the orthonormal string basis, are unitary, satisfy UπUσ=Uπσ and have Uπ∗=Uπ−1. There are n! permutations: prescribing the images of 0,…,n−1 successively leaves n,n−1,…,1 possibilities. In the finite average P+=(1/n!)∑πUπ, inversion and left or right multiplication by a permutation reindex the summands. Thus P+∗=P+, UπP+=P+Uπ=P+ and P+2=P+; its range consists exactly of the invariant vectors, since every invariant vector v satisfies P+v=v. For n=1, π^(1)=1 and Uπ=IV; for n=0, all identities follow from the separate empty-tensor conventions.

[F3]

On the spanning string basis, expansion of matrix entries gives (S⊗T)(S′⊗T′)=(SS′)⊗(TT′), and evaluation of the inner products of [F1] gives (S⊗T)∗=S∗⊗T∗. Thus a tensor product of unitary matrices is unitary. Applying the same W to every factor commutes with any place permutation, since both rules send v1⊗⋯⊗vn to the tensor of the same permuted list Wvj; this identity is also valid for n=0 (Commuting symmetric-group and linear actions on a tensor power).

[F4]

For a bounded operator L on a finite-dimensional Hilbert space, the operator norm is ∥L∥=sup⁡∥z∥≤1∥Lz∥, and ∥L∥≤∑j∥Lej∥ for every orthonormal basis (ej); a linear map with finite-dimensional domain is bounded, and ∥Lz∥≤∥L∥ ∥z∥ (The operator norm as the least bound and as the unit-sphere or unit-ball supremum, A bounded linear operator between normed spaces). In orthonormal product coordinates, writing z=∑bzb⊗eb gives ∥(S⊗I)z∥2=∑b∥Szb∥2≤∥S∥2∑b∥zb∥2; similarly ∥I⊗T∥≤∥T∥. Since S⊗T=(S⊗I)(I⊗T) and ∥LMz∥≤∥L∥∥M∥∥z∥, this proves ∥S⊗T∥≤∥S∥∥T∥, and iteration supplies the tensor estimates below.

[F6]

The Pauli matrices satisfy σi∗=σi, σx2=σy2=σz2=I2, σxσy=iσz, σyσz=iσx, σzσx=iσy; also σ++σ−=2σx and σ+−σ−=2iσy for σ+:=σx+iσy, σ−:=σx−iσy, and σ+e−=2e+, σ+e+=0, σ−e+=2e−, σ−e−=0, σze±=±e± (Finite spin matrices).

[F8]

For integers 0≤k≤n, (nk)=n!/(k!(n−k)!)>0 and (nk)=(nn−k) and (nk−1)/(nk)=k/(n−k+1) for 1≤k≤n; the ratio follows by cancelling the factorials using k!=k(k−1)! and (n−k+1)!=(n−k+1)(n−k)!. Each string with k minus signs corresponds bijectively to its set of minus positions, so there are (nk) such strings and the string sums split accordingly ((nk) k! (n−k)!=n! for k≤n; hence (nk) k!=nk‾, the quotient n!/(k!(n−k)!) is a natural number, and (nk)=(nn−k), The factorial n! and the falling factorial nk‾, defined by recursion in N, The set [A]k of k-element subsets and the binomial coefficient (nk):=∣[n]k∣, A finite sum in a commutative monoid indexed by an arbitrary finite set).

[F10]

For s=n2 the operators Jz,J±,Jx,Jy,J2 of Finite spin matrices act on their orthonormal basis em, −s≤m≤s, by Jzem=mem, J+em=(s−m)(s+m+1)em+1, J−em=(s+m)(s−m+1)em−1, with endpoint conventions es+1=e−s−1=0, have self-adjoint Jx,Jy,Jz and J−=J+∗, satisfy [Jz,J±]=±J±, [J+,J−]=2Jz and J2=s(s+1)I, and have ∥J+em∥2=(s−m)(s+m+1).

Proof

1.1givenF10algebra

[The zero tensor power] For n=0, the conventions give P+=I, S0=Cv0, U0(W)=I, and Ji=0. Clauses (a)-(c) follow immediately; all ladder and Casimir formulas in (d) vanish, and v0↦e0 identifies the zero operators with the spin-zero model [F10]. Each Vi(t)=I has derivative zero and is strongly continuous, and U0 is constant, proving (e)-(f); the subgroup and Ai assertions, independent of n, are proved below. For the tensor and basis calculations in the remaining steps assume n≥1.

1.2F2F3algebra

[The symmetric subspace and invariance] By [F2] P+ is an orthogonal projection and UπP+=P+; hence P+v=v exactly for v∈Sn, and Sn=ker⁡(I−P+) is closed because the finite matrix I−P+ is continuous. For W∈M2(C) and π∈Sn, [F3] gives W⊗nP+=1n!∑πW⊗nUπ=1n!∑πUπW⊗n=P+W⊗n, so v∈Sn implies P+(W⊗nv)=W⊗nP+v=W⊗nv, i.e. W⊗n(Sn)⊆Sn; therefore Un(W):=W⊗n∣Sn is a well-defined linear map on Sn for every W.

1.3F6F7F9algebra

[The curves Ai] Using I2∗=I2, σi∗=σi, σi2=I2, [F7] and [F9]: Ai(t)∗Ai(t)=cos⁡(t/2)2I2+sin⁡(t/2)2σi2−icos⁡(t/2)sin⁡(t/2)(σi−σi∗)=(cos⁡2+sin⁡2)(t/2)I2=I2, and direct computation of the three 2×2 determinants from the displayed matrices gives det⁡Ai(t)=cos⁡2(t/2)+sin⁡2(t/2)=1; so Ai(t)∈SU(2). The identities cos⁡(s+t)=cos⁡scos⁡t−sin⁡ssin⁡t and sin⁡(s+t)=sin⁡scos⁡t+cos⁡ssin⁡t, applied with σi2=I2, give Ai(s)Ai(t)=Ai(s+t). Differentiating the displayed formula at t=0 with [F7] and [F9] gives Ai′(0)=−i2σi.

2.1step 1.2F3F5algebra

[Representation law and unitarity] For W,W′∈M2(C), [F3] gives (W′)⊗nW⊗n=(W′W)⊗n on V⊗n, hence on Sn by step 1.2 and restriction: Un(W′)Un(W)=Un(W′W). If W∈SU(2), then W is unitary by [F5], [F6]; by [F3] W⊗n is unitary on V⊗n, and for v∈Sn, ∥Un(W)v∥=∥W⊗nv∥=∥v∥ with W⊗n mapping Sn bijectively onto Sn (its inverse is (W−1)⊗n restricted, where W−1=W∗∈SU(2)); hence Un(W) is unitary with (Un(W))∗=Un(W−1). This gives the operator statements of (a) and (b).

2.2step 1.2F1F2F8algebra

[The basis vk and the dimension] Fix k and let w=eε1⊗⋯⊗eεn have exactly k minus entries. In n!P+w=∑πUπw every string with exactly k minus entries occurs once for each of the k!(n−k)! permutations carrying the pattern of w to it, so n!P+w=k!(n−k)! sk and P+w=(nk)−1sk by [F8]; since the 2n elementary tensors span V⊗n [F1], this shows Sn=ran⁡P+⊆span⁡{s0,…,sn}, while sk∈Sn because each Uπ permutes the strings with k minus entries, so P+sk=sk. Distinct elementary tensors are orthonormal [F1], and two summands from sk and sl with k≠l are distinct elementary tensors, so ⟨sk,sl⟩=δkl(nk) by [F8]; hence v0,…,vn are orthonormal and span Sn, an orthonormal basis of the (n+1)-dimensional space Sn.

2.3step 1.2F2F3F6algebra

[The operators Δ(σi) and their commutation] For T,T′∈M2(C) the place sums satisfy [Δ(T),Δ(T′)]=Δ([T,T′]) and Δ(T)∗=Δ(T∗): expanding the definitions, the summands at different places commute and the summand at a common place contributes [T,T′], while the adjoint of a place sum is the place sum of the adjoints [F3]. Consequently [Δ(σi),Δ(σj)]=Δ([σi,σj]); and Δ(σi) commutes with every Uπ because permuting places preserves the multiset of places, so by [F2] each Δ(σi) preserves Sn, and its restriction Ji=12Δ(σi)∣Sn is self-adjoint.

3.1step 2.1F5F9algebra

[SU(2) is a subgroup] If W,W′∈SU(2), then (WW′)∗(WW′)=W′∗W∗WW′=W′∗W′=I2 and det⁡(WW′)=det⁡Wdet⁡W′=1, so WW′∈SU(2); I2∗I2=I2 and det⁡I2=1, so I2∈SU(2); and for W∈SU(2) one has W∗W=WW∗=I2 (multiplying W∗W=I2 by W−1 on the right gives W∗=W−1, whence WW∗=I2) and det⁡W∗=det⁡W‾=1, so W−1=W∗∈SU(2). Hence SU(2) is a subgroup, and Un is a group homomorphism by step 2.1.

3.2step 2.1step 2.3step 1.3F1F3F4F6F7algebra

[Generator of t↦Un(Ai(t))] By step 2.1, Vi(t):=Un(Ai(t)) satisfies Vi(s)Vi(t)=Vi(s+t) and Vi(0)=I, and each Vi(t) is unitary. The telescoping identity W⊗k−W0⊗k=∑j=1k(W⊗(j−1)⊗(W−W0)⊗W0⊗(k−j)), valid by [F3] with W⊗0:=I, gives ∥W⊗k−W0⊗k∥≤kmax⁡(∥W∥,∥W0∥)k−1∥W−W0∥ by [F4]; with W=Ai(t), W0=Ai(t0) and ∥Ai(⋅)∥=1, the explicit formula for Ai and [F7] show that ∥Vi(t)v−Vi(t0)v∥ is bounded by n∥Ai(t)−Ai(t0)∥ ∥v∥ for all v∈Sn by the operator-norm bound; so Vi is a strongly continuous one-parameter unitary group. For the derivative, write Ai(t)=I2+tAi′(0)+r(t) with ∥r(t)∥≤∣t∣η(t) and η(t)→0 by [F6], [F7] and step 1.3; the same telescoping expansion with W=Ai(t), W0=I2, using ∥Ai(t)∥≤1 and ∥r(t)∥≤∣t∣η(t), gives Ai(t)⊗n−I2⊗n=t∑j=1n1⊗(j−1)⊗Ai′(0)⊗1⊗(n−j)+R(t) with ∥R(t)∥≤n∣t∣ η′(t) and η′(t)→0 as t→0; restricting to Sn and using step 2.3 and Ai′(0)=−i2σi of step 1.3 gives 1t(Vi(t)v−v)→−i2Δ(σi)v=−iJiv for every v∈Sn. Hence the infinitesimal generator is Gi=−iJi with domain D(Gi)=Sn.

3.3step 2.2step 2.3F1F6F8algebra

[Actions on the basis vectors] Fix 0≤k≤n. On an elementary tensor with k minus entries, Δ(σz) acts by #+−#−=n−2k, and by linearity Δ(σz)sk=(n−2k)sk, so Jzvk=(n/2−k)vk. For σ+=σx+iσy one has σ+e+=0 and σ+e−=2e+ by [F6]. If 1≤k≤n, each target string with k−1 minus entries in Δ(σ+)sk occurs with coefficient 2 from each of its n−k+1 plus positions, so Δ(σ+)sk=2(n−k+1)sk−1 and J+vk=(n−k+1)(nk−1)/(nk)vk−1=k(n−k+1)vk−1 by [F8]. For k=0, every factor is e+, so J+v0=0, agreeing with the stated v−1=0 convention. Similarly, if 0≤k<n, σ−e−=0 and σ−e+=2e− give Δ(σ−)sk=2(k+1)sk+1, hence J−vk=(k+1)(nk+1)/(nk)vk+1=(k+1)(n−k)vk+1 by [F8]. For k=n, every factor is e−, so J−vn=0, agreeing with vn+1=0.

4.1step 3.3F10algebra

[Commutators and identification with the spin-n/2 model] By steps 2.3 and 3.3, Jz,J± satisfy [Jz,J±]=±J± and [J+,J−]=2Jz on the basis (vk); substituting Jx=12(J++J−), Jy=12i(J+−J−) and expanding gives [Jx,Jy]=iJz, [Jy,Jz]=iJx, [Jz,Jx]=iJy. Since vk↦en/2−k maps the basis of Sn onto the basis of the spin-n2 space of Finite spin matrices, and by step 3.3 the three operators have there exactly the action prescribed in that proposition, the described map is a unitary equivalence intertwining (Jx,Jy,Jz) with the spin-n2 matrices and carrying J2 to s(s+1)I with s=n/2; in particular J2=n2(n2+1)ISn. This proves (a)-(d) except the strong-continuity clause.

5.1step 1.1step 1.2step 2.1step 3.1step 3.2F3F4algebra∎

[Strong continuity] Let Wm→W0 in M2(C). By the telescoping identity of step 3.2 with [F3], ∥Wm⊗n−W0⊗n∥≤nmax⁡(∥Wm∥,∥W0∥)n−1∥Wm−W0∥, and max⁡(∥Wm∥,∥W0∥) is bounded because ∥Wm∥→∥W0∥ by [F4] and convergence of the entries; hence ∥Wm⊗n−W0⊗n∥→0. Restricting the operators to the invariant subspace Sn can only decrease the operator norm, so ∥Un(Wm)−Un(W0)∥→0, proving (f); in particular Un is strongly continuous on SU(2), since norm continuity implies strong continuity by [F4]. This completes (a)-(f).

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Finite spin-block decomposition

Statement

Let V be a nonzero finite-dimensional complex inner product space, with inner product linear in the first argument (Real and complex inner product spaces, with the inner product linear in the first argument, Finite-dimensional vector space, and its dimension dim⁡FV; infinite-dimensional means having no finite basis), and let Jx,Jy,Jz∈End⁡(V) be self-adjoint operators satisfying (Self-adjoint and normal endomorphisms of a finite-dimensional real or complex inner product space)

[Jx,Jy]=iJz,[Jy,Jz]=iJx,[Jz,Jx]=iJy.

Put J+:=Jx+iJy, J−:=Jx−iJy and J2:=Jx2+Jy2+Jz2. Then:

(a) J+∗=J−, J−∗=J+, [Jz,J±]=±J±, [J+,J−]=2Jz, J2=J−J++Jz2+Jz=J+J−+Jz2−Jz, and J2 is self-adjoint, commutes with Jx,Jy,Jz,J±, and satisfies ⟨J2v,v⟩=∥Jxv∥2+∥Jyv∥2+∥Jzv∥2≥0 for every v∈V;

(b) there are an integer r≥1, numbers j1,…,jr∈12N={0,12,1,… } with ∑ν=1r(2jν+1)=dim⁡V, and an orthonormal basis (vν,k), ν∈{1,…,r}, k∈{0,1,…,2jν}, of V such that

Jzvν,k=(jν−k)vν,k,J+vν,k=k(2jν−k+1) vν,k−1,J−vν,k=(k+1)(2jν−k) vν,k+1,

J2vν,k=jν(jν+1)vν,k,Jxvν,k=12(J++J−)vν,k,Jyvν,k=12i(J+−J−)vν,k,

with vν,−1:=0 and vν,2jν+1:=0; each summand Vν:=span⁡{vν,0,…,vν,2jν} is Jx,Jy,Jz-invariant, the decomposition V=⨁ν=1rVν is orthogonal, and the map vν,k↦ejν−k is a unitary equivalence of (Jx,Jy,Jz) on Vν with the spin-jν matrices of Finite spin matrices;

(c) let λ be the largest eigenvalue of Jz (it exists and is real), and let u≠0 satisfy Jzu=λu. Then λ≥0, J+u=0, and with wk:=J−ku one has wk≠0 for 0≤k≤2λ, wk=0 for k=2λ+1, 2λ∈N, and ∥wk∥2=k!(2λ)(2λ−1)⋯(2λ−k+1)∥u∥2 for 0≤k≤2λ+1; the vectors wk/∥wk∥, 0≤k≤2λ, form the block of (b) with j=λ; consequently λ∈12N and the largest eigenvalue of Jz on each block is one of the numbers j1,…,jr;

(d) for every M∈12N the number of blocks with jν∈M+Z and jν≥M equals dim⁡ker⁡(Jz−MI), and the number of blocks with jν=M equals dim⁡ker⁡(Jz−MI)−dim⁡ker⁡(Jz−(M+1)I); in particular the multiset {j1,…,jr} is determined by the pair (V,Jz) together with the relations, and the spectrum of Jz is the multiset union of the numbers jν−k, 0≤k≤2jν (a block with parameter j contributes its weights only to the lattice coset j+Z).

No physical interpretation of the operators as angular momentum, spin or observables is asserted; the theorem is a classification of finite-dimensional triples satisfying the displayed commutation relations.

Facts & Assumptions

Given: a nonzero finite-dimensional complex inner product space V; self-adjoint Jx,Jy,Jz satisfying the three commutation relations; the derived operators J± and J2.

[F1]

Adjoints are conjugate-linear, (T∗)∗=T, (T+T′)∗=T∗+T′, (TT′)∗=T′∗T∗, and a restriction of a self-adjoint operator to a closed invariant subspace is self-adjoint; self-adjoint operators are normal (The adjoint T∗:W→V is characterised by ⟨Tv,w⟩W=⟨v,T∗w⟩V, Self-adjoint and normal endomorphisms of a finite-dimensional real or complex inner product space).

[F2]

Every self-adjoint endomorphism of a nonzero finite-dimensional complex inner product space has a real eigenvalue and a unit eigenvector, so its set of eigenvalues is a nonempty finite subset of R and has a largest element (Complex spectral theorem: a normal endomorphism of a finite-dimensional complex inner product space has an orthonormal eigenbasis, and conversely).

[F3]

For a self-adjoint T and eigenvectors u,v with Tu=λu, Tv=μv and λ≠μ one has ⟨u,v⟩=0: (λ−μ)⟨u,v⟩=⟨Tu,v⟩−⟨u,Tv⟩=0; the orthogonal complement of an invariant subspace is invariant for a self-adjoint operator; and ⟨T2v,v⟩=∥Tv∥2 for self-adjoint T (Orthogonality and the orthogonal complement, Self-adjoint and normal endomorphisms of a finite-dimensional real or complex inner product space, Real and complex inner product spaces, with the inner product linear in the first argument).

[F4]

The norm satisfies ∥v∥2=⟨v,v⟩≥0, ∥v∥=0 iff v=0, and ⟨Tv,w⟩=⟨v,T∗w⟩ for all v,w (The induced length is a norm, The adjoint T∗:W→V is characterised by ⟨Tv,w⟩W=⟨v,T∗w⟩V).

[F5]

An orthogonal decomposition V=W⊕W⊥ has dim⁡V=dim⁡W+dim⁡W⊥, and a finite orthogonal family of nonzero vectors is linearly independent (Orthogonality and the orthogonal complement, Finite-dimensional vector space, and its dimension dim⁡FV; infinite-dimensional means having no finite basis, Linear subspace of a vector space).

[F6]

The matrices of Finite spin matrices for s=j act on an orthonormal basis em, −j≤m≤j, by the formulas displayed in the Statement with s=j, are self-adjoint, satisfy the same commutation relations, and have J2=j(j+1)I (Finite spin matrices).

[F7]

For integers 0≤k≤m the products k!(2λ)(2λ−1)⋯(2λ−k+1) vanish if and only if 2λ is a nonnegative integer strictly less than k; sums of dimensions over a finite orthogonal direct sum add (The natural numbers N (von Neumann), A finite sum in a commutative monoid indexed by an arbitrary finite set, Finite-dimensional vector space, and its dimension dim⁡FV; infinite-dimensional means having no finite basis).

Proof

1.1F1F3F4algebra

[The derived algebra] Taking adjoints of J±=Jx±iJy with Jx∗=Jx, Jy∗=Jy and the conjugation laws gives J+∗=J− and J−∗=J+. The commutation relations give [Jz,J±]=[Jz,Jx]±i[Jz,Jy]=iJy±i(−iJx)=±(Jx±iJy)=±J± and [J+,J−]=[Jx+iJy,Jx−iJy]=−i[Jx,Jy]+i[Jy,Jx]=2Jz. Direct expansion gives J−J+=Jx2+Jy2−Jz=J2−Jz2−Jz and J+J−=Jx2+Jy2+Jz=J2−Jz2+Jz, hence J2=J−J++Jz2+Jz=J+J−+Jz2−Jz. Since J± and J2 are polynomials in the self-adjoint Ji, [F1] gives J2∗=J2. Finally [J2,Jx]=[Jy2,Jx]+[Jz2,Jx]=Jy[Jy,Jx]+[Jy,Jx]Jy+Jz[Jz,Jx]+[Jz,Jx]Jz=−i(JyJz+JzJy)+i(JzJy+JyJz)=0, and cyclically [J2,Jy]=[J2,Jz]=0, so J2 commutes with Jx,Jy,Jz and hence with J±; and ⟨J2v,v⟩=⟨Jx2v,v⟩+⟨Jy2v,v⟩+⟨Jz2v,v⟩=∥Jxv∥2+∥Jyv∥2+∥Jzv∥2≥0 by [F3], [F4]. This proves (a).

2.1step 1.1F2F3F4algebra

[A highest vector] By [F2] Jz has a largest eigenvalue λ∈R and a unit eigenvector u with Jzu=λu, where u≠0 because ∥u∥=1. From [Jz,J+]=J+ of step 1.1, Jz(J+u)=J+Jzu+J+u=(λ+1)J+u, so J+u is either 0 or an eigenvector of Jz with eigenvalue λ+1; maximality of λ forces J+u=0. Applying step 1.1 to u gives ⟨J2u,u⟩=⟨(J+J−+Jz2−Jz)u,u⟩=⟨J+J−u,u⟩+λ2−λ, while ⟨J2u,u⟩≥λ2 by (a); equivalently, since J+u=0, J+J−u=J+J−u−J−J+u=[J+,J−]u=2Jzu=2λu and therefore ⟨J2u,u⟩=2λ+λ2−λ=λ(λ+1)≥∥Jzu∥2=λ2, so λ≥0.

3.1step 1.1step 2.1F3F4F5algebra

[The string and its length] Define wk:=J−ku for k∈N, so w0=u. Inductively, Jzwk=(J−Jz−J−)J−k−1u=(λ−k)wk by step 1.1, and using J+J−=J2−Jz2+Jz, the fact that J2 commutes with J− and J2u=λ(λ+1)u from step 2.1, one computes J+wk=J+J−wk−1=(J2−Jz2+Jz)wk−1=[λ(λ+1)−(λ−k+1)2+(λ−k+1)]wk−1=k(2λ−k+1)wk−1 for k≥1. Pairing with wk−1 and using [F4] gives the recursion ∥wk∥2=k(2λ−k+1)∥wk−1∥2, hence ∥wk∥2=k!(2λ)(2λ−1)⋯(2λ−k+1)∥u∥2 by induction, the case k=0 being the empty product. If wk≠0 for every k, the vectors wk are nonzero eigenvectors of Jz for the pairwise distinct eigenvalues λ,λ−1,λ−2,…; any dim⁡V+1 of them form an orthogonal family of nonzero vectors, hence are linearly independent by [F3] and [F5], which is impossible in the finite-dimensional space V. So there is a least m≥1 with wm=0, and then 0=∥wm∥2=m(2λ−m+1)∥wm−1∥2 forces 2λ=m−1∈N.

4.1step 1.1step 2.1step 3.1F3F6algebra

[The normalized block] Put 2j:=m−1=2λ with j=λ∈12N, and for 0≤k≤2j set vk:=wk/∥wk∥, which is defined because wk≠0 for k<m=2j+1 and unit by construction. The eigenvectors wk have distinct eigenvalues j,j−1,…,−j of the self-adjoint Jz, so they are pairwise orthogonal by [F3]; hence v0,…,v2j is an orthonormal family. The action formulas follow from steps 1.1, 2.1 and 3.1: Jzvk=(j−k)vk; J−vk=∥wk+1∥/∥wk∥ vk+1=(k+1)(2j−k) vk+1 for k<2j and J−v2j=0; J+vk=k(2j−k+1)∥wk−1∥/∥wk∥ vk−1=k(2j−k+1) vk−1 for k≥1 and J+v0=0; and J2vk=λ(λ+1)vk=j(j+1)vk because J2 commutes with J−. The span Vλ:=span⁡{v0,…,v2j} is therefore invariant under J± and Jz, hence under Jx,Jy,Jz, and [F6] identifies the action with the spin-j model, giving Jxvk=12(J++J−)vk and Jyvk=12i(J+−J−)vk.

5.1step 4.1F1F3F5algebra

[The invariant complement] Let W:=Vλ. If w∈W⊥ and w′∈W, then for each i the invariance of W and self-adjointness of Ji give ⟨Jiw,w′⟩=⟨w,Jiw′⟩=0, so Jiw∈W⊥; thus W⊥ is invariant under Jx,Jy,Jz. The restriction of each Ji to W⊥ is self-adjoint by [F1], the commutation relations hold there, and V=W⊕W⊥ is an orthogonal direct sum with dim⁡V=dim⁡W+dim⁡W⊥ by [F5]; if W⊥≠{0}, the restricted triple again satisfies the hypotheses of the Statement.

6.1step 2.1step 3.1step 4.1step 5.1F5F6F7algebra

[Exhaustion and the decomposition] Iterate step 5.1: as long as the orthogonal complement of the blocks already produced is nonzero, step 2.1 supplies a largest eigenvalue and step 4.1 a further nonzero invariant block of dimension 2jν+1≥1, orthogonally split off, decreasing the dimension of the complement by 2jν+1. Since the dimension is a nonnegative integer, the process terminates after some finite number r≥1 of blocks with complement {0}, yielding an orthogonal decomposition V=⨁ν=1rVν and an orthonormal basis (vν,k) as in (b); the basis is orthonormal as a union of the block bases because the blocks are pairwise orthogonal, and ∑ν(2jν+1)=dim⁡V by additivity of dimensions in orthogonal direct sums [F5], [F7]. Each block action is the spin-jν action by step 4.1 and [F6], which proves (b); step 3.1 applied to a block identifies the largest eigenvalue of Jz on that block as jν, completing (c) except for λ≥0, proved in step 2.1.

7.1step 6.1algebra∎

[Multiplicity count] Let M∈12N and consider the eigenvalue-M eigenspace of Jz. A block with parameter j contributes the basis vector with k=j−M exactly when j∈M+Z and 0≤M≤j, and contributes nothing in all other cases (its weights are j,j−1,…,−j, all in the coset j+Z); hence dim⁡ker⁡(Jz−MI)=#{ν:jν∈M+Z, jν≥M} by the orthonormal basis of step 6.1. Applying this at M and at M+1, both in the same coset M+Z, and subtracting gives #{ν:jν=M}=dim⁡ker⁡(Jz−MI)−dim⁡ker⁡(Jz−(M+1)I). The multiset {j1,…,jr} is therefore determined by Jz (using the fixed relations to know that the decomposition exists), and the spectrum statement follows by collecting the eigenvalues jν−k of the blocks; the coset qualifier in (d) excludes nothing, since blocks outside M+Z contribute no vector to either eigenspace. This proves (d).

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Tensor spin composition

Statement

Let s,t∈12N, write Ls:={−s,−s+1,…,s} and Lt:={−t,−t+1,…,t}, and let (Vs,Jx(s),Jy(s),Jz(s)) and (Vt,Jx(t),Jy(t),Jz(t)) be the spin models of Finite spin matrices with orthonormal bases (em(s))m∈Ls and (en(t))n∈Lt and the ladders J±(s):=Jx(s)±iJy(s), J±(t):=Jx(t)±iJy(t) and Casimirs J(s)2, J(t)2 (The elementary tensors of two bases form the product basis of the tensor product). On the tensor product Vs⊗Vt put

Jitot:=Ji(s)⊗IVt+IVs⊗Ji(t)(i=x,y,z),J±tot:=Jxtot±iJytot,Jtot 2:=(Jxtot)2+(Jytot)2+(Jztot)2.

Then:

(a) each Jitot is self-adjoint, the triple Jxtot,Jytot,Jztot satisfies the commutation relations of Finite spin-block decomposition, and Jtot 2=J(s)2⊗IVt+IVs⊗J(t)2+2∑iJi(s)⊗Ji(t) is self-adjoint and commutes with Jxtot,Jytot,Jztot;

(b) there is an orthogonal decomposition Vs⊗Vt=⨁j=∣s−t∣s+tWj into Jtot-invariant subspaces, one for each j∈{∣s−t∣,∣s−t∣+1,…,s+t}, with dim⁡Wj=2j+1, on each of which (Jxtot,Jytot,Jztot) is unitarily equivalent to the spin-j model of Finite spin matrices; consequently ∑j=∣s−t∣s+t(2j+1)=(2s+1)(2t+1);

(c) for every M∈12Z the M-eigenspace of Jztot has dimension N(M):=#{m∈Ls: M−m∈Lt}={max⁡(0, min⁡(s,M+t)−max⁡(−s,M−t)+1),M−(s+t)∈Z,0,M−(s+t)∉Z., and the identity (em(s)⊗en(t))m,n is an orthonormal basis of Vs⊗Vt consisting of eigenvectors of Jztot with eigenvalues m+n;

(d) the number of summands with parameter equal to M in the decomposition of (b) is N(M)−N(M+1) if M≥0, and this number equals 1 when M−(s+t)∈Z and ∣s−t∣≤M≤s+t, and equals 0 otherwise; hence the multiset of parameters in (b) is exactly {∣s−t∣,∣s−t∣+1,…,s+t}, each occurring once;

(e) for each j∈{∣s−t∣,…,s+t} the space of vectors of Jztot-eigenvalue j orthogonal to all summands Wj′ with j′>j is one-dimensional; for a unit vector uj in it the vectors vj,k:=(k!(2j)(2j−1)⋯(2j−k+1))−1/2(J−tot)kuj(0≤k≤2j) form an orthonormal basis of Wj with Jztotvj,k=(j−k)vj,k, J+totvj,k=k(2j−k+1) vj,k−1 and J−totvj,k=(k+1)(2j−k) vj,k+1; the coefficients of vj,k in the product basis are therefore determined by the phase of uj, and choosing that phase fixes the Clebsch-Gordan numbers explicitly; in particular us+t=es(s)⊗et(t) is a permissible choice.

The theorem is a purely mathematical statement about tensor products of the explicit matrices of Finite spin matrices; it asserts no physical coupling, conservation law or measurement fact.

Facts & Assumptions

Given: s,t∈12N; the spin models Vs,Vt with their operators and orthonormal bases; the tensor product Vs⊗Vt; the operators Jitot, J±tot, Jtot 2.

[F1]

By The elementary tensors of two bases form the product basis of the tensor product, the algebraic tensor product Vs⊗Vt has the finite basis em(s)⊗en(t). Give its coordinate vectors the inner product ⟨∑cmnem⊗en,∑dmnem⊗en⟩=∑cmndmn‾. This is positive definite, makes that basis orthonormal, and expansion of two elementary tensors gives ⟨u⊗v,u′⊗v′⟩=⟨u,u′⟩⟨v,v′⟩. This is a finite construction with explicitly supplied bases; no completion or countable choice is invoked.

[F2]

Define S⊗T by its values Sem⊗Ten on the finite product basis. Bilinear expansion gives (S⊗T)(u⊗v)=Su⊗Tv, and hence (S⊗T)(S′⊗T′)=(SS′)⊗(TT′). The inner product of [F1] gives ⟨(S⊗T)(u⊗v),u′⊗v′⟩=⟨u⊗v,(S∗⊗T∗)(u′⊗v′)⟩, proving (S⊗T)∗=S∗⊗T∗ on the spanning product basis. Also I⊗I=I, so [S⊗I,S′⊗I]=[S,S′]⊗I, [I⊗T,I⊗T′]=I⊗[T,T′], and operators acting on different factors commute.

[F3]

On Vs and Vt the operators are self-adjoint, satisfy [Jx,Jy]=iJz and its cyclic companions, and have J2=s(s+1)IVs, respectively t(t+1)IVt; the ladder actions on the bases are J±(s)em(s)=(s∓m)(s±m+1)em±1(s) and Jz(s)em(s)=mem(s) (Finite spin matrices).

[F4]

Every finite-dimensional triple of self-adjoint operators satisfying the commutation relations decomposes orthogonally into spin blocks with ∑ν(2jν+1)=dim⁡V, the block actions being the spin-jν actions; the number of blocks with parameter M≥0 equals dim⁡ker⁡(Jz−MI)−dim⁡ker⁡(Jz−(M+1)I), and M must lie in 12N (Finite spin-block decomposition (b), (d)).

[F6]

For a self-adjoint operator the eigenspaces for distinct eigenvalues are orthogonal, and for a unit vector w and self-adjoint T, ⟨Tw,w⟩=⟨w,Tw⟩ is real; ∥w∥2=⟨w,w⟩ (Self-adjoint and normal endomorphisms of a finite-dimensional real or complex inner product space, Orthogonality and the orthogonal complement, The induced length is a norm).

[F7]

Factorials and products: k!=∏i=1ki≥1 for k≥0, and a finite product of positive numbers is positive (The factorial n! and the falling factorial nk‾, defined by recursion in N, The natural numbers N (von Neumann)).

Proof

1.1F2F3algebra

[Self-adjointness and the algebra] For i=x,y,z, [F2] and [F3] give (Ji(s)⊗I)∗=Ji(s)⊗I and (I⊗Ji(t))∗=I⊗Ji(t), so Jitot is self-adjoint. Adding the two-sided commutators, [Jxtot,Jytot]=[Jx(s),Jy(s)]⊗I+I⊗[Jx(t),Jy(t)]=iJz(s)⊗I+iI⊗Jz(t)=iJztot, and cyclically; so the triple satisfies the commutation relations.

1.2F1F3algebra

[Weights] By [F1] and [F3], Jztot(em(s)⊗en(t))=(m+n)em(s)⊗en(t), so the elementary tensors form an orthonormal eigenbasis of Jztot with eigenvalues m+n, and the number of basis vectors of eigenvalue M is the number of pairs (m,n)∈Ls×Lt with m+n=M, which is #{m∈Ls:M−m∈Lt}=N(M); all sums m+n lie in s+t+Z, so N(M)=0 off that lattice. On the lattice, the interval endpoints lie in the same unit-spaced lattice Ls; the displayed closed form follows by counting the interval [max⁡(−s,M−t),min⁡(s,M+t)], which is empty exactly when the maximum exceeds the minimum, and otherwise contributes one element per unit lattice step.

2.1step 1.1F2F3algebra

[The Casimir] Expanding Jtot 2=∑i(Ji(s)⊗I+I⊗Ji(t))2 gives J(s)2⊗I+I⊗J(t)2+2C with C:=∑iJi(s)⊗Ji(t), which is self-adjoint by [F2] and [F3]; moreover [Jtot 2,Jitot]=[J(s)2,Ji(s)]⊗I+I⊗[J(t)2,Ji(t)]+2[C,Ji(s)⊗I+I⊗Ji(t)] vanishes because the first two commutators are 0 by [F3] and 2[C,Ji(s)⊗I+I⊗Ji(t)]=2i∑k,lεkil(Jl(s)⊗Jk(t)+Jk(s)⊗Jl(t))=0, the parenthesised tensor being symmetric in k,l while εkil is antisymmetric.

2.2step 1.2F1algebra

[The parity-qualified difference] Restrict now to M≥0. If M−(s+t)∉Z, neither M nor M+1 is a weight and both counts are zero by step 1.2. On the compatible lattice, exchange the factors if necessary to suppose s≥t. Counting the interval in step 1.2 gives N(M)=min⁡(2t+1,s+t−M+1) for 0≤M≤s+t, and N(M)=0 for M>s+t. Therefore N(M)−N(M+1)=0 on the plateau 0≤M<s−t, equals one for s−t≤M≤s+t, and is zero above s+t, including the upper endpoint where N(s+t)=1 and N(s+t+1)=0. Restoring factor symmetry gives precisely the admissible lattice points in [∣s−t∣,s+t].

3.1step 1.1step 2.1step 2.2F4F5algebra

[Existence and parameter list] By steps 1.1, 2.1 and [F4] the triple (Jxtot,Jytot,Jztot) decomposes Vs⊗Vt into an orthogonal direct sum of spin blocks with parameters jν∈12N, and by F4 the number of blocks with parameter M equals dim⁡ker⁡(Jztot−MI)−dim⁡ker⁡(Jztot−(M+1)I)=N(M)−N(M+1) for M≥0, which by step 2.2 is 1 precisely for M∈{∣s−t∣,∣s−t∣+1,…,s+t} and 0 otherwise. Hence the block parameters are exactly the numbers ∣s−t∣,∣s−t∣+1,…,s+t, each once; each block has dimension 2j+1 and its operators are unitarily equivalent to the spin-j model by F4. The dimension identity ∑j=∣s−t∣s+t(2j+1)=(2s+1)(2t+1) follows from ∑ν(2jν+1)=dim⁡(Vs⊗Vt) and [F5]. This proves (b) and (d).

4.1step 3.1F1F3F4F6F7algebra∎

[The constructive bases and Clebsch-Gordan phase] Fix j in the parameter list and let Ej be the space of Jztot-eigenvectors of eigenvalue j orthogonal to all blocks Wj′ with j′>j. Counting as in step 3.1, dim⁡ker⁡(Jztot−jI)=#{ν:jν−j∈Z, jν≥j} (a block with parameter jν≥j≥0 and jν−j∈Z contributes exactly one eigenvector of eigenvalue j), and the higher blocks with parameters j′>j contribute one each, so dim⁡Ej=#{ν:jν−j∈Z, jν≥j}−#{ν:jν−j∈Z, jν>j}=1. Pick a unit uj∈Ej, which is determined up to multiplication by a phase; applying the ladder identities of F4 in the block normal form, the vectors vj,k of the Statement are well defined and form the orthonormal basis of Wj, with the stated actions, because ∥(J−tot)kuj∥2=k!(2j)(2j−1)⋯(2j−k+1) by [F7] and the recursion of F4; expanding vj,k in the product basis [F1] exhibits the Clebsch-Gordan coefficients, and replacing uj by eiθuj multiplies them all by the same phase. For j=s+t, es(s)⊗et(t) has weight s+t and is annihilated by J+tot because J+(s)es(s)=J+(t)et(t)=0 by [F3]; this is the top vector of the unique block of parameter s+t, so it is a permissible choice. This proves (e) and completes the proof.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Tsirelson bound

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let H,K be complex Hilbert spaces (Hilbert space) and let A0,A1∈B(H), B0,B1∈B(K) be bounded self-adjoint operators (A bounded linear operator between normed spaces, Self-adjoint, positive, unitary and normal operators) with Ai2=IH and Bj2=IK (The Hilbert-space adjoint of a bounded operator). Let H⊗^K be the completed tensor product of Hilbert tensor construction: inner product, completion and product bases and put

C:=A0⊗(B0+B1)+A1⊗(B0−B1)∈B(H⊗^K).

Then:

(a) C is self-adjoint and C2=4 IH⊗^K−[A0,A1]⊗[B0,B1];

(b) ∥C2∥≤8 and ∥C∥≤22 in the operator norm (The operator norm as the least bound and as the unit-sphere or unit-ball supremum, Hilbert-adjoint identities);

(c) for all ψ,ψ′∈H⊗^K, ∣⟨Cψ,ψ′⟩∣≤22 ∥ψ∥ ∥ψ′∥ (this is the Tsirelson bound for the stated pair of binary sharp models); consequently the absolute value of the CHSH combination of the four expectations ⟨(Ai⊗Bj)ψ,ψ⟩ of a unit vector ψ is at most 22;

(d) the bound is attained on H=K=C2: let σx,σy,σz be the Pauli matrices of Finite spin matrices, let Ψ:=12(e+⊗e−−e−⊗e+) with the orthonormal basis (e+,e−) of C2, put σa:=axσx+ayσy+azσz for a unit vector a∈R3, and take A0:=σz, A1:=σx, B0:=σ(z+x)/2, B1:=σ(z−x)/2; then ⟨σa⊗σbΨ,Ψ⟩=−a⋅b for all unit a,b, and ⟨CΨ,Ψ⟩=−22, so ∣⟨CΨ,Ψ⟩∣=22;

(e) a comparison with the stated local model: if in addition A0,A1,B0,B1 are four ±1-valued random variables on one probability space, then the absolute value of their CHSH combination is at most 2 (The CHSH bound and the finite-sampling law for independent trials (a)); this is a statement about that explicitly specified common-outcome model, not about all theories, and it is not compared here with any reported experiment.

The theorem is a mathematical statement about bounded operators; no measurement, preparation or empirical claim is asserted.

Facts & Assumptions

Given: Countable Choice; complex Hilbert spaces H,K; self-adjoint involutions A0,A1∈B(H), B0,B1∈B(K); the completed tensor product H⊗^K; the operator C.

[F1]

The completed tensor product is a complex Hilbert space containing the algebraic tensor product densely, and the elementary tensors u⊗v have inner products ⟨u⊗v,u′⊗v′⟩=⟨u,u′⟩⟨v,v′⟩; ∥u⊗v∥=∥u∥ ∥v∥ (Hilbert tensor construction: inner product, completion and product bases, Hilbert space).

[F2]

For bounded S,S′∈B(H), T,T′∈B(K): (S⊗T)(S′⊗T′)=(SS′)⊗(TT′), S⊗IK and IH⊗T commute, (S⊗T)∗=S∗⊗T∗, and ∥S⊗T∥≤∥S∥ ∥T∥; the identity operator on H⊗^K is IH⊗IK (Bounded tensor maps: extension, norm estimate, adjoint and commutation).

[F3]

For bounded operators on Hilbert spaces, the adjoint is additive and conjugate-linear, (TS)∗=S∗T∗, T∗∗=T, ∥T∗∥=∥T∥ and ∥T∗T∥=∥T∥2 under Countable Choice (Hilbert-adjoint identities). From The operator norm as the least bound and as the unit-sphere or unit-ball supremum, ∥TSx∥≤∥T∥∥Sx∥≤∥T∥∥S∥∥x∥; taking the supremum over the unit ball gives ∥TS∥≤∥T∥∥S∥. Likewise the vector triangle inequality gives ∥T+S∥≤∥T∥+∥S∥. These estimates include zero spaces.

[F4]

For a bounded operator T and all ψ,ψ′, ∣⟨Tψ,ψ′⟩∣≤∥T∥ ∥ψ∥ ∥ψ′∥ by Cauchy–Schwarz and the norm definition (Cauchy–Schwarz: ∣⟨u,v⟩∣≤∥u∥∥v∥, with equality exactly for linearly dependent vectors, The operator norm as the least bound and as the unit-sphere or unit-ball supremum). A self-adjoint involution satisfies T∗T=TT∗=I, so it is unitary in the Hilbert-space sense (Self-adjoint, positive, unitary and normal operators).

[F5]

The Pauli matrices are self-adjoint, σi2=I2, σxσy=iσz and cyclically, and σa for a unit vector a is self-adjoint with σa2=I2, because σa2=∑i,jaiajσiσj=∑iai2I2+∑i≠jaiajσiσj=∥a∥2I2 by the anticommutation relations σiσj+σjσi=2δijI2 (Finite spin matrices).

[F6]

If A0,A1,B0,B1:Ω→{−1,1} are random variables on a probability space, then ∣E[A0B0]+E[A0B1]+E[A1B0]−E[A1B1]∣≤2, with no choice principle used (The CHSH bound and the finite-sampling law for independent trials (a)).

[F7]

A bounded self-adjoint involution satisfies ∥Tx∥2=⟨Tx,Tx⟩=⟨x,T∗Tx⟩=∥x∥2. Thus ∥Ai∥,∥Bj∥≤1 (with equality on nonzero factor spaces), and the estimate proved in [F3] gives ∥A0A1∥≤∥A0∥∥A1∥≤1 (The Hilbert-space adjoint of a bounded operator, The operator norm as the least bound and as the unit-sphere or unit-ball supremum).

Proof

1.1F1F2F3F7algebra

[Self-adjointness and the square identity] By [F2] and [F3], (S⊗T)∗=S∗⊗T∗ for the four tensor products entering C, so every summand is self-adjoint and C∗=C. Expanding the square with [F2] and [F7] (all operators below are bounded, so products may be reassociated), C2=(A0⊗(B0+B1))2+(A1⊗(B0−B1))2+(A0⊗(B0+B1))(A1⊗(B0−B1))+(A1⊗(B0−B1))(A0⊗(B0+B1)); using Ai2=I and Bj2=I together with (B0+B1)2=2I+{B0,B1}, (B0−B1)2=2I−{B0,B1}, (B0+B1)(B0−B1)=[B1,B0] and (B0−B1)(B0+B1)=[B0,B1], the first two summands total IH⊗(2IK+{B0,B1})+IH⊗(2IK−{B0,B1})=4IH⊗^K, while the cross terms total −A0A1⊗[B0,B1]+A1A0⊗[B0,B1]=−[A0,A1]⊗[B0,B1]; hence C2=4I−[A0,A1]⊗[B0,B1].

1.2F1F5algebra

[The singlet correlations] On C2 with orthonormal basis (e+,e−) define Ψ:=12(e+⊗e−−e−⊗e+); then ∥Ψ∥=1 by [F1]. Direct computation with the displayed Pauli matrices gives σxe+=e−, σxe−=e+, σye+=ie−, σye−=−ie+, σze+=e+ and σze−=−e−. Hence σi⊗σiΨ=−Ψ for i=x,y,z (the two elementary tensors of the difference receive the same sign change), For clarity put q=(e+⊗e−+e−⊗e+)/2 and d±=(e+⊗e+±e−⊗e−)/2. Direct Pauli actions give (σx⊗σy)Ψ=−iq, (σy⊗σx)Ψ=iq, (σx⊗σz)Ψ=−d+, (σz⊗σx)Ψ=d+, (σy⊗σz)Ψ=id− and (σz⊗σy)Ψ=−id−. Each of q,d+,d− is orthogonal to Ψ by the elementary product-basis pairings, so every off-diagonal correlation is zero. therefore ⟨σi⊗σjΨ,Ψ⟩=−δij, and for unit vectors a,b, ⟨σa⊗σbΨ,Ψ⟩=∑i,jaibj(−δij)=−a⋅b, with σa,σb self-adjoint involutions by [F5].

2.1step 1.1F2F3F7algebra

[The norm bound] By [F3] and [F7], ∥[A0,A1]∥≤∥A0A1∥+∥A1A0∥≤2 and ∥[B0,B1]∥≤2, so [F2] gives ∥[A0,A1]⊗[B0,B1]∥≤4; hence ∥C2∥≤4+4=8 from step 1.1 and the triangle inequality. Since C is self-adjoint by step 1.1, [F3] gives ∥C∥2=∥C∗C∥=∥C2∥≤8, so ∥C∥≤22.

2.2step 1.2F5algebra

[Attainment] Take A0=σz, A1=σx, B0=σ(z+x)/2, B1=σ(z−x)/2, all self-adjoint involutions by [F5]. By step 1.2 the four expectations are E00=−z⋅z+x2=−12, E01=−z⋅z−x2=−12, E10=−x⋅z+x2=−12 and E11=−x⋅z−x2=+12, because x⋅z=0 and x⋅x=z⋅z=1. The CHSH combination of the operator C is E00+E01+E10−E11=−42=−22, so ∣⟨CΨ,Ψ⟩∣=22 and the bound (c) is attained.

3.1step 2.1F4algebra

[The quadratic-form bound] For all ψ,ψ′, [F4] and step 2.1 give ∣⟨Cψ,ψ′⟩∣≤∥C∥ ∥ψ∥ ∥ψ′∥≤22∥ψ∥∥ψ′∥; for a unit vector this is clause (c), and writing Ai⊗Bj as self-adjoint operators with the CHSH combination equal to C expresses the four-term combination as ⟨Cψ,ψ⟩.

4.1F6∎

[Local comparison] Clause (e) is exactly [F6] applied to the four ±1-valued random variables of the stated common-outcome model; no relation between that model and the operators of this theorem is asserted.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Legendre Rodrigues equation

Statement

For l∈N={0,1,2,… } define the function Pl:R→R by the Rodrigues formula (The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set)

Pl(t):=12ll!dldtl[(t2−1)l].

Then:

(a) each Pl is the polynomial function of the polynomial pl(t)=12ll!∑k: 2k≥ll(−1)l−k(lk)(2k)!(2k−l)!t2k−l ∈R[t] (so the two expressions agree at every real t, Evaluation and roots of a polynomial in a commutative target ring), pl has degree l with leading coefficient (2l)!2l(l!)2≠0 (Degree, leading coefficient and monic polynomial, with the zero polynomial having no degree);

(b) Pl(1)=1 and Pl(−1)=(−1)l;

(c) Pl satisfies Legendre's differential equation on R: (1−t2)Pl′′(t)−2tPl′(t)+l(l+1)Pl(t)=0(t∈R);

(d) P0=1, P1(t)=t and P2(t)=12(3t2−1).

No probability, measure or physical interpretation is attached to Pl; this item constructs the polynomials that later items use as angular functions.

Facts & Assumptions

Given: l∈N; the function Fl(t):=(t2−1)l; the function Pl defined by the Rodrigues formula; the polynomial pl displayed in (a).

[F2]

The general Leibniz rule: for m-times differentiable f,g, (fg)(m)=∑j=0m(mj)f(j)g(m−j) (The general Leibniz rule for the n-th derivative of a product).

[F4]

Polynomial evaluation is a ring homomorphism R[t]→R, p↦p(t), compatible with addition and multiplication, and a nonzero polynomial of degree d has leading coefficient the coefficient of td (The polynomial ring over a commutative ring as finitely supported coefficient sequences with convolution, Evaluation and roots of a polynomial in a commutative target ring, Degree, leading coefficient and monic polynomial, with the zero polynomial having no degree).

[F5]

The binomial theorem over a commutative ring: (x+y)n=∑k=0n(nk)xkyn−k in every commutative ring and for all elements x,y (The binomial theorem over an arbitrary commutative ring).

Proof

1.1givenF1F3F4F5algebra

[(a), expansion] By [F5] applied in the polynomial ring to x=t2 and y=−1, Fl(t)=∑k=0l(lk)(−1)l−kt2k as polynomials; differentiating l times termwise with [F1] and [F4] (the sum is finite, so linearity applies), every term with 2k<l differentiates to 0 and every term with 2k≥l contributes (lk)(−1)l−k(2k)!(2k−l)!t2k−l; hence Pl agrees with the polynomial pl of the Statement. All exponents 2k−l with 2k≥l are nonnegative, and the largest is 2l−l=l, attained only at k=l with coefficient (2l)!/l!≠0 by [F3]; so pl has degree l and leading coefficient (2l)!/(l!⋅2ll!)=(2l)!/(2l(l!)2), which is nonzero.

1.2givenF1F2F3algebra

[(b), endpoint values] By the factorisation Fl=(t−1)l(t+1)l and [F2], the l-th derivative at t=1 is ∑j=0l(lj)djdtj(t−1)l∣t=1⋅dl−jdtl−j(t+1)l∣t=1, and every term with j<l vanishes because djdtj(t−1)l still contains a factor (t−1); the surviving term is l!⋅(1+1)l=l!2l by [F1], [F3], so Pl(1)=l!2l2ll!=1. The same computation at t=−1 gives l!⋅(−1−1)l=(−1)ll!2l, so Pl(−1)=(−1)ll!2l2ll!=(−1)l.

1.3givenF1F2F3F4algebra

[(c), Legendre's equation] If l=0, then P0=1 and the equation is immediate. Suppose l≥1. Differentiating Fl=(t2−1)l with [F1] gives Fl′=2lt(t2−1)l−1, hence (t2−1)Fl′=2ltFl. Apply dl+1dtl+1 to both sides using [F2]; on the left the third and higher derivatives of t2−1 vanish, so (t2−1)Fl(l+2)+(l+1)(2t)Fl(l+1)+(l+12)(2)Fl(l), and on the right 2l(tFl(l+1)+(l+1)Fl(l)), where (l+12)=l(l+1)2 by [F3]; subtracting gives (t2−1)Fl(l+2)+2tFl(l+1)−l(l+1)Fl(l)=0. Since Pl=12ll!Fl(l) is exactly the l-fold derivative, dividing by the nonzero constant and multiplying by −1 yields (1−t2)Pl′′−2tPl′+l(l+1)Pl=0 at every real t: the displayed polynomial identity holds for the polynomial pl by [F4].

2.1givenF1algebra∎

[(d)] For l=0, P0=1; for l=1, P1(t)=12ddt(t2−1)=t; for l=2, P2(t)=18d2dt2(t4−2t2+1)=18(12t2−4)=3t2−12 by [F1].

The identities proved here are polynomial identities in one real variable; the item makes no angular, measure or spectral claim.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

Sphere integrals in angle coordinates

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let S2:={x∈R3:∣x∣=1} and let σ be the polar surface measure on S2 of The polar surface set function on the unit sphere, a finite Borel measure (Measures on sigma-algebras, The Borel sigma-algebra of a topological space). Let

Φ:(0,π)×(0,2π)→S2,Φ(θ,φ):=(sin⁡θcos⁡φ, sin⁡θsin⁡φ, cos⁡θ),

and let H:={(x,y,0)∈R3:x≥0} denote the closed half-plane, a Lebesgue-null Borel set (Every affine hyperplane of Rn, and hence every proper linear subspace, is Lebesgue null). Then:

(a) for every Borel function F:S2→[0,+∞], ∫S2F dσ=∫(0,π)×(0,2π)F(Φ(θ,φ)) sin⁡θ d(θ,φ), the right-hand integral being the Lebesgue integral over the open rectangle;

(b) σ(S2)=4π;

(c) the identity of (a) holds with F∈L1(S2,σ;C) in the sense that F∘Φ is Lebesgue integrable on the rectangle against sin⁡θ d(θ,φ) and the two integrals agree (real and imaginary parts, positive and negative parts).

Facts & Assumptions

Given: Countable Choice; the sphere S2; the polar surface measure σ; the parametrisation Φ; the half-plane H; the map T(r,θ,φ):=rΦ(θ,φ) on (0,∞)×(0,π)×(0,2π).

[F1]

The polar formula: for every Borel f:R3→[0,∞], ∫R3f dλ3=∫0∞r2∫S2f(rω) dσ(ω) dr, and σ is the unique Borel measure on S2 with this property (Polar coordinates decompose Lebesgue measure into r^{n-1} dr d sigma, The polar surface set function on the unit sphere).

[F2]

Change of variables for C1 diffeomorphisms of open sets: if U,V⊆Rm are open and T:U→V is a C1 bijection with C1 inverse, then for every nonnegative Borel h, ∫Vh(y) dy=∫Uh(T(x))∣det⁡DT(x)∣ dx (Borel change of variables from the compact-support formula and Radon uniqueness).

[F3]

Inverse function theorem: if T is C1 near a and DT(a) is invertible, then T restricts to a C1 bijection of a neighbourhood of a onto an open set, with C1 inverse; and C1 means total differentiability with continuous derivative (The Euclidean inverse function theorem, Continuously differentiable maps, local inverses, and local diffeomorphisms, The total (Fréchet) derivative Df(a) as the linear first-order approximation with o(∥h∥2) remainder).

[F4]

The slit K:={(x,0,z):x≥0} is closed and contained in the coordinate hyperplane {y=0}, which is Lebesgue null; complements of null sets have full measure, and λ3(K)=0 (Every affine hyperplane of Rn, and hence every proper linear subspace, is Lebesgue null, Measures on sigma-algebras).

[F5]

Tonelli's theorem for nonnegative Borel functions and Fubini's theorem for integrable functions identify multiple Lebesgue integrals with iterated integrals over σ-finite products; on a rectangle with product measure, integrals of continuous functions may be computed by iterated integrals, and Jordan–Fubini computes Riemann integrals over rectangles as iterated integrals (Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Fubini's theorem for L^1 functions on a sigma-finite product, Fubini over a bounded Jordan set when all but a content-zero family of sections are integrable).

[F6]

Sine and cosine are differentiable with sin⁡′=cos⁡, cos⁡′=−sin⁡, satisfy the addition formulas and Pythagoras' identity sin⁡2+cos⁡2=1, sin⁡>0 on (0,π), and cos⁡ is strictly decreasing on [0,π] with cos⁡0=1, cos⁡π=−1; the derivative of sine integrates by the fundamental theorem (The derivatives of sine and cosine are cosine and minus sine, The addition formulas for sine and cosine, t↦(cos⁡t,sin⁡t) is a bijection from [0,2π) onto the real unit circle, Signs, monotonicity intervals, and ranges of sine and cosine, Parity and the Pythagorean identity for sine and cosine, The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a)).

[F7]

Continuous maps pull Borel sets back to Borel sets, and indicators of Borel sets are measurable; the monotone convergence theorem computes integrals of sums of nonnegative measurable functions (A continuous map has Borel preimages of Borel sets, Monotone convergence for the integral, The Borel sigma-algebra of a topological space).

Proof

1.1F5F6F7algebra

[The angular measure τ] For Borel E⊆S2 put τ(E):=∫(0,π)×(0,2π)1E(Φ(θ,φ))sin⁡θ d(θ,φ), the integral existing in [0,∞] because the integrand is nonnegative and measurable: Φ is continuous, hence Borel measurable, and 1E∘Φ is measurable by [F7]. The set function τ is countably additive by the monotone convergence theorem applied to the partial sums of ∑k1Ek for pairwise disjoint Borel Ek, and τ(S2)=∫0π∫02πsin⁡θ dφ dθ=2π⋅2=4π by [F5] and [F6]; hence τ is a finite Borel measure.

1.2F6algebra

[T is a bijection onto the slit space] Ω:=R3∖K={y≠0}∪{y=0, x<0} is open. The image of T is contained in Ω: for φ∈(0,π) the second coordinate rsin⁡θsin⁡φ is positive, for φ∈(π,2π) it is negative, and for φ=π one has y=0 with first coordinate rsin⁡θcos⁡π=−rsin⁡θ<0. Conversely, let x∈Ω and put r:=∣x∣>0 and u:=x/r∈S2. Since u∉K, the case u=(0,0,±1) is excluded (it has y=0, x-coordinate 0≥0), so ∣u3∣<1; by [F6] there is a unique θ∈(0,π) with cos⁡θ=u3, and then sin⁡θ>0 and v:=(u1,u2)/sin⁡θ satisfies v12+v22=1 by Pythagoras [F6]. By the parametrisation of the unit circle [F6] there is a unique φ∈[0,2π) with v=(cos⁡φ,sin⁡φ); φ=0 would give u=(sin⁡θ,0,u3) with y=0 and x-coordinate sin⁡θ>0, i.e. u∈K, which is excluded, so φ∈(0,2π). Then x=T(r,θ,φ) with (r,θ,φ) uniquely determined, so T is a bijection of (0,∞)×(0,π)×(0,2π) onto Ω.

1.3F6algebra

[T is C1 with positive Jacobian] The components of T are products and composites of the differentiable functions sin⁡,cos⁡ and the coordinate functions by [F6], so T is C1 with total derivative the matrix of partial derivatives (The total (Fréchet) derivative Df(a) as the linear first-order approximation with o(∥h∥2) remainder, Total differentiability gives a local O(∥h∥2) increment bound and therefore continuity), namely the columns ∂rT=(sin⁡θcos⁡φ,sin⁡θsin⁡φ,cos⁡θ), ∂θT=(rcos⁡θcos⁡φ,rcos⁡θsin⁡φ,−rsin⁡θ) and ∂φT=(−rsin⁡θsin⁡φ,rsin⁡θcos⁡φ,0); expanding the determinant along the third column and using cos⁡2φ+sin⁡2φ=1 and sin⁡2θ+cos⁡2θ=1 [F6] gives det⁡DT=r2sin⁡θ(cos⁡2θ+sin⁡2θ)=r2sin⁡θ>0 (The Jacobian determinant of a square-dimensional C1 map is the determinant of its Jacobian matrix).

2.1step 1.2step 1.3F2F3

[T is a C1 diffeomorphism onto Ω] At each point of the domain DT is invertible by step 1.3, so by the inverse function theorem [F3] T restricts to a C1 bijection of a neighbourhood of that point onto an open subset of Ω with C1 inverse. Since T is injective with image Ω by step 1.2, these local inverses agree on overlaps with the unique set-theoretic inverse T−1:Ω→(0,∞)×(0,π)×(0,2π), which is therefore C1: near any point of Ω it coincides with a local C1 inverse. Hence T is a C1 diffeomorphism of open sets, and [F2] applies to it.

3.1step 2.1step 1.1F2F4F5F7algebra

[The polar formula for τ] Let f:R3→[0,∞] be Borel. By [F4], λ3(K)=0 and ∫R3f=∫Ωf. Applying the change of variables formula [F2] to the diffeomorphism T of step 2.1 and the nonnegative Borel function f∣Ω gives ∫Ωf=∫(0,∞)×(0,π)×(0,2π)f(T(r,θ,φ)) r2sin⁡θ d(r,θ,φ). By Tonelli's theorem [F5] the triple integral factors as ∫0∞r2(∫(0,π)×(0,2π)f(rΦ(θ,φ))sin⁡θ d(θ,φ))dr=∫0∞r2∫S2f(rω) dτ(ω) dr: for each fixed r the function ω↦f(rω) is nonnegative and Borel, and its τ-integral equals the rectangle integral of f(rΦ(⋅))sin⁡(⋅) by the defining formula for τ, extended from indicators to nonnegative Borel integrands by the monotone convergence theorem [F7]. Hence τ satisfies the polar formula of [F1].

4.1step 1.1step 3.1F1F5F7algebra∎

[Identification and consequences] By the uniqueness clause of [F1] the measure satisfying the polar formula is σ; hence τ=σ, which is exactly the identity of (a) for indicators and, by the monotone convergence theorem, for all nonnegative Borel F [F7]; clause (b) is the computation τ(S2)=4π of step 1.1. For (c), if F∈L1(S2,σ;C) then applying (a) to ∣F∣ shows F∘Φ is integrable against sin⁡θ d(θ,φ), because the integral of ∣F∣ over S2 is finite; splitting F into real and imaginary parts and each into positive and negative parts, and applying (a) to the four nonnegative functions (or Fubini, [F5]) gives the identity for F.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Associated Legendre function normalization

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let Pl be the Legendre polynomials of Legendre Rodrigues equation, and for l∈N and 0≤m≤l define the function wlm:(−1,1)→R by

wlm(t):=(−1)m(1−t2)m/2 dmdtmPl(t),

where (1−t2)1/2>0 denotes the nonnegative square root. Let σ be the polar surface measure on S2 and let the angle coordinates be Φ(θ,φ)=(sin⁡θcos⁡φ,sin⁡θsin⁡φ,cos⁡θ) (Sphere integrals in angle coordinates, The polar surface set function on the unit sphere). The image of Φ is S2∖C, where C:={(x,0,z)∈S2:x≥0}; C is closed and σ(C)=0 by applying the integral identity of Sphere integrals in angle coordinates to its indicator. For l≥m≥0, define the representative Y~lm:S2→C on S2∖C by Y~lm(Φ(θ,φ)):=blmwlm(cos⁡θ)eimφ2π, and on the cut C by Y~lm(sin⁡θ,0,cos⁡θ):=blmwlm(cos⁡θ)/2π(0≤θ≤π). The limits at φ=0 and 2π agree because m is an integer, and the endpoint values of wlm from (a) give the limits at the poles; hence Y~lm is continuous. For 1≤m≤l set Y~l,−m:=(−1)mY~lm‾, and write Ylm for the class of Y~lm in L2(S2,σ;C). Then:

(a) each wlm extends continuously to [−1,1], with wlm(±1)=0 when m≥1, and on (−1,1) it satisfies the associated Legendre equation −((1−t2)wlm′)′+m21−t2wlm=l(l+1)wlm(−1<t<1), equivalently (1−t2)wlm′′−2twlm′+[l(l+1)−m21−t2]wlm=0;

(b) for fixed m≥0 the functions wlm, l≥m, are pairwise orthogonal in L2((−1,1),dt): ∫−11wlmwkm dt=0 for k≠l; moreover 0<∫−11wlm2 dt<+∞;

(c) with blm:=(∫−11wlm2 dt)−1/2>0 and the representatives just defined, one has ∣Yl,−m∣=∣Ylm∣ and Ylm‾=(−1)mYl,−m for all l and all m with ∣m∣≤l;

(d) the family {Ylm:l≥0, −l≤m≤l} is orthonormal in L2(S2,σ;C) (Real and complex inner product spaces, with the inner product linear in the first argument, The space Lp(μ) as the quotient by null functions): for all (l,m)≠(l′,m′) with ∣m∣≤l, ∣m′∣≤l′, the integral ∫S2YlmYl′m′‾ dσ equals 0, and it equals 1 when (l,m)=(l′,m′);

(e) wl0=Pl, and wll(t)=(−1)l(2l)!2ll!(1−t2)l/2.

The functions Ylm are complex-valued L2-functions on the sphere; no physical interpretation as angular momentum eigenfunctions or spherical harmonics of an atom is asserted.

Facts & Assumptions

Given: the Legendre polynomials Pl; integers l≥0, 0≤m≤l; the functions wlm; the measure σ and the parametrisation Φ; the constants blm and functions Ylm.

[F1]

Pl satisfies (1−t2)Pl′′−2tPl′+l(l+1)Pl=0 on R, is a polynomial of degree l, and Pl(m) is therefore the polynomial function of a degree-(l−m) polynomial (Legendre Rodrigues equation).

[F3]

Sphere integrals in angle coordinates: for Borel F≥0, ∫S2F dσ=∫(0,π)×(0,2π)F(Φ(θ,φ))sin⁡θ d(θ,φ), and the complex L1 version; σ(S2)=4π (Sphere integrals in angle coordinates).

[F4]

Characters of the circle: ∫02πeikφ dφ=2π if k=0 and =0 for every nonzero k∈Z (The trigonometric characters are orthonormal in L2 of the torus, The complex exponential by its power series).

[F5]

A continuous real function on the compact interval [−1,1] is (Riemann and Lebesgue) integrable; products of continuous functions are continuous, and ∫−11g>0 for a continuous g≥0 that is not identically zero (A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion, Complex Lp classes and Euclidean test-function conventions).

[F6]

If G:[a,b]→R is differentiable at every point of [a,b] as a function on [a,b] (with one-sided derivatives at a,b) and G′ is integrable on [a,b], then ∫abG′=G(b)−G(a) (The second fundamental theorem: if G is differentiable on [a,b] with G′=f and f is integrable, then ∫abf=G(b)−G(a)).

[F7]

For f,g∈L2 the product fg‾ is in L1 with ∫∣fg‾∣≤∥f∥2∥g∥2; the inner product in L2(S2,σ;C) is linear in the first argument; conjugation satisfies zw‾=z‾ w‾ and ∣z∣2=zz‾ (Cauchy-Schwarz inequality for L2, Real and complex inner product spaces, with the inner product linear in the first argument, Conjugation is an involutive real-field automorphism, zz‾=∣z∣2, and modulus is definite, multiplicative, and subadditive, Complex Lp classes and Euclidean test-function conventions).

[F8]

For the given 0≤m≤l, dmPl/dtm is a polynomial of degree l−m with nonzero leading coefficient, so wlm is not the zero function on (−1,1), because (1−t2)m/2>0 there (Legendre Rodrigues equation, The factorial n! and the falling factorial nk‾, defined by recursion in N).

[F9]

If G:[a,b]→R is continuous on [a,b] and differentiable at every point of (a,b) as a function on [a,b], the mean value theorem gives a c∈(a,b) with G(b)−G(a)=G′(c)(b−a) (The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a)).

Proof

1.1F1F2algebra

[Differentiated equation] Let v:=Pl(m). Differentiating the Legendre equation of [F1] m times with the general Leibniz rule [F2] gives, since the third and higher derivatives of 1−t2 vanish, (1−t2)v′′−2t(m+1)v′+[l(l+1)−m(m+1)]v=0: the coefficient of v′′ is 1−t2; that of v′ collects −2t from differentiating 1−t2 once times the derivative count (m1) and −2t from the term −2tPl′, total −2t(m+1); and that of v collects (m2)(−2)+m(−2)+l(l+1)=−m(m−1)−2m+l(l+1)=l(l+1)−m(m+1).

1.2F1F2algebra

[Continuity at the endpoints] For m=0, wl0=Pl is a polynomial; for m≥1, wlm=(−1)m(1−t2)m/2Q(t) with Q=Pl(m) a polynomial, so ∣wlm(t)∣≤(1−t2)m/2sup⁡[−1,1]∣Q∣→0 as t→±1, and the function extended by wlm(±1):=0 is continuous on [−1,1].

2.1step 1.1F2algebra

[Substitution] Put f(t):=(1−t2)m/2>0 on (−1,1), so that wlm=(−1)mfv; the sign is constant and will be carried along. Put u:=(−1)mv, so w=fu and u satisfies the same homogeneous equation as v in step 1.1. From u=f−1w and [F2], u′=f−1w′−f−2f′w and u′′=f−1w′′−2f−2f′w′+(2f−3(f′)2−f−2f′′)w. Substituting into step 1.1 and multiplying by f gives (1−t2)w′′+[−2f′f(1−t2)−2t(m+1)]w′+[−f′′f(1−t2)+2(f′f)2(1−t2)+2t(m+1)f′f−m(m+1)+l(l+1)]w=0. Direct differentiation gives f′f=−mt1−t2 and therefore f′′f=(f′f)′+(f′f)2=−m−mt2(1−t2)2+m2t2(1−t2)2=m(−1+(m−1)t2)(1−t2)2. Hence the coefficient of w′ is 2mt−2t(m+1)=−2t, and the coefficient of w equals m−m(m+1)t21−t2−m(m+1)+l(l+1)=l(l+1)−m21−t2, using m−m(m+1)t21−t2−m(m+1)=−m21−t2 (expand the left side over the common denominator 1−t2). This is the associated equation.

3.1step 2.1step 1.2F2F5F6F9algebra

[Orthogonality] Fix m and let k≠l, k≥m. Multiply the associated equation for wlm by wkm and that for wkm by wlm, subtract, and use −((1−t2)wlm′)′wkm+((1−t2)wkm′)′wlm=−ddt[(1−t2)(wlm′wkm−wlmwkm′)]: (l(l+1)−k(k+1))wlmwkm=−ddtG(t) with G(t):=(1−t2)(wlm′(t)wkm(t)−wlm(t)wkm′(t)) on (−1,1). By step 1.2, put Ql=Pl(m) and Qk=Pk(m). If m=0, then both w functions are polynomials and the factor 1−t2 gives G(t)→0 at both endpoints. If m≥1, put f(t)=(1−t2)m/2, so f′=−mt(1−t2)m/2−1 and wlm=(−1)mfQl, wkm=(−1)mfQk; since Ql,Qk and their derivatives are bounded on [−1,1], each term in (1−t2)wlm′wkm and (1−t2)wlmwkm′, after the product rule, is bounded by a constant times (1−t2)m or (1−t2)m+1, so again G(t)→0 at both endpoints. Extend G by G(±1)=0; these limits make it continuous on [−1,1]. The subtracted equations give G′(t)=q(t):=−(l(l+1)−k(k+1))wlm(t)wkm(t) on (−1,1), and step 1.2 makes q continuous on [−1,1]. For each x∈(−1,1), [F9] on [x,1] gives (G(1)−G(x))/(1−x)=q(ξx) for some ξx∈(x,1), so as x↑1 this quotient tends to q(1). On [−1,x], [F9] gives (G(x)−G(−1))/(x+1)=q(ηx) for some ηx∈(−1,x); as x↓−1 this quotient tends to q(−1). Thus G is differentiable at every point of [−1,1] with the endpoint convention of [F6], and G′=q is continuous, hence integrable by [F5]. Now [F6] gives ∫−11G′=G(1)−G(−1)=0, so (l(l+1)−k(k+1))∫−11wlmwkm dt=0 and, since l≠k, ∫−11wlmwkm dt=0.

3.2step 2.1step 1.2F5F8algebra

[Positivity and finiteness of the norms] By steps 2.1 and 1.2, wlm2 is continuous on [−1,1], so its integral is finite and nonnegative by [F5]; it is not zero because wlm is not the zero function on (−1,1): Pl(m) is a nonzero polynomial by [F8] and (1−t2)m/2>0 in the interior, so wlm is nonzero at every point where the polynomial does not vanish, and a continuous nonzero function has positive integral by [F5]. Hence 0<∫−11wlm2 dt<+∞ and blm of the Statement is a well-defined positive real.

3.3step 2.1F1F7algebra

[The functions Ylm and their conjugation] The formula in the Statement defines continuous representatives on S2, so each class Ylm belongs to L2(S2,σ) because σ is finite; their chart moduli satisfy ∣Ylm(Φ(θ,φ))∣=blm∣wlm(cos⁡θ)∣/2π. By definition Yl,−m=(−1)mYlm‾, so ∣Yl,−m∣=∣Ylm∣ and, conjugating twice, Yl,−m‾=(−1)mYlm, which is the displayed identity for every m≥0. This establishes (c), and the identities in (e) follow from Pl(0)=Pl and Pl(l)=(2l)!2ll! (the leading coefficient of Pl times l!, by [F1]).

4.1step 3.1step 3.2step 3.3F3F4F7algebra∎

[Orthonormality] Let m,m′≥0. The extensions on C do not affect integrals because σ(C)=0. Thus [F3] applies to the L1 function YlmYl′m′‾, whose integral against σ is the inner product; with [F4], this gives ⟨Ylm,Yl′m′⟩=blmbl′m′12π(∫02πei(m−m′)φdφ)(∫−11wlm(t)wl′m′(t) dt)=δmm′ blmbl′m∫−11wlmwl′m dt, which by steps 3.1 and 3.2 equals 0 for l≠l′ and 1 for l=l′. For entries with negative lower index, write the inner product as a product of conjugation operations using [F7]: ⟨Yl,−m,Yl′,−m′⟩=(−1)m+m′⟨Ylm‾,Yl′m′‾⟩=(−1)m+m′⟨Ylm,Yl′m′⟩‾, and ⟨Ylm,Yl′,−m′⟩=(−1)m′∫S2YlmYl′m′ dσ, whose rectangle form has φ-factor ∫02πei(m+m′)φdφ, vanishing unless m=m′=0 by [F4], in which case it reduces to the already treated case; collecting the cases gives orthonormality for all pairs. This proves (d), and (a), (b) were proved in steps 2.1, 3.1 and 3.2.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Weighted angular completeness

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let Pl be the Legendre polynomials of Legendre Rodrigues equation and wlm the associated functions of Associated Legendre function normalization, and let L2=L2((−1,1),dt;C) with inner product linear in the first argument (The space Lp(μ) as the quotient by null functions, Complex Lp classes and Euclidean test-function conventions). Then, for each fixed m∈N:

(a) for every l≥m the polynomial Pl(m) has degree l−m and nonzero leading coefficient, and for every d∈N the polynomials Pm(m),Pm+1(m),…,Pm+d(m) form a basis of the space of polynomial functions of degree at most d on [−1,1];

(b) the family (wlm)l≥m has dense linear span in L2: if f∈L2 satisfies ∫−11f(t)wlm(t) dt=0 for every l≥m, then f=0 almost everywhere; consequently the orthogonal family (wlm)l≥m is complete in L2 and its normalization (wlm/∥wlm∥2)l≥m is a Hilbert basis of L2, so every f∈L2 has f=∑l≥m(∫−11f wlm dt)∥wlm∥2−2wlm with the finite-subset-net convergence of Orthonormal families, complete orthonormal systems and Hilbert bases;

(c) for every fixed m the orthogonal projection of L2 onto the one-dimensional span of wmm is f↦∥wmm∥2−2(∫−11f wmm dt)wmm, and wmm(t)=(−1)m(2m)!2mm!(1−t2)m/2.

Facts & Assumptions

Given: m∈N; the Legendre polynomials Pl; the functions wlm; the complex Hilbert space L2((−1,1),dt).

[F1]

Pl is a polynomial of degree l with nonzero leading coefficient al=(2l)!2l(l!)2 (Legendre Rodrigues equation). Iterating the monomial derivative rule (For a natural n≥1 the function x↦xn is differentiable everywhere with derivative ι(n) x n−1; for n=0 it is the constant 1, with derivative 0; for a natural n≥1 the function x↦x−n is differentiable at every x≠0 with derivative −ι(n) x−n−1; consequently every polynomial function is differentiable at every real, with the derivative computed term by term) and factorial recursion (The factorial n! and the falling factorial nk‾, defined by recursion in N) gives, for l≥m, leading term all!/(l−m)! tl−m in Pl(m); lower-degree terms either vanish or have smaller degree. Thus its degree is l−m, and Pm(m)=amm!.

[F2]

wlm=(−1)m(1−t2)m/2Pl(m) is continuous on [−1,1], the family (wlm)l≥m is orthogonal in L2, and ∥wlm∥22=∫−11wlm2 dt∈(0,∞) (Associated Legendre function normalization).

[F3]

The real polynomial functions are uniformly dense in the real continuous functions on [−1,1]: for every continuous g:[−1,1]→R and every ε>0 there is a polynomial p with sup⁡[−1,1]∣g−p∣<ε (Real Stone–Weierstrass theorem for compact Hausdorff spaces).

[F4]

Under Countable Choice, for every bounded interval (a,b) the complex continuous functions on [a,b] are dense in complex Lp((a,b)) for 1≤p<+∞, in particular in L2 (Continuous functions are dense in Lp of finite tori and of bounded intervals (2)). For continuous g, integrating ∣g(t)∣2≤∥g∥∞2 over (a,b) gives ∥g∥22≤(b−a)∥g∥∞2, hence ∥g∥2≤b−a∥g∥∞ by the defining L2 norm (The space Lp(μ) as the quotient by null functions, Complex Lp classes and Euclidean test-function conventions).

[F5]

∥fg∥1≤∥f∥2∥g∥2 and ∣∫fg‾∣≤∥f∥2∥g∥2 in L2 (Cauchy-Schwarz inequality for L2, Complex Lp classes and Euclidean test-function conventions).

[F6]

For a family that is orthonormal and complete in a Hilbert space, the expansion, Parseval identity and coefficient formulas hold with finite-subset-net convergence (Orthonormal families, complete orthonormal systems and Hilbert bases, Parseval equivalences for an orthonormal family).

Proof

1.1F1algebra

[(a)] By [F1], Pl(m) has degree l−m and nonzero leading coefficient for each l≥m. Hence in the ordered list Pm(m),Pm+1(m),…,Pm+d(m) the j-th polynomial has degree j, 0≤j≤d: the leading-degree matrix is triangular with nonzero diagonal, so the polynomials are linearly independent; they are d+1 vectors in the (d+1)-dimensional space of polynomials of degree at most d, hence a basis.

2.1step 1.1F2algebra

[(b), first reduction] Let f∈L2 satisfy ∫−11fwlm dt=0 for all l≥m, and put g(t):=f(t)(1−t2)m/2. Since t↦(1−t2)m/2 is continuous and bounded on [−1,1], g∈L2 and ∥g∥2≤∥f∥2; moreover for every l≥m, ∫−11g Pl(m) dt=(−1)m∫−11fwlm dt=0 by [F2]. By step 1.1, every polynomial q of degree at most d is a linear combination of Pm(m),…,Pm+d(m), so ∫−11gq dt=0 for every polynomial q.

3.1step 2.1F3F4F5algebra

[(b), density step] Polynomials are dense in L2((−1,1)): given φ∈L2 and ε>0, [F4] supplies a continuous h on [−1,1] with ∥φ−h∥2<ε/2, and apply [F3] separately to the real and imaginary parts of h, obtaining real polynomials p1,p2 whose respective uniform errors are below ε/(42), and put p=p1+ip2; then sup⁡[−1,1]∣h−p∣<ε/(22), so ∥h−p∥2<ε/2 by [F4] (applied to the interval length 2), whence ∥φ−p∥2<ε. Applying this to the complex-valued g componentwise (real and imaginary parts) gives polynomials qj with ∥g−qj∥2→0; since ∫gq=0 for all polynomials, Cauchy–Schwarz [F5] gives ∥g∥22=∫gg‾=lim⁡j∫gqj‾=0, so g=0 almost everywhere. On (−1,1) the factor (1−t2)m/2 is strictly positive, so f=g/(1−t2)m/2=0 almost everywhere.

4.1step 3.1F2F6algebra

[(b), completeness and expansion] By [F2] the family (wlm)l≥m is orthogonal with all norms nonzero, so (wlm/∥wlm∥2)l≥m is orthonormal; if a vector of L2 is orthogonal to the whole family, it is orthogonal to every wlm and vanishes by step 3.1, so the orthonormal family is complete in the sense of [F6]; the expansion, normalization and coefficient formulas follow from [F6].

5.1step 4.1F1F2F5algebra∎

[(c)] The orthogonal projection onto the span of a unit vector u is x↦⟨x,u⟩u, and for an unnormalized nonzero vector w it is x↦∥w∥2−2⟨x,w⟩w; with w=wmm and the inner product linear in the first argument this is the displayed formula. Indeed its value lies in the span of w, and its residual has zero inner product with w, uniquely determining the orthogonal projection; Cauchy–Schwarz bounds its norm by ∥x∥2. The explicit value of wmm follows from Pm(m)=(2m)!2mm! (the leading coefficient of Pm times m!) inserted into the definition of wlm.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Sphere ONB

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let σ be the polar surface measure on S2 (The polar surface set function on the unit sphere) and let (Ylm)l≥0, −l≤m≤l be the functions of Associated Legendre function normalization (The complex exponential by its power series). Then (Ylm) is a complete orthonormal family of the complex Hilbert space L2(S2,σ) (Orthonormal families, complete orthonormal systems and Hilbert bases, The space Lp(μ) as the quotient by null functions): the family is orthonormal, and if F∈L2(S2,σ) satisfies ∫S2FYlm‾ dσ=0 for all l≥0 and −l≤m≤l, then F=0 almost everywhere. Consequently every F has the expansion F=∑l=0∞∑m=−ll⟨F,Ylm⟩Ylm,∥F∥2=∑l,m∣⟨F,Ylm⟩∣2, with ⟨F,Ylm⟩=∫S2FYlm‾ dσ and finite-subset-net convergence (Parseval equivalences for an orthonormal family).

Facts & Assumptions

Given: the measure σ on S2; the functions Ylm=blmwl∣m∣(cos⁡θ)eimφ/2π with the convention Yl,−m=(−1)mYlm‾; a function F∈L2(S2,σ).

[F1]

Orthonormality of the family (Ylm) in L2(S2,σ) (Associated Legendre function normalization (d)).

[F2]

Angle coordinates: ∫S2G dσ=∫(0,π)×(0,2π)G(Φ(θ,φ))sin⁡θ d(θ,φ) for G∈L1, with Φ(θ,φ)=(sin⁡θcos⁡φ,sin⁡θsin⁡φ,cos⁡θ) (Sphere integrals in angle coordinates).

[F3]

Weighted completeness: for each k∈N, if h∈L2((−1,1),dt) satisfies ∫−11h(t)wlk(t) dt=0 for all l≥k, then h=0 a.e. (Weighted angular completeness (b)).

[F4]

The normalized characters ek(φ)=eikφ/2π, k∈Z, are orthonormal in L2((0,2π),dφ), and their linear span is dense: trigonometric polynomials are uniformly dense in the continuous periodic functions (Trigonometric polynomials are uniformly dense in continuous functions on the torus), and continuous functions on [0,2π] with equal endpoint values are dense in L2((0,2π)) (Continuous functions are dense in Lp of finite tori and of bounded intervals, The trigonometric characters are orthonormal in L2 of the torus); hence if h∈L2((0,2π)) has ∫02πh(φ)e−ikφ dφ=0 for every k∈Z, then h=0 a.e.

[F5]

Tonelli's and Fubini's theorems identify the rectangle integral in [F2] with iterated integrals, for nonnegative and for L1 integrands (Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Fubini's theorem for L^1 functions on a sigma-finite product).

[F6]

If ∫S2∣F∣2 dσ<∞, then by [F2] and [F5] the iterated integral ∫0π∫02π∣F(Φ(θ,φ))∣2sin⁡θ dφ dθ is finite, so for almost every θ the function φ↦F(Φ(θ,φ)) lies in L2((0,2π)) (The space Lp(μ) as the quotient by null functions, Complex Lp classes and Euclidean test-function conventions).

[F7]

∣∫gh‾∣≤∥g∥2∥h∥2 and vanishing of all pairings against a dense family forces vanishing (Cauchy-Schwarz inequality for L2).

Proof

1.1givenF2F3F5F6F7algebra

[Completeness: Fourier coefficients along the circles] Suppose F∈L2(S2,σ) is orthogonal to every Ylm, and for k∈Z and θ∈(0,π) define Gk(θ):=12π∫02πF(Φ(θ,φ))e−ikφ dφ; by [F6] and Cauchy–Schwarz [F7] each Gk is defined for almost every θ, and ∫0π∣Gk(θ)∣2sin⁡θ dθ≤∫0π∫02π∣F(Φ(θ,φ))∣2 dφ sin⁡θ dθ<∞, so Gk∈L2(sin⁡θ dθ) and, after the reparametrisation t=cos⁡θ, Gk∈L2((−1,1),dt). For l≥k≥0, [F2] and [F5] compute ⟨F,Ylk⟩=blk∫0πGk(θ)wlk(cos⁡θ)sin⁡θ dθ=blk∫−11Gk(t)wlk(t) dt; hence ∫−11Gk(t)wlk(t) dt=0 for every l≥k, and by [F3] Gk=0 almost everywhere for every k≥0. For k<0 write m:=−k>0; the convention gives Yl,−m‾=(−1)mYlm, so orthogonality of F to Yl,−m gives 0=(−1)m∫S2F Ylm dσ, whose rectangle form is blm∫0πG−m(θ)wlm(cos⁡θ)sin⁡θ dθ with 0≤m≤l; by [F3] again, G−m=Gk=0 almost everywhere. So all circle Fourier coefficients of F vanish for almost every θ.

2.1step 1.1F2F4F6algebra

[Completeness: vanishing on the spheres] By [F4] applied to the L2((0,2π))-function φ↦F(Φ(θ,φ)), which exists for almost every θ by [F6] and has all Fourier coefficients zero by step 1.1, one has F(Φ(θ,φ))=0 for almost every φ and almost every θ; hence F∘Φ=0 almost everywhere on the rectangle, and by [F2] (which for the nonnegative function ∣F∣ identifies ∫S2∣F∣ dσ with the rectangle integral) F=0 almost everywhere on S2.

3.1step 2.1F1F7algebra∎

[The expansion] By [F1] the family is orthonormal, and by steps 1.1 and 2.1 it is complete; the expansion, Parseval identity and coefficient formula are the Hilbert-basis statements of Orthonormal families, complete orthonormal systems and Hilbert bases and Parseval equivalences for an orthonormal family, with finite-subset-net convergence.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Sphere Laplacian diagonal domain

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Fix a real parameter ℏ>0. Let (Ylm)l≥0,−l≤m≤l be the complete orthonormal family of L2(S2,σ) constructed in Associated Legendre function normalization and Sphere ONB, and let H:=span⁡{Ylm:−l≤m≤l, l∈N} be its algebraic (finite) span. Define operators A and Lz on L2(S2,σ) by the diagonal lemma Diagonal adjoint and graph-core lemma:

D(A):={F:∑l,ml2(l+1)2 ∣⟨F,Ylm⟩∣2<∞},AF:=∑l,ml(l+1)⟨F,Ylm⟩Ylm, D(Lz):={F:∑l,mm2 ∣⟨F,Ylm⟩∣2<∞},LzF:=∑l,mℏm ⟨F,Ylm⟩Ylm,

the sums being in the Hilbert-space norm (Hilbert space, Orthonormal families, complete orthonormal systems and Hilbert bases). Finally, for F in the finite span H, use the canonical continuous representatives Y~lm from Associated Legendre function normalization: set F~:=∑l,m⟨F,Ylm⟩Y~lm, a finite sum uniquely determined by the orthonormal coefficients. The same convention applies to AF,LzF∈H. On the chart Φ(θ,φ)=(sin⁡θcos⁡φ,sin⁡θsin⁡φ,cos⁡θ) with (θ,φ)∈(0,π)×(0,2π), write F(θ,φ):=F~(Φ(θ,φ)) and put

(SF)(θ,φ):=−1sin⁡θ∂θ(sin⁡θ ∂θF)(θ,φ)−1sin⁡2θ∂φ2F(θ,φ),(DF)(θ,φ):=−iℏ ∂φF(θ,φ),

where the partial derivatives are those of Directional derivatives and partial derivatives of a map U⊆Rm→Rn. Then:

(a) A and Lz are self-adjoint operators on L2(S2,σ) with the displayed domains; on Ylm they act by AYlm=l(l+1)Ylm and LzYlm=ℏmYlm (Symmetric, self-adjoint and essentially self-adjoint operators); A is nonnegative, ⟨AF,F⟩≥0, and A1=0;

(b) H is a core of A and of Lz: the finite truncations FN:=∑l≤N∑m⟨F,Ylm⟩Ylm of an element F of D(A) satisfy FN→F and AFN→AF in L2, and analogously with Lz;

(c) for every F∈H and every chart point with sin⁡θ≠0, SF(θ,φ)=(AF)(θ,φ),DF(θ,φ)=(LzF)(θ,φ);

(d) for every l the orthogonal projection Pl of L2(S2,σ) onto the (2l+1)-dimensional span of Yl,−l,…,Yll, defined by PlF:=∑m=−ll⟨F,Ylm⟩Ylm (the unique vector in that span with orthogonal residual), reduces both A and Lz, i.e. PlA⊆APl and PlLz⊆LzPl; consequently A and Lz commute on the finite span H, and for every F∈D(A) with ∑l,m(l(l+1)ℏm)2∣⟨F,Ylm⟩∣2<∞ one has ALzF=LzAF=∑l,ml(l+1)ℏm⟨F,Ylm⟩Ylm. Writing clm=⟨F,Ylm⟩ and λl=l(l+1), the maximal product domains are D(ALz)={F∈L2:∑l,m(λlℏm)2∣clm∣2<∞} and D(LzA)=D(A)∩D(ALz). Thus this is exactly the common domain D(ALz)∩D(LzA)=D(LzA), while the two maximal product domains need not coincide. The product-weight condition alone does not imply F∈D(A): the vector with coefficients cl0=1/l for l≥1 and all other coefficients zero belongs to L2, has zero product weight, but does not belong to D(A). The product-weight condition is not automatic for all F∈D(A)∩D(Lz).

The differential expression S is evaluated only away from the poles, where the coordinate formulas are defined; no claim is made about a maximal differential operator on all of S2.

Facts & Assumptions

Given: the complete orthonormal family (Ylm); the finite span H; the diagonal operators A,Lz; the differential expressions S,D; the canonical representatives Y~lm=blmwlm(cos⁡θ)eimφ/2π for m≥0, with Y~l,−m=(−1)mY~lm‾ for m>0, and their finite coefficient sums F~.

[F1]

Diagonal lemma: for a complete orthonormal family (ej) and real eigenvalues λj the diagonal operator S with Sej=λjej and D(S)={x:∑jλj2∣⟨x,ej⟩∣2<∞} is self-adjoint and is the unique self-adjoint extension of its finite span; the finite span is a core, the finite truncations converge in the graph norm, and ⟨Sx,x⟩=∑jλj∣⟨x,ej⟩∣2 (Diagonal adjoint and graph-core lemma).

[F2]

(Ylm) is a complete orthonormal family of L2(S2,σ) (Sphere ONB).

[F3]

wlm is continuous on [−1,1], differentiable on (−1,1), and satisfies the associated Legendre equation −((1−t2)wlm′)′+m21−t2wlm=l(l+1)wlm (Associated Legendre function normalization (a)).

[F4]

Y~lm(Φ(θ,φ))=blmwlm(cos⁡θ)eimφ/2π for m≥0 and Y~l,−m=(−1)mY~lm‾; their L2 classes are Ylm (Associated Legendre function normalization (c)).

[F5]

Sine and cosine are differentiable with sin⁡′=cos⁡, cos⁡′=−sin⁡; partial derivatives are computed by the usual rules, including the chain rule for t=cos⁡θ and the product rule (The derivatives of sine and cosine are cosine and minus sine, The chain rule, in one line from Carathéodory: if g is differentiable at c and f is differentiable at g(c), then f∘g is differentiable at c with (f∘g)′(c)=f′(g(c)) g′(c), Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, Directional derivatives and partial derivatives of a map U⊆Rm→Rn). The complex exponential has complex derivative itself (The complex exponential is entire and its complex derivative is itself). For m≠0, its defining difference quotient along increments imh gives ddφeimφ=imeimφ; for m=0 this is the constant derivative 0. Repeating gives d2dφ2eimφ=−m2eimφ.

[F6]

For a finite orthonormal family (uj) in a Hilbert space, set Px:=∑j⟨x,uj⟩uj. Its value lies in their span and ⟨x−Px,uj⟩=0 for each j, so it is uniquely determined by orthogonal residual. The finite pairings give P2=P and ⟨Px,y⟩=⟨x,Py⟩; also ∥x∥2=∥Px∥2+∥x−Px∥2, so ∥P∥≤1. Thus it is the bounded self-adjoint orthogonal projection onto that span; no finite-dimensional ambient hypothesis is needed (Real and complex inner product spaces, with the inner product linear in the first argument, The induced length is a norm).

Proof

1.1F1F2F4algebra

[(a) and (b)] Enumerate the indices by j=l2+l+m, which bijects the block −l≤m≤l to l2≤j<(l+1)2 and all blocks to N. Apply [F1] to this enumeration of the complete orthonormal family (Ylm) of [F2] with the real eigenvalues λlm=l(l+1): the resulting operator is exactly A as defined, is self-adjoint, has the displayed domain, satisfies ⟨AF,F⟩=∑lml(l+1)∣⟨F,Ylm⟩∣2≥0, and has the finite span H as a core with the stated truncation property; the eigenvalues are nonnegative, so A1=0 because 1=4πY00 by [F4] (with P0=1, b00=1/2), whose eigenvalue is zero. The same application with the real eigenvalues ℏm, m∈Z (whose square summability domain is the displayed one because ℏ>0), gives clause (a) and the second half of (b) for Lz.

1.2F3F4F5algebra

[(c) for A] For F=Ylm with m≥0, evaluate its canonical representative; by [F4], F(θ,φ)=c w(cos⁡θ)eimφ with w=wlm and c=blm/2π. Computing the partial derivatives with [F5]: ∂φ2F=−m2F and, writing t=cos⁡θ, ∂θF=−csin⁡θ w′(t)eimφ, so sin⁡θ ∂θF=−csin⁡2θ w′(t)eimφ and ∂θ(sin⁡θ ∂θF)=c(−2sin⁡θcos⁡θ w′+sin⁡3θ w′′)eimφ; dividing by −sin⁡θ gives c(2tw′−(1−t2)w′′)eimφ=−c ((1−t2)w′)′eimφ, where in the last step ((1−t2)w′)′=−2tw′+(1−t2)w′′. Also −sin⁡−2θ ∂φ2F=+m2(1−t2)−1F. Adding, SF=c[−((1−t2)w′)′+m21−t2w]eimφ, which by [F3] equals c l(l+1)w eimφ=l(l+1)F=AF. For negative indices, Yl,−m=(−1)mYlm‾ and the real-coefficient differential expression S commutes with conjugation, so the same eigenvalue identity holds. Both sides are linear over finite sums, so the identity extends to all of H.

2.1step 1.2F4F5algebra

[(c) for Lz] With F=Ylm, ∂φF=imF, so DF=−iℏ(im)F=ℏmF=LzF by the definition of Lz. For negative indices [F4] gives Yl,−m=(−1)mYlm‾ and ∂φYl,−m=(−1)m∂φYlm‾=−im Yl,−m, so again DYl,−m=ℏ(−m)Yl,−m=LzYl,−m. Extending by linearity to H completes the proof of (c).

3.1step 1.1step 2.1F1F6algebra∎

[(d)] On the basis ALzYlm=l(l+1)ℏm Ylm=LzAYlm, and both operators are linear, so A and Lz commute on the finite span H. For each l, the span Vl of Yl,−l,…,Yll consists of eigenvectors of A with eigenvalue l(l+1) and of Lz with eigenvalues ℏ(−l),…,ℏl; by [F6] the projection is Plx=∑m⟨x,Ylm⟩Ylm, and grouping the diagonal expansions by l shows APlx=PlAx for x∈D(A) and LzPlx=PlLzx for x∈D(Lz). Finally, let F∈D(A) satisfy the displayed product-weight condition. Since ∣m∣≤l and l2≤l2(l+1)2, F∈D(Lz). The product-weight condition gives LzF∈D(A) and AF∈D(Lz), so both compositions are defined: ALzF=∑l,ml(l+1)ℏm clmYlm with ∑∣l(l+1)ℏm∣2∣clm∣2<∞ by hypothesis, and symmetrically for LzAF, so the two series converge to the same element. For the maximal domains, the product-weight condition implies F∈D(Lz), since λl≥2 whenever m≠0; it is then exactly the condition LzF∈D(A). Hence it describes D(ALz) alone. The other composition additionally requires F∈D(A), giving D(LzA)=D(A)∩D(ALz). The stated m=0 counterexample shows these maximal domains can differ. To verify that the product condition is not automatic on D(A)∩D(Lz), take cll=l−3 for l≥1 and every other coefficient zero: the separate domain sums are bounded by constants times ∑l−2 and ∑l−4, respectively, but the product sum has terms ℏ2(l+1)2/l2≥ℏ2 and diverges.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Harmonic polynomial dimension

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). For k∈N let Pk be the real vector space of homogeneous polynomial functions of degree k on R3, with Pk={0} for k<0. Put r2=x12+x22+x32, Δ=∂12+∂22+∂32, and Hk=ker⁡(Δ∣Pk). Define their explicit complexifications PkC:={p+iq:p,q∈Pk} and HkC:={p+iq:p,q∈Hk}. The Laplacian and real rotation pullbacks extend complex-linearly to these spaces. Then:

(a) for every l≥2, Δ:Pl→Pl−2 is onto and Pl=Hl⊕r2Pl−2 as real vector spaces; the same decomposition holds over C after complexification;

(b) for every l≥0, dim⁡RPl=(l+22) and dim⁡RHl=2l+1, and these equal dim⁡CPlC and dim⁡CHlC, respectively;

(c) for l≥0 and 0≤m≤l, using the exact positive constants blm and Legendre polynomials of Associated Legendre function normalization, set, initially for x≠0, hlm(x):=(−1)mblm(2π)−1/2(x1+ix2)mrl−mPl(m)(x3/r). This expression has a unique polynomial extension to R3, with value 0 at the origin for l>0 and the nonzero constant b00/2π=1/4π for l=0. It is a nonzero element of HlC, and hlm∣S2=Ylm. For 1≤m≤l define hl,−m=(−1)mhlm‾. The 2l+1 functions (hlm)−l≤m≤l form a complex basis of HlC, and their restrictions form a complex basis of its restriction space. The restriction maps on Hl and HlC are injective; a corresponding real basis of Hl is hl0 together with 2Re⁡hlm and 2Im⁡hlm for 1≤m≤l.

(d) let ω,m0,ℏ>0, ℓ=ℏ/(m0ω), and let (ψN)N∈N3 be the complex Hilbert-space eigenbasis of the realised isotropic oscillator in Multidimensional oscillator completeness, with HψN=ℏω(∣N∣+32)ψN, ∣N∣=N1+N2+N3. Put EM:=span⁡C{ψN:∣N∣=M}. Then:

(d1) each f∈EM has a unique representation f=Pe−r2/(2ℓ2) with a complex polynomial P of degree at most M, and the map ΦM:EM→PMC assigning its degree-M homogeneous part (zero if its degree is lower) is a complex-linear isomorphism. It restricts to a real-linear isomorphism from the real span of the same ψN onto PM;

(d2) every rotation R∈SO(3) preserves EM. With π(R)f=f(R−1⋅) on EM and ρ(R)p=p(R−1⋅) on PMC, one has ΦMπ(R)=ρ(R)ΦM; the same identity holds on the real subspaces;

(d3) one has the decomposition into rotation-invariant complex subspaces EM=⨁j=0⌊M/2⌋ΦM−1(r2jHM−2jC), and the analogous real decomposition. The harmonic orbital degrees in this level are exactly l=M,M−2,M−4,⋯≥0, each occurring once as a harmonic-degree summand. In particular not every l≤M occurs, and ∑j=0⌊M/2⌋(2(M−2j)+1)=(M+22) is the oscillator-level multiplicity.

All assertions concern polynomials and the realised oscillator, with no physical interpretation asserted.

Facts & Assumptions

Given: The real and complex polynomial spaces, r2, Δ, the exact Ylm convention, and the oscillator parameters and basis.

[F1]

Product differentiation gives Δ(fg)=fΔg+2∇f⋅∇g+gΔf, Δr2=6, and Euler's identity x⋅∇q=dq for a degree-d homogeneous polynomial (Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, homogeneous polynomial and homogeneous ideal). Sine and cosine satisfy sin⁡′=cos⁡ and cos⁡′=−sin⁡ (published result The derivatives of sine and cosine are cosine and minus sine). For totally differentiable maps f:U→V and g:V→Rp, the total derivative satisfies D(g∘f)(a)=Dg(f(a))∘Df(a) when the composite is defined (published result The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). The Euclidean inverse function theorem states that a C1 map with invertible total derivative at a point has a local C1 inverse g satisfying Dg(y)=Df(g(y))−1 (published result The Euclidean inverse function theorem). The polar-chart map T(r,θ,φ)=rΦ(θ,φ) from (0,∞)×(0,π)×(0,2π) onto R3∖K, where K:={(x,0,z):x≥0}, is a C1 diffeomorphism: this is proved in step 2.1 of Sphere integrals in angle coordinates; its determinant r2sin⁡θ is nonzero on that open domain. Thus its inverse coordinates r,θ,φ are C1 there.

[F3]

Pl(t)=(2ll!)−1(d/dt)l(t2−1)l has degree l and nonzero leading coefficient (Legendre Rodrigues equation). The supplied associated function is wlm=(−1)m(1−t2)m/2Pl(m), and Ylm=blmwlm(cos⁡θ)eimφ/2π, with Yl,−m=(−1)mYlm‾ and the full family orthonormal in complex L2(S2) (Associated Legendre function normalization, Sphere ONB). Its angular differential equation gives ΔS2Ylm=−l(l+1)Ylm (Sphere Laplacian diagonal domain, Sphere integrals in angle coordinates, The derivatives of sine and cosine are cosine and minus sine, The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a), The Euclidean inverse function theorem).

[F4]

ψN=QNe−r2/(2ℓ2) with QN a real product of scaled Hermite polynomials, degree ∣N∣, and degree-∣N∣ part a nonzero real multiple of xN. The oscillator is self-adjoint and diagonal in the complete product basis, with eigenvalue λ∣N∣=ℏω(∣N∣+3/2); the energies differ exactly when the total degrees differ, not whenever the multi-indices differ. On the finite product span its action is H=−ℏ2Δ/(2m0)+m0ω2r2/2 (Multidimensional oscillator completeness, Oscillator scaling and exact realization, Hermite recurrence, Hermite orthogonality).

Proof

1.1F1givenalgebra

[An explicit Laplacian right inverse with the correct degree] Fix l≥2 and q∈Pl−2. Let J=⌊(l−2)/2⌋ and Aj=2(j+1)(2l−2j−1) for 0≤j≤J. Every Aj>0. Define c0=(4l−2)−1 and cj=−cj−1/Aj for 1≤j≤J, and put Slq=∑j=0Jcjr2j+2Δjq. Every term has degree 2j+2+(l−2−2j)=l, and the whole polynomial lies in r2Pl−2. More generally [F1] gives Δ(r2tp)=2t(2d+2t+1)r2t−2p+r2tΔp for p∈Pd: compute ∇r2t=2tr2t−2x and Δr2t=2t(2t+1)r2t−2 and substitute Euler's identity. For t=j+1, p=Δjq this is Δ(r2j+2Δjq)=Ajr2jΔjq+r2j+2Δj+1q. In ΔSlq, the coefficient of q is c0A0=1, each coefficient of r2jΔjq for j≥1 is cjAj+cj−1=0, and the last leftover vanishes since ΔJ+1q=0. Hence ΔSlq=q.

1.2F3algebra

[Legendre parity makes the solid expression a polynomial] Rodrigues' formula in [F3] gives Pl(−t)=(−1)lPl(t) by differentiating an even polynomial l times. After m more derivatives, Pl(m)(−t)=(−1)l−mPl(m)(t), and its degree is l−m with nonzero leading coefficient for 0≤m≤l. Thus if Pl(m)(t)=∑kdktk, every nonzero term has 0≤k≤l−m and k≡l−m(mod2). Consequently rl−mPl(m)(x3/r)=∑kdkx3k(r2)(l−m−k)/2 is a polynomial: every exponent of r2 is a nonnegative integer. Multiplication by (−1)mblm(x1+ix2)m/2π makes the displayed hlm homogeneous of degree l. This polynomial determines the unique extension from x≠0, taking value zero for l>0 and the nonzero constant b00/2π for l=0; since P0=1 and ∫−11P02=2, this constant is 1/4π. In angle coordinates x1+ix2=rsin⁡θeiφ, so the factor (−1)m gives exactly hlm=rlYlm, including the required convention.

1.3F2F4algebra

[Polynomial Gaussian core and the top-degree isomorphism] Write G(x)=e−r2/(2ℓ2). The leading terms in [F4] show by induction on degree that every complex polynomial P of degree at most D is a finite linear combination of the QN with ∣N∣≤D: subtract the matching nonzero multiples of all leading monomials of degree D and apply the induction to the polynomial of lower degree. The same argument uses real coefficients when P is real. Hence every PG belongs to the oscillator's finite product core. An element of EM is uniquely PG with degree at most M, because G never vanishes; equality of the continuous functions as L2 classes implies equality everywhere and thus polynomial uniqueness. On the basis (ψN)∣N∣=M, the degree-M map ΦM takes each vector to a nonzero real multiple of its distinct monomial xN. These monomials are a basis of PMC, so ΦM is directly an isomorphism, and its restriction to real spans is the stated real isomorphism. This avoids the incorrect assertion that different multi-indices necessarily have different energies.

2.1step 1.1F2algebra

[Direct sum and all dimensions] Multiplication by r2 is injective on polynomials, so r2Pl−2 and Pl−2 have the same finite dimension. Step 1.1 proves that Δ:r2Pl−2→Pl−2 is onto, hence bijective by [F2]. For p∈Pl, p−Sl(Δp)∈Hl and Sl(Δp)∈r2Pl−2, proving the sum. If a vector in the intersection is r2q and harmonic, bijectivity makes q=0, proving directness. Thus dim⁡Hl=(l+22)−(l2)=2l+1 for l≥2; for l=0,1, every polynomial is harmonic and the dimensions are 1,3. Applying the real decomposition separately to real and imaginary parts proves the complex decomposition. A real basis remains a complex basis of its explicit complexification, since taking real and imaginary parts verifies both spanning and independence, so the complex dimensions agree.

2.2step 1.3F1F4algebra

[Rotations preserve oscillator levels before intertwining] Let f=PG∈EM and R∈SO(3). Since ∣R−1x∣=∣x∣, its pullback is (P∘R−1)G, a polynomial Gaussian of degree at most M; step 1.3 places it in the finite product core, hence in D(H). The chain rule and R−1(R−1)T=I give Δ(f∘R−1)=(Δf)∘R−1, and multiplication by r2 also commutes with the pullback. The exact differential expression of [F4] on the core therefore gives Hπ(R)f=π(R)Hf=λMπ(R)f. Expand the rotated vector in its finite product basis. Since the energies are λK=ℏω(K+3/2) for total degree K and are distinct when K≠M, the eigen-equation kills every coefficient outside total degree M. Thus π(R)f∈EM, and applying R−1 proves equality of the preserved space. Only after establishing this invariance define π(R) on EM. An invertible linear substitution preserves each homogeneous degree, so the degree-M part of P∘R−1 is (ΦMf)∘R−1. Hence ΦMπ(R)=ρ(R)ΦM, also on the real subspaces.

3.1step 2.1step 1.2F1F3algebra

[Harmonicity, complex basis and its real counterpart] On the open polar-coordinate set r>0, 0<θ<π, [F1] gives that the inverse coordinates are C1 and that their derivative is the inverse Jacobian of T. Differentiating the coordinate map with the sine and cosine derivatives in [F1], the Jacobian of T has columns ∂rT=(sin⁡θcos⁡φ,sin⁡θsin⁡φ,cos⁡θ), ∂θT=(rcos⁡θcos⁡φ,rcos⁡θsin⁡φ,−rsin⁡θ), and ∂φT=(−rsin⁡θsin⁡φ,rsin⁡θcos⁡φ,0); its determinant is r2sin⁡θ≠0. Matrix inversion gives the inverse Jacobian rows, hence ∂x1=sin⁡θcos⁡φ ∂r+cos⁡θcos⁡φr ∂θ−sin⁡φrsin⁡θ ∂φ, ∂x2=sin⁡θsin⁡φ ∂r+cos⁡θsin⁡φr ∂θ+cos⁡φrsin⁡θ ∂φ, and ∂x3=cos⁡θ ∂r−sin⁡θr ∂θ. Squaring these operators and adding with the product rule, the mixed terms cancel, the second-derivative coefficients are 1, r−2, and (r2sin⁡2θ)−1, and the first-derivative coefficients are 2r−1, r−2cot⁡θ, and 0. Thus Δ=∂r2+2r∂r+1r2(∂θ2+cot⁡θ ∂θ+1sin⁡2θ∂φ2). Applying this identity to hlm(x)=rlYlm(θ,φ) gives Δhlm=rl−2(l(l+1)Ylm−SYlm)=0 by [F3], where S is the angular differential expression from Sphere Laplacian diagonal domain, Sphere integrals in angle coordinates, The derivatives of sine and cosine are cosine and minus sine, The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a), The Euclidean inverse function theorem and SYlm=l(l+1)Ylm. The polynomial Δhlm is continuous and vanishes on this dense set, hence vanishes everywhere. Its real and imaginary parts are therefore real harmonic polynomials. Negative indices follow by complex conjugation, and nonvanishing follows from the norm-one restriction Ylm. Restriction is injective: if a homogeneous polynomial vanishes on S2, then p(x)=∣x∣lp(x/∣x∣)=0 for x≠0, and continuity also gives its value at zero. Thus the 2l+1 orthonormal restrictions are linearly independent in a complex space of dimension 2l+1 by step 2.1 and hence form a basis. The invertible pairwise transformation from hlm,hl,−m to 2Re⁡hlm,2Im⁡hlm, together with real hl0, produces 2l+1 real harmonic polynomials which are complex-linearly independent, hence real-linearly independent and a real basis.

3.2step 2.1step 1.3step 2.2F1F2F4algebra

[Full oscillator-level decomposition and multiplicities] Iterating the proved direct decomposition of step 2.1 yields PM=⨁j=0⌊M/2⌋r2jHM−2j: at each substitution directness is preserved because multiplication by r2 is injective. Complexifying gives the same result for PMC. A rotation commutes with the Laplacian, preserves degree, and preserves r2, so each displayed harmonic-degree summand is invariant. Transfer it through the intertwining isomorphism of steps 1.3 and 2.2 to obtain exactly the stated decomposition of EM, with one copy of each degree M,M−2,…. All listed degrees occur because their dimensions 2(M−2j)+1 are positive; no other degree is part of the exhaustive direct decomposition. Their sum is (M+22) by the polynomial-space dimension, and equals the number of oscillator basis vectors with ∣N∣=M by [F2] and [F4]. The same transfer over R gives the claimed real decomposition.

4.1step 1.1step 2.1step 1.2step 3.1step 1.3step 2.2step 3.2given∎

[Conclusion] Steps 1.1 and 2.1 prove (a) and (b), steps 1.2 and 3.1 prove (c), and steps 1.3, 2.2 and 3.2 prove (d1)–(d3) with the full real and complex spaces explicit. Countable Choice is retained for the invoked sphere and realised-oscillator suppliers; the finite polynomial recursions, finite basis changes and fixed top-degree maps require no additional choice principle.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Strongly continuous rotation action

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let S2:={x∈R3:∣x∣=1}, let σ be the polar surface measure on S2 (The polar surface set function on the unit sphere), a finite Borel measure with σ(S2)=4π (Sphere integrals in angle coordinates), and let L2(S2,σ) be the complex Hilbert space of equivalence classes of square-integrable functions with inner product linear in the first argument (The space Lp(μ) as the quotient by null functions, Complex Lp classes and Euclidean test-function conventions, Hilbert space, Real and complex inner product spaces, with the inner product linear in the first argument). Let (Ylm)l≥0, −l≤m≤l be the complete orthonormal family of Associated Legendre function normalization and Sphere ONB, put Vl:=span⁡{Yl,−l,…,Yll} and let H be the algebraic span of all Ylm. For l≥0 let Hlpol:={ p:R3→C: p is a homogeneous polynomial of degree l with Δp=0 } be the complex vector space of harmonic homogeneous polynomials of degree l (homogeneous polynomial and homogeneous ideal, Monomials, coefficients, degree in each variable and total degree in F[x1,…,xn], The Laplacian of a C2 function and of a C2 vector field), let resl:Hlpol→L2(S2,σ) send p to the class of the restriction p∣S2, and let R∈SO(3) range over the special orthogonal group, a subgroup of the general linear group with R−1=RT and unit operator norm, the set M3(R) being identified with R9 and given the Euclidean metric (The orthogonal group O(3) and the rotation group SO(3), SO(3) is a subgroup of the general linear group, with R−1=RT). For F∈L2(S2,σ) define U(R)F to be the class of the function Ω↦F~(R−1Ω), where F~ is any representative of F. Then:

(a) U(R) is a well-defined C-linear isometry of L2(S2,σ) onto itself with inverse U(R−1), hence a unitary operator (Linear isometries, and orthogonal or unitary operators on finite-dimensional inner product spaces); U(I3)=id and U(RS)=U(R)U(S) for all R,S∈SO(3), so R↦U(R) is a group homomorphism from SO(3) into the group of unitary operators on L2(S2,σ);

(b) for every R∈SO(3) and p∈Hlpol one has p∘R−1∈Hlpol; the map ρ(R):p↦p∘R−1 is a C-linear automorphism of Hlpol with inverse ρ(R−1), and resl(ρ(R)p)=U(R) resl(p); consequently U(R)Vl=Vl for every l≥0, and U(R)H=H;

(c) (strong continuity) for every F∈L2(S2,σ) and every R0∈SO(3), for every ε>0 there is δ>0 such that ∥R−R0∥F<δ for the Euclidean (Frobenius) norm on M3(R) implies ∥U(R)F−U(R0)F∥2<ε; thus every orbit map R↦U(R)F is continuous (Continuity of a map between metric spaces, at a point and globally, in the ε-δ form). More precisely, if PNF:=∑l≤N∑m⟨F,Ylm⟩Ylm denotes the orthogonal truncation of F to ⨁l≤NVl (Orthonormal families, complete orthonormal systems and Hilbert bases), then for all R ∥U(R)F−F∥2≤2 ∥F−PNF∥2+∥U(R)PNF−PNF∥2.

All assertions are mathematical statements about the action of SO(3) on L2(S2,σ); no physical interpretation of the rotation action is asserted.

Facts & Assumptions

Given: the sphere S2, the polar measure σ, the complete orthonormal family (Ylm), the blocks Vl, the finite span H, the spaces Hlpol, the restriction maps resl, and the group SO(3).

[F1]

The polar surface measure equals the chart surface measure on S2, is finite with σ(S2)=4π, and is preserved by orthogonal transformations: σ(OE)=σ(E) for every orthogonal O and every Borel E⊆S2 (Agreement with the existing polar sphere measure, The polar surface set function on the unit sphere, Sphere integrals in angle coordinates).

[F2]

(Ylm) is a complete orthonormal family of L2(S2,σ); each Vl has the orthonormal basis (Yl,−l,…,Yll); Vl⊥Vk for l≠k; H is dense; and the orthogonal truncations PNF converge to F in norm (Parseval) (Sphere ONB, Orthonormal families, complete orthonormal systems and Hilbert bases, Parseval equivalences for an orthonormal family).

[F3]

Each Ylm has a continuous polynomial representative on S2 by [F4] and Associated Legendre function normalization. For 0≤m≤l, its chart formula is Ylm(θ,φ)=blmwlm(cos⁡θ)eimφ/2π; for 1≤m≤l, Yl,−m=(−1)mYlm‾. Thus the negative-index formula includes the factor (−1)m, not the phase-free positive-index expression. Continuity at chart boundaries follows from the polynomial representative, and conjugation preserves it. Every F∈H is consequently continuous on the whole sphere.

[F4]

The restriction map resl is a C-linear isomorphism Hlpol→Vl, so every F∈Vl has a unique representative pF∈Hlpol with resl(pF)=F (Harmonic polynomial dimension (b),(c)).

[F5]

Polynomial maps on R3 are C∞; for a C2 function p and the linear map L(x)=R−1x the chain rule gives J(p∘L)(x)=Jp(Lx) R−1 and therefore Δ(p∘L)(x)=tr⁡((R−1)TD2p(Lx)R−1)=Δp(Lx), because R−1(R−1)T=I3; the top-degree part of p∘L is the composition of the top-degree part of p with L, so invertibility of L preserves the degree of a nonzero polynomial and preserves homogeneity; sums and scalar multiples of harmonic homogeneous polynomials of one degree are again harmonic homogeneous of that degree (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a), A total derivative computes every directional derivative, and its matrix is the Jacobian, The Laplacian of a C2 function and of a C2 vector field, Directional derivatives and partial derivatives of a map U⊆Rm→Rn, Ck Euclidean maps and diffeomorphisms, Ck Euclidean maps are closed under componentwise algebra and composition, homogeneous polynomial and homogeneous ideal, Monomials, coefficients, degree in each variable and total degree in F[x1,…,xn]).

[F6]

Let Φ:S2→S2 be a measurable bijection with σ(Φ−1E)=σ(E) for all Borel E; then for F,G∈L2(S2,σ) the functions F∘Φ, G∘Φ lie in L2(S2,σ), with ⟨F∘Φ,G∘Φ⟩=⟨F,G⟩ (substitution formula for the measure-preserving bijection), and F=G σ-a.e. implies F∘Φ=G∘Φ σ-a.e.; these identities are those of the change-of-variables formula applied to indicators and to nonnegative Borel functions and then to the four positive and negative parts (Agreement with the existing polar sphere measure, The space Lp(μ) as the quotient by null functions, Complex Lp classes and Euclidean test-function conventions).

[F7]

A continuous function on the compact metric space S2 is uniformly continuous (Every continuous map from a nonempty compact Hausdorff space to a uniform space is uniformly continuous); if g:S2→C is continuous and bounded then ∫S2∣g∣2 dσ≤4π ∥g∥∞2, because the nonnegative integral is monotone and σ(S2)=4π (The nonnegative Lebesgue integral, [F1]).

[F8]

Euclidean norms on M3(R)≅R9 and on R3: for the Frobenius norm ∥A∥F=(∑i,jAij2)1/2, ∣Ax∣≤∥A∥F∣x∣ for all x∈R3 by the Cauchy-Schwarz inequality; for R∈SO(3), ∥R∥F=3 and ∥R∥=1, and R−1−R0−1=R−1(R0−R)R0−1 for R,R0∈SO(3) (Cauchy–Schwarz: ∣⟨u,v⟩∣≤∥u∥∥v∥, with equality exactly for linearly dependent vectors, SO(3) is a subgroup of the general linear group, with R−1=RT, The orthogonal group O(3) and the rotation group SO(3)).

[F9]

SO(3) is a subgroup of the general linear group with R−1=RT, so (RS)−1=S−1R−1 and I3−1=I3 for all R,S∈SO(3) (SO(3) is a subgroup of the general linear group, with R−1=RT, The orthogonal group O(3) and the rotation group SO(3)).

Proof

1.1F1F6algebra

[(a), substitution map] Let R∈SO(3) and Φ:=R−1. Then Φ maps S2 bijectively onto S2, since ∣R−1y∣=∣y∣, and by [F1] with the orthogonal matrix R one has σ(Φ−1E)=σ(RE)=σ(E) for every Borel E. If F~=F~′ is any representative pair of a class F∈L2, then {F≠F′} has σ-measure zero, hence so does R{F≠F′} by [F1], so all representatives produce the same class: U(R) is well defined; and U(R)F∈L2 with ∥U(R)F∥2=∥F∥2, ⟨U(R)F,U(R)G⟩=⟨F,G⟩, by [F6].

1.2F5algebra

[(b), harmonicity] Let p∈Hlpol and R∈SO(3); by [F5] with L(x)=R−1x, p∘R−1 is C∞ with Δ(p∘R−1)(x)=Δp(R−1x)=0 for all x, and since R−1 is an invertible linear substitution, p∘R−1 is a homogeneous polynomial of degree l (for p=0 it is 0); hence p∘R−1∈Hlpol. The assignment ρ(R) is C-linear because composition with a fixed map is linear, and ρ(R)ρ(R−1)=ρ(R−1)ρ(R)=id because p∘R−1∘R=p and p∘R∘R−1=p; so ρ(R) is a linear automorphism of Hlpol.

1.3F3F7F8algebra

[(c), continuity on the finite span] Let F∈H and R0∈SO(3). By [F3] F is continuous on S2, hence uniformly continuous by [F7]: given ε there is η>0 with ∣x−y∣<η, x,y∈S2, implying ∣F(x)−F(y)∣<ε/(24π). For Ω∈S2 and R∈SO(3), U(R)F(Ω)−U(R0)F(Ω)=F(R−1Ω)−F(R0−1Ω) and, by [F8], ∣R−1Ω−R0−1Ω∣=∣R−1(R0−R)R0−1Ω∣≤∥R−1∥op∥R0−R∥F∥R0−1Ω∥=∥R−R0∥F. Choosing δ:=η one gets sup⁡S2∣U(R)F−U(R0)F∣≤ε/(24π) whenever ∥R−R0∥F<δ, and then [F7] gives ∥U(R)F−U(R0)F∥2<ε/2<ε.

2.1step 1.1F9algebra

[(a), group law and unitarity] For F∈L2 and Ω∈S2, (U(R−1)U(R)F)(Ω)=(U(R)F)(RΩ)=F~(R−1RΩ)=F~(Ω) and symmetrically U(R)U(R−1)=id, so U(R) is a bijective linear isometry, i.e. unitary, with inverse U(R−1) (Linear isometries, and orthogonal or unitary operators on finite-dimensional inner product spaces); U(I3)=id since R−1=RT=Ω-fixing for R=I3; (U(R)U(S)F)(Ω)=(U(S)F)(R−1Ω)=F~(S−1R−1Ω)=F~((RS)−1Ω)=(U(RS)F)(Ω) by [F9] the inverse identity (RS)−1=S−1R−1.

3.1step 2.1step 1.2F4algebra

[(b), intertwining and blocks] Let F∈Vl with representative p=pF∈Hlpol ([F4]) and let R∈SO(3). For Ω∈S2, U(R)F(Ω)=F~(R−1Ω)=p(R−1Ω)=(p∘R−1)(Ω): here R−1Ω lies in S2, so the identity of classes resl(ρ(R)p)=U(R)resl(p) follows. Hence U(R)Vl⊆Vl for every l by step 1.2 and [F4], and applying U(R−1) (step 2.1) gives equality; since H is the span of the Vl and U(R) is linear, U(R)H=H.

4.1step 2.1step 1.3F2algebra∎

[(c), the estimate and general F] For arbitrary F∈L2(S2,σ) and any R,R0, the unitarity of step 2.1 gives ∥U(R)F−U(R0)F∥2=∥U(R)(F−PNF)+U(R)PNF−U(R0)PNF+U(R0)(PNF−F)∥2≤2∥F−PNF∥2+∥U(R)PNF−U(R0)PNF∥2, with PNF∈H by [F2]. Given ε>0, choose N with ∥F−PNF∥2<ε/4 by [F2], then by step 1.3 a δ>0 with ∥U(R)PNF−U(R0)PNF∥2<ε/2 when ∥R−R0∥F<δ; then ∥U(R)F−U(R0)F∥2<ε. Taking R0=I3 and keeping the terms separately yields the displayed estimate ∥U(R)F−F∥2≤2∥F−PNF∥2+∥U(R)PNF−PNF∥2.

  • The orbit continuity proved in (c) is the strong-continuity condition of the definition of a unitary representation (Strongly continuous unitary representations, invariant linear subspaces and intertwiners) for the group SO(3) equipped with the Euclidean-metric topology inherited from M3(R)≅R9; this identification of vocabulary is not used in the proof of the present batch, whose mathematical content is exactly clauses (a)-(c).
  • Clause (b) is the statement that rotations permute the angular-momentum blocks Vl and act on each of them by the finite-dimensional automorphism ρ(R) transported by the restriction isomorphism resl; it is the algebraic input for the self-adjoint generators of the following item.
TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Orbital self-adjoint generators

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Use the notation of Strongly continuous rotation action: S2, the polar measure σ with σ(S2)=4π, the Hilbert space L2(S2,σ), the complete orthonormal family (Ylm), the blocks Vl, the finite span H, the spaces Hlpol of complex harmonic homogeneous polynomials of degree l, the restriction maps resl, and the unitary action U(R)F=class of F~(R−1⋅) of SO(3). Fix ℏ>0 and for i=1,2,3 define the rotation matrices R1(t):=(1000cos⁡t−sin⁡t0sin⁡tcos⁡t),R2(t):=(cos⁡t0sin⁡t010−sin⁡t0cos⁡t),R3(t):=(cos⁡t−sin⁡t0sin⁡tcos⁡t0001), put Ui(t):=U(Ri(t)). Let Lz denote the diagonal operator of Sphere Laplacian diagonal domain, set S:=R2(−π/2), W:=R1(π/2), and define L3:=Lz,L1:=U(S)−1LzU(S),L2:=U(W)−1LzU(W), with domains D(L1)=U(S)−1D(Lz), D(L2)=U(W)−1D(Lz) and D(L3)=D(Lz). Here e−itLz/ℏ is the explicit diagonal unitary sending Ylm to e−itmYlm; e−itL1/ℏ and e−itL2/ℏ mean its conjugates by U(S) and U(W), respectively. The proof identifies these operators with the complete strong derivative generators, so this construction has the same generator meaning without an appeal to an unrestricted spectral or Stone theorem. For a polynomial p put Dip:=(x×∇p)i=∑j,k=13εijk xj ∂kp, where εijk is the sign of the permutation (1,2,3)↦(i,j,k) and is 0 when two indices coincide. Then:

(a) each Li is a self-adjoint operator on L2(S2,σ) with U(Ri(t))=e−itLi/ℏ for all real t; the finite span H is contained in D(Li) and is a core of Li (Densely defined, closed and closable operators, and cores, Symmetric, self-adjoint and essentially self-adjoint operators); and for F∈Vl with polynomial representative pF∈Hlpol one has LiF=resl(−iℏ DipF), that is, (LiF)(Ω)=−iℏ DipF(Ω) for σ-almost every Ω∈S2; in particular LiYlm∈Vl for all l,m and LiH⊆H;

(b) on H the operators satisfy the commutation relations LiLj−LjLi=iℏ∑k=13εijkLk;

(c) L3=Lz as self-adjoint operators (equal domains and equal action); in particular L3Ylm=ℏmYlm, so that L3 is diagonal in the basis (Ylm), and for every F∈H and chart point (θ,φ) one has (L3F)(θ,φ)=−iℏ ∂φF(θ,φ);

(d) rotational conjugation identifies the domains: with the two quarter-turn rotations S:=R2(−π/2) and W:=R1(π/2) one has U(S)L1U(S)−1=L3,U(W)L2U(W)−1=L3, as self-adjoint operators, and consequently U(S)D(L1)=D(L3) and U(W)D(L2)=D(L3): the three domains are unitarily identified.

All assertions are mathematical statements about operators on L2(S2,σ); no physical interpretation of the generators as angular momentum is asserted, and no measurement postulate is used.

Facts & Assumptions

Given: the sphere, its measure, the harmonic basis, the blocks Vl, their finite span H, the rotation action, ℏ>0, and the matrices and explicitly defined operators in the Statement. Inner products are linear in the first argument.

[F1]

Strongly continuous rotation action supplies the unitary group law, strong continuity, U(R)Vl=Vl, U(R)H=H and polynomial intertwining. Harmonic polynomial dimension supplies the complex restriction isomorphism resl:Hlpol→Vl.

[F2]

Sphere Laplacian diagonal domain supplies LzYlm=ℏmYlm with full domain {F:∑l,mm2∣⟨F,Ylm⟩∣2<∞}. Sphere ONB and Parseval give the complete countable basis, mutually orthogonal blocks and convergence of P≤N:=∑l≤NPl to the identity. Associated Legendre function normalization gives the azimuthal factor eimφ, also for negative m by Yl,−m=(−1)mYlm‾.

[F3]

Direct multiplication of the displayed matrices, the sine/cosine addition and derivative formulas, and the quarter-turn values give Ri′(s)=JiRi(s)=Ri(s)Ji, where Jix=ei×x; Ri(−s)′=−JiRi(−s); SR1(t)S−1=R3(t) and WR2(t)W−1=R3(t). These are finite matrix identities, not generator identities.

Proof

1.1F2algebra

[(a), the complete diagonal domain and adjoint] Write clm=⟨F,Ylm⟩. On the domain of [F2], Parseval defines AF=∑l,mℏmclmYlm with convergence in norm, and finite sums prove symmetry. If G∈D(A∗) and K=A∗G, testing F=Ylm in ⟨AF,G⟩=⟨F,K⟩ gives ⟨K,Ylm⟩=ℏm⟨G,Ylm⟩. Parseval then implies ∑ℏ2m2∣⟨G,Ylm⟩∣2=∥K∥2<∞: G∈D(A) and K=AG. Conversely, for such G, the equality holds first on finite sums and then for all F∈D(A) by their coefficient truncations. Thus A=A∗ on exactly the full domain of [F2], and A=Lz. For F∈D(A), both P≤NF→F and AP≤NF→AF by the two convergent coefficient square sums; hence H is a graph core. All sums use the fixed countable index family; no fiber selection or uncountable basis choice occurs.

2.1step 1.1F1algebra

[(a),(d), adjoints and cores of the conjugates] If Q is either U(S) or U(W), define B=Q−1AQ on Q−1D(A). For G,K, the adjoint identity ⟨BF,G⟩=⟨F,K⟩ for every F∈Q−1D(A) is, on putting X=QF, exactly ⟨AX,QG⟩=⟨X,QK⟩ for every X∈D(A). By step 1.1 this holds exactly when QG∈D(A) and QK=AQG. Thus B∗=B with this full domain. Since QH=H, unitary transport of the graph approximations in step 1.1 gives H as a core of B. With L3=A, this proves self-adjointness, cores and the two stated conjugation/domain identities.

2.2step 1.1F2algebra

[(a), explicit groups and the full derivative domains] Define E(t)F=∑e−itmclmYlm. Parseval proves unitarity and the group law. For each F, a finite coefficient truncation controls the tail by 2∥F−P≤NF∥ and its finite head is continuous, proving strong continuity. If F∈D(A), ∣(e−itm−1)/t∣≤∣m∣ and convergence on finite heads with square-summable weighted tails give (E(t)F−F)/t→−iAF/ℏ. Conversely, if the quotient converges in norm to G, each coefficient tends to −imclm; Parseval for G gives ∑m2∣clm∣2=∥G∥2<∞ and G=−iAF/ℏ. Thus this is exactly the full strong derivative domain, not merely a core calculation. The same argument transported by Q gives exactly D(Q−1AQ) for Q−1E(t)Q, with derivative −iQ−1AQF/ℏ.

3.1step 2.2F1F2F3algebra

[(a),(c), identification with the actual rotation groups] Pullback by R3(t)−1 replaces φ by φ−t. Therefore U3(t)Ylm=e−itmYlm for every integer m by [F2], including negative m with its conjugation phase. Equality on the dense finite span, and boundedness, imply U3(t)=E(t) on all L2. The two finite matrix identities in [F3], transported by the group law, imply U1(t)=U(S)−1E(t)U(S) and U2(t)=U(W)−1E(t)U(W). These are precisely the explicitly defined exponentials in the Statement. Step 2.2 gives the exact derivative domains for all three groups and their actions −iLi/ℏ; it also proves uniqueness if an operator is defined by the strong derivative. In particular L3=Lz as full operators, without an AC Stone premise.

4.1step 2.2step 3.1F1F3F4algebra

[(a), polynomial differentiation with the correct inverse-curve sign] For F=resl(p)∈Vl, set gi(s,Ω)=∇p(Ri(−s)Ω)⋅JiRi(−s)Ω. The correct chain derivative is ddsp(Ri(−s)Ω)=−gi(s,Ω), directly from Ri(−s)′=−JiRi(−s) in [F3]. The fundamental theorem gives (Ui(t)F−F)/t=−t−1∫0tgi(s,⋅) ds. Uniform continuity on [−δ,δ]×S2 proves convergence uniformly, and hence in L2, to −gi(0,⋅)=−Dip∣S2. By step 2.2, F∈D(Li) and −iLiF/ℏ=−res(Dip), so LiF=res(−iℏDip). The curve p∘Ri(−t) lies in the finite-dimensional polynomial space Hlpol by [F1]; its coefficientwise derivative at zero also lies there, because homogeneity and the harmonic equations are finite linear coefficient constraints. This derivative is −Dip, proving that Dip is harmonic homogeneous of degree l, including the zero polynomial. Thus LiVl⊆Vl and LiH⊆H.

5.1step 1.1step 2.1step 3.1step 4.1F1F2algebra

[(a), finite blocks and their ODE identification] Put Ai(l)=Li∣Vl. This is an everywhere-defined Hermitian matrix: symmetry follows from self-adjointness on the whole space and Vl⊆D(Li). The quarter turns preserve Vl, so A1(l) and A2(l) are unitarily conjugate to the diagonal A3(l) with entries ℏm. Consequently ∥Ai(l)∥≤ℏl. Polynomial differentiation at time t, or the finite group law and step 4.1, gives (Ui(t)∣Vl)′=(−i/ℏ)Ai(l)Ui(t)∣Vl with initial value the identity. The absolutely convergent matrix exponential power series solves this equation; uniqueness follows by differentiating its inverse times any other solution. Hence the finite polynomial action equals exp⁡(−itAi(l)/ℏ). Orthogonal block truncation then identifies its all-vector extension with the groups of step 3.1. Explicitly D(Li)={F:∑l∥Ai(l)PlF∥2<∞} and LiF=∑lAi(l)PlF: this is the diagonal domain of step 1.1 transported by a block-preserving quarter turn. Thus no hidden narrower domain or additional time-continuity hypothesis is introduced.

5.2step 4.1F4algebra

[(b), commutators] The polynomial product rule gives [xa∂b,xc∂d]=xaδbc∂d−xcδad∂b. Substituting D1=x2∂3−x3∂2, D2=x3∂1−x1∂3 gives [D1,D2]=−D3, and cyclic permutation and antisymmetry give [Di,Dj]=−∑kεijkDk. All products are defined on H by step 4.1. Therefore [Li,Lj]=(−iℏ)2[Di,Dj]=ℏ2∑kεijkDk=iℏ∑kεijkLk.

6.1step 2.1step 3.1step 4.1F4algebra∎

[(c), the chart action] The chain rule gives ∂φp(sin⁡θcos⁡φ,sin⁡θsin⁡φ,cos⁡θ)=∇p(Ω)⋅(e3×Ω)=D3p(Ω). Apply step 4.1 and linearity to finite harmonic sums. Together with the full equality L3=Lz of step 3.1 and the domain conjugations of step 2.1, this proves every clause (a)-(d).

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Dipole angular selection rules

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Use the notation of Harmonic polynomial dimension: S2={x∈R3:∣x∣=1} with polar measure σ (Sphere integrals in angle coordinates, Complex Lp classes and Euclidean test-function conventions), the Hilbert space L2(S2,σ) with inner product linear in the first argument, the complete orthonormal family (Ylm) (Associated Legendre function normalization, Sphere ONB), the blocks Vl=span⁡{Yl,−l,…,Yll} and the finite span H, the spaces Hlpol of complex harmonic homogeneous polynomials of degree l, and the polynomial representatives hlm∈Hlpol of the Ylm, characterized by resl(hlm)=Ylm where resl(p) is the class of p∣S2 (homogeneous polynomial and homogeneous ideal, Monomials, coefficients, degree in each variable and total degree in F[x1,…,xn], The Laplacian of a C2 function and of a C2 vector field). Put H−1pol:={0}, fix l≥0 and m∈{−l,…,l}, write H:=hlm, and define the three polynomial multipliers Q0:=x3,Q±:=x1±ix22, so that Qq=∑i=13cq,ixi with c0,i=δi3 and c±,i=(δi1±i δi2)/2. Let Di and Lz,L3 be as in Orbital self-adjoint generators, where L3=Lz and LzYlm=ℏmYlm (Sphere Laplacian diagonal domain). Then, for each q∈{−1,0,1}:

(a) (harmonic decomposition) the polynomial Gq:=∑icq,i ∂iH2l+1 lies in Hl−1pol and Qq H=Rq+r2Gq,Rq:=QqH−r2Gq∈Hl+1pol, where r2(x)=x12+x22+x32; for l=0 one has Gq=0 and Rq=QqH;

(b) (degree selection) the function Qq∣S2⋅Ylm lies in Vl+1⊕Vl−1 (with V−1:={0}); consequently, for all l′≥0 and m′∈{−l′,…,l′}, ⟨Qq∣S2 Ylm, Yl′m′⟩=0whenever l′∉{l−1, l+1};

(c) (weight selection) as elements of H⊂L2(S2,σ) one has Lz(Qq∣S2 Ylm)=ℏ(m+q) (Qq∣S2 Ylm), because ∂φQq=iq Qq on S2 (Sine and cosine defined by their real power series, The complex exponential by its power series); hence in addition ⟨Qq∣S2 Ylm, Yl′m′⟩=0whenever m′≠m+q;

(d) (selection rules) combining (b) and (c): the angular matrix element ⟨Yl′m′, Qq∣S2 Ylm⟩ can be nonzero only if l′=l+1 and m′=m+q,orl≥1, l′=l−1 and m′=m+q; in particular only Δl=±1 and Δm∈{0,±1} occur with the component-dependent shift Δm=q for the multiplier Qq. These are necessary conditions on the angular integrals only: no assertion is made that an allowed pair has a nonzero integral, and no absorption, emission or transition probability is asserted or computed.

Facts & Assumptions

Given: the family (Ylm) with representatives hlm, the multipliers Qq, the fixed indices l,m, and the polynomial H=hlm.

[F1]

The restriction map resl:Hlpol→Vl is a C-linear isomorphism; (Ylm)m is an orthonormal basis of Vl; Vl⊥Vk for l≠k; (Ylm) is complete in L2(S2,σ) and H is dense (Harmonic polynomial dimension (b),(c), Sphere ONB, Orthonormal families, complete orthonormal systems and Hilbert bases, Parseval equivalences for an orthonormal family, Real and complex inner product spaces, with the inner product linear in the first argument).

[F2]

Δ is linear, satisfies the product rule Δ(fg)=fΔg+2∇f⋅∇g+gΔf, commutes with each ∂i, and annihilates the coordinate functions xi; r2 has Δr2=6 and ∇r2=2x; Euler's identity x⋅∇q=kq holds for homogeneous q of degree k; hence for homogeneous p of degree l and q of degree l−1: Δ(xip)=2∂ip and Δ(r2q)=(4l+2)q+r2Δq (The Laplacian of a C2 function and of a C2 vector field, Directional derivatives and partial derivatives of a map U⊆Rm→Rn, Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0, Ck Euclidean maps are closed under componentwise algebra and composition, A total derivative computes every directional derivative, and its matrix is the Jacobian, homogeneous polynomial and homogeneous ideal).

[F3]

Ylm(θ,φ)=blmwl∣m∣(cos⁡θ)eimφ/2π for m≥0 and Yl,−m=(−1)mYlm‾; in particular ∂φYlm=im Ylm in the angular chart, and Ylm is continuous (Associated Legendre function normalization (a),(c), The complex exponential by its power series, Sine and cosine defined by their real power series).

[F5]

Lz is self-adjoint on L2(S2,σ) with LzYlm=ℏmYlm, and for F∈H the chart values of LzF are −iℏ∂φF; eigenvectors of Lz for distinct eigenvalues are orthogonal, because ⟨LzF,G⟩=⟨F,LzG⟩ for F,G in the domain of the self-adjoint Lz (Sphere Laplacian diagonal domain (a),(c), Self-adjoint and normal endomorphisms of a finite-dimensional real or complex inner product space, Orthogonality and the orthogonal complement).

[F6]

L3=Lz (equal self-adjoint operators) by Orbital self-adjoint generators (c).

Proof

1.1F2algebra

[(a), the two harmonic pieces] By [F2] and ΔH=0, Δ(QqH)=∑icq,iΔ(xiH)=∑icq,i 2∂iH. Since ∂i commutes with Δ and ΔH=0, each ∂iH lies in Hl−1pol, so Gq=∑icq,i∂iH/(2l+1)∈Hl−1pol; by [F2] and Euler's identity for the homogeneous polynomial Gq of degree l−1, Δ(r2Gq)=6Gq+4(l−1)Gq+r2ΔGq=(4l+2)Gq+r2ΔGq, and here ΔGq=0, so Δ(r2Gq)=(4l+2)Gq=2∑icq,i∂iH. Hence ΔRq=Δ(QqH)−Δ(r2Gq)=0: Rq∈Hl+1pol, and Rq is homogeneous of degree l+1 as the difference of two such polynomials; for l=0 the derivative terms vanish because H is constant, so Gq=0 and Rq=QqH.

2.1step 1.1F1algebra

[(b), degree selection] Restricting the identity of step 1.1 to S2, where r2=1, gives Qq∣S2⋅Ylm=resl+1(Rq)+resl−1(Gq), where for l=0 the second term is omitted because Gq=0; for l≥1 both restriction maps are the defined maps in [F1]. This is an element of Vl+1+Vl−1, using H∣S2=Ylm ([F1]) and QqH∣S2=(Qq∣S2)Ylm; the sum is direct and each summand lies in the corresponding block, so orthogonality Vl+1⊥Vl′ and Vl−1⊥Vl′ for l′∉{l−1,l+1} ([F1]) gives ⟨Qq∣S2Ylm,Yl′m′⟩=0 for such l′.

3.1step 2.1F3F4F5algebra

[(c), weight selection] Let F:=Qq∣S2⋅Ylm∈H by step 2.1. On the chart, by [F3] and [F4], ∂φF=(∂φQq)Ylm+Qq ∂φYlm=iq F+im F=i(m+q)F; hence by [F5] LzF=−iℏ∂φF=ℏ(m+q)F in H. If now m′≠m+q, then using the self-adjointness of Lz ([F5]) and LzYl′m′=ℏm′Yl′m′, ℏ(m+q)⟨F,Yl′m′⟩=⟨LzF,Yl′m′⟩=⟨F,LzYl′m′⟩=ℏm′⟨F,Yl′m′⟩, so ⟨F,Yl′m′⟩=0.

4.1step 2.1step 3.1F6∎

[(d), the rules] Clause (d) restates steps 2.1 and 3.1 for the multiplier Qq: a nonzero angular integral requires l′∈{l−1,l+1} and m′=m+q; the case l′=l−1 presupposes l≥1 because −l′≤m′≤l′ with l−1≥0; and by [F6] the weight m is the L3-eigenvalue in units ℏ, so the multiplier Qq shifts the L3-eigenvalue by ℏq.

5 · Examples, experiments and counterexamples

None yet.

Sources