Masteries · Jason GlowneyArithmoí Book 1

Masteries

1.1.22 Hilbert Spaces

A Hilbert space, named after the math­e­mati­cian David Hilbert, ex­tends the meth­ods of vector al­ge­bra and cal­cu­lus from two and three-di­men­sion­al Eu­clid­ean spaces to spaces ∈ℝ,ℂ\in \mathbb{R},\mathbb{C}, of finite and in­fi­nite di­men­sions. A Hilbert space is an ab­stract vector space that is com­plete and has the struc­ture of an inner prod­uct to make pos­si­ble the mea­sure­ment of lengths and angles. In this sec­tion, we use many of the de­f­i­n­i­tions from Ap­plied Analy­sis by Hunter and Nachter­gaele 1, which is well-writ­ten with a straight­for­ward ap­proach to many of the analy­sis topics con­sid­ered in this book. To start, we must define what the inner prod­uct of a com­plex linear space XX is.

De­f­i­n­i­tion (1.1.137) – An inner prod­uct on a com­plex linear space XX is a map

⟨∙,∙⟩:X×X→ℂ\begin{array}{r} \left\langle \bullet , \bullet \right\rangle:X \times X\rightarrow \mathbb{C} \end{array}

such that, for all xj∈Xx_{j} \in X and a,b∈ℂa,b\in \mathbb{C}:

  • ⟨x1,ax2+bx3⟩=a⟨x1,x2⟩+b⟨x1,x3⟩\left\langle x_{1},ax_{2} + bx_{3} \right\rangle = a\left\langle x_{1},x_{2} \right\rangle + b\left\langle x_{1},x_{3} \right\rangle – Linear in 2nd ar­gu­ment
  • ⟨x1,x2⟩=⟨x1,x2⟩\left\langle x_{1},x_{2} \right\rangle = \overline{\left\langle x_{1},x_{2} \right\rangle} – Her­mit­ian Sym­met­ric
  • ⟨x1,x1⟩≥0\left\langle x_{1},x_{1} \right\rangle \geq 0 – Non-neg­a­tive
  • ⟨x1,x1⟩=0\left\langle x_{1},x_{1} \right\rangle = 0, iff x1=0x_{1} = 0 – Pos­i­tive-def­i­nite

A linear space with an inner prod­uct space is called a pre-Hilbert space. For XX, a pre-Hilbert space, we can define a norm on XX to be

‖X‖=⟨x1,x1⟩,\begin{array}{r} \left\| X \right\| = \sqrt{\left\langle x_{1},x_{1} \right\rangle}, \end{array}
(1.1.138)

where XX is also a normed linear space, lead­ing us to define a Hilbert space,

De­f­i­n­i­tion (1.1.139) – A Hilbert space is a com­plete inner prod­uct space.

Com­plete, in layman's terms, means no gaps. In math­e­mat­i­cal analy­sis, a metric space XX is com­plete if every Cauchy se­quence of points ∈X\in X has a limit that is also ∈X\in X, or if every Cauchy se­quence ∈X\in X con­verges ∈X\in X , stated as

Cauchy Sequence: {xj}: limm,n→∞⁡|xm−xn|=0.\begin{array}{r} Cauchy\ Sequence:\ \left\{ x_{j} \right\}:\ \lim_{m,n \rightarrow \infty}{\left| x_{m} - x_{n} \right| = 0}. \end{array}
(1.1.140)

The stan­dard inner prod­uct on ℂn\mathbb{C}^{n}, is given by

⟨x1,x2⟩=∑n=1jx1n∙x2n,\begin{array}{r} \left\langle x_{1},x_{2} \right\rangle = \sum_{n = 1}^{j}{\overline{x_{1_{n}}} \bullet}{x_{2}}_{n}, \end{array}
(1.1.141)

for x1,x2∈ℂx_{1},x_{2}\in \mathbb{C}, which rep­re­sents a finite-di­men­sion­al Hilbert space as it is com­plete. An in­fi­nite-di­men­sion­al Hilbert space ∈ℂ\in \mathbb{C}, where C([a,b])C\left( \lbrack a,b\rbrack \right)is the space of all com­plex-valued func­tions, with a de­fined inner prod­uct

⟨f,g⟩=∫abf(x)g(x) dx,\begin{array}{r} \left\langle f,g \right\rangle = \int_{a}^{b}{f\overline{(x)}}g(x)\ dx, \end{array}
(1.1.142)

for f,g[a,b]:→ℂf,g\lbrack a,b\rbrack: \rightarrow \mathbb{C}, being con­tin­u­ous func­tions. This space may not be com­plete, as it is bound on the in­ter­val [a,b]\lbrack a,b\rbrack, and the com­ple­tion of the space con­cern­ing its norm

‖f‖=(∫ab|f(x)|2 dx)12.\begin{array}{r} \left\| f \right\| = \left( \int_{a}^{b}{\left| f(x) \right|^{2}\ dx} \right)^{\frac{1}{2}}. \end{array}
(1.1.143)

This norm is often de­not­ed as L2([a,b])L^{2}\left( \lbrack a,b\rbrack \right), where the 2-su­per­script­ed rep­re­sents the pp-norm of two for Hilbert spaces. In­fi­nite-di­men­sion­al Hilbert spaces re­quire square-in­te­grable func­tions to be con­ver­gent, and these func­tions would have a finite result for

∫−∞∞|f(x)|2 dx.\begin{array}{r} \int_{- \infty}^{\infty}{\left| f(x) \right|^{2}\ dx.} \end{array}
(1.1.144)

For the circle space, bound on the in­ter­val 𝕋=[0,2π]\mathbb{T} =\lbrack 0,2\pi\rbrack, L2(𝕋)L^{2}\left( \mathbb{T} \right) is Hilbert space of square-in­te­grable func­tions f:𝕋→Cf:\mathbb{T} \rightarrow C, with the inner prod­uct

⟨f,g⟩=∮f(x)g(x) dx.\begin{array}{r} \left\langle f,g \right\rangle = \oint_{}^{}\overline{f(x)}g(x)\ dx. \end{array}
(1.1.145)

Going back to the matrix-vector prod­uct ex­am­ple G(X1(iπ8))∙ϵˆG\left( X_{1}\left( \frac{i\pi}{8} \right) \right) \bullet \widehat{\epsilon}, we denote its result vˆ\widehat{v}, as

G(X1(iπ8))∙ϵˆ=(Acos⁡(π8)+iAsin(π8)Bcos(π8)+iCsin(π8)Ccos(π8)+iBsin(π8))=(A′B′C′)=vˆ,\begin{array}{r} G\left( X_{1}\left( \frac{i\pi}{8} \right) \right) \bullet \widehat{\epsilon} = \begin{pmatrix} A\cos{\left( \frac{\pi}{8} \right) + iAsin\left( \frac{\pi}{8} \right)} \\ Bcos\left( \frac{\pi}{8} \right) + iCsin\left( \frac{\pi}{8} \right) \\ Ccos\left( \frac{\pi}{8} \right) + iBsin\left( \frac{\pi}{8} \right) \end{pmatrix} = \begin{pmatrix} A^{'} \\ B^{'} \\ C^{'} \end{pmatrix} = \widehat{v}, \end{array}
(1.1.146)

and taking the inner prod­uct vˆ∈ℂ3\widehat{v} \in \mathbb{C}^{3} with itself, this is de­ter­mined by eval­u­at­ing

⟨vˆ,vˆ⟩=∑j=13vjvj=|A|2+|B|2+|C|2,\begin{array}{r} \left\langle \widehat{v},\widehat{v} \right\rangle = \sum_{j = 1}^{3}\overline{v_{j}}v_{j} = |A|^{2} + |B|^{2} + |C|^{2}, \end{array}
(1.1.147)

where vj\overline{v_{j}} de­notes the com­plex con­ju­gate of vjv_{j}. The result of ⟨vˆ,vˆ⟩\left\langle \widehat{v},\widehat{v} \right\rangle is pos­i­tive def­i­nite, and the inner prod­uct would only be equal to zero if vˆ=0.\widehat{v} = 0.For ex­am­ple, we take a vector uˆ\widehat{u} equal to the matrix-vector prod­uct G(X1(iπ3))∙ϵˆG\left( X_{1}\left( \frac{i\pi}{3} \right) \right) \bullet \widehat{\epsilon},

G(X1(iπ3))∙ϵˆ=(Acos⁡(π3)+iAsin(π3)Bcos(π3)+iCsin(π3)Ccos(π3)+iBsin(π3))=(A′′B′′C′′)=uˆ,\begin{array}{r} G\left( X_{1}\left( \frac{i\pi}{3} \right) \right) \bullet \widehat{\epsilon} = \begin{pmatrix} A\cos{\left( \frac{\pi}{3} \right) + iAsin\left( \frac{\pi}{3} \right)} \\ Bcos\left( \frac{\pi}{3} \right) + iCsin\left( \frac{\pi}{3} \right) \\ Ccos\left( \frac{\pi}{3} \right) + iBsin\left( \frac{\pi}{3} \right) \end{pmatrix} = \begin{pmatrix} A^{''} \\ B^{''} \\ C^{''} \end{pmatrix} = \widehat{u}, \end{array}
(1.1.148)

and con­sid­er the fol­low­ing inner prod­uct ⟨vˆ,uˆ⟩\left\langle \widehat{v},\widehat{u} \right\rangle

⟨vˆ,uˆ⟩=|A|2e(5iπ24)+(|B|2+|C|2)(cos⁡(π8)cos⁡(π3)+sin⁡(π8)sin⁡(π3))+2iBC(cos⁡(π8)sin⁡(π3)−sin⁡(π8)cos⁡(π3))\begin{array}{r} \left\langle \widehat{v},\widehat{u} \right\rangle = |A|^{2}e^{\left( \frac{5i\pi}{24} \right)} + \left( |B|^{2} + |C|^{2} \right)\left( \cos{\left( \frac{\pi}{8} \right)\cos\left( \frac{\pi}{3} \right)} + \sin{\left( \frac{\pi}{8} \right)\sin\left( \frac{\pi}{3} \right)} \right) \\ + 2iBC\left( \cos{\left( \frac{\pi}{8} \right)\sin\left( \frac{\pi}{3} \right)} - \sin{\left( \frac{\pi}{8} \right)\cos\left( \frac{\pi}{3} \right)} \right) \end{array}
(1.1.149)

which is equal to its con­ju­gate

⟨uˆ,vˆ⟩=|A|2e(5iπ24)+(|B|2+|C|2)(cos⁡(π8)cos⁡(π3)+sin⁡(π8)sin⁡(π3))+2iBC(cos⁡(π8)sin⁡(π3)−sin⁡(π8)cos⁡(π3)).\begin{array}{r} \overline{\left\langle \widehat{u},\widehat{v} \right\rangle} = |A|^{2}e^{\left( \frac{5i\pi}{24} \right)} + \left( |B|^{2} + |C|^{2} \right)\left( \cos{\left( \frac{\pi}{8} \right)\cos\left( \frac{\pi}{3} \right)} + \sin{\left( \frac{\pi}{8} \right)\sin\left( \frac{\pi}{3} \right)} \right) \\ + 2iBC\left( \cos{\left( \frac{\pi}{8} \right)\sin\left( \frac{\pi}{3} \right)} - \sin{\left( \frac{\pi}{8} \right)\cos\left( \frac{\pi}{3} \right)} \right). \end{array}
(1.1.150)

Moving from vec­tors to ma­tri­ces ∈ℂ\in \mathbb{C}, let ℂm×n\mathbb{C}^{m \times n} define the space of all com­plex m×nm \times n ma­tri­ces, and then define the inner prod­uct on ℂm×n\mathbb{C}^{m \times n}

⟨A,B⟩=trace(A∗B),\begin{array}{r} \left\langle A,B \right\rangle = trace\left( A^{*}B \right), \end{array}
(1.1.151)

for which, this inner prod­uct can be rep­re­sent­ed as

⟨A,B⟩=∑j=1m∑k=1naj,kbj,k,\begin{array}{r} \left\langle A,B \right\rangle = \sum_{j = 1}^{m}{\sum_{k = 1}^{n}{{\overline{a}}_{j,k}b_{j,k}}}, \end{array}
(1.1.152)

giving us the same inner prod­uct for a matrix in ℂm×n\mathbb{C}^{m \times n}, with a vector in ℂmn.\mathbb{C}^{mn}.The norm for the cor­re­spond­ing ℂm×n\mathbb{C}^{m \times n} space is called the Hilbert-Schmidt norm, rep­re­sent­ed by

‖A‖=(∑j=1m∑k=1n|aj,k|2)12,\begin{array}{r} \left\| A \right\| = \left( \sum_{j = 1}^{m}{\sum_{k = 1}^{n}\left| a_{j,k} \right|^{2}} \right)^{\frac{1}{2}}, \end{array}
(1.1.153)

which is just the square root of the Frobe­nius norm

‖A‖2=∑j=1m∑k=1n|aj,k|2.\begin{array}{r} \left\| A \right\|^{2} = \sum_{j = 1}^{m}{\sum_{k = 1}^{n}\left| a_{j,k} \right|^{2}}. \end{array}
(1.1.154)

The gen­er­a­tors of U(3)U(3) are com­plex square ma­tri­ces of the form ℂm×n\mathbb{C}^{m \times n}, where m=n.m = n.With the inner prod­uct on ℂm×n\mathbb{C}^{m \times n} de­fined, we en­deav­or to show that gen­er­a­tors and their prod­ucts are an in­fi­nite-di­men­sion­al space on the sphere. To do so, we must show that all the gen­er­a­tors and their prod­ucts are square-in­te­grable func­tions over the in­ter­val [0,2π]\lbrack 0,2\pi\rbrack, verify that our inner prod­uct space is a pre-Hilbert space, and show that it is com­plete.

De­f­i­n­i­tion (1.1.155) – Let 𝒢𝐦([𝟎,𝟐𝛑])\mathcal{G}_{\mathbf{m}}\left( \left\lbrack \mathbf{0},\mathbf{2}\mathbf{\pi} \right\rbrack \right)denote the space of m×mm \times m com­plex matrix func­tion gen­er­a­tors from G(Xj)G\left( X_{j} \right) with an inner prod­uct on 𝒢m\mathcal{G}_{m}, as

⟨G(Xj),G(Xk)⟩=tr(G(Xj)∗G(Xk))\begin{array}{r} \left\langle G\left( X_{j} \right),G\left( X_{k} \right) \right\rangle = tr\left( G\left( X_{j} \right)^{*}G\left( X_{k} \right) \right) \end{array}

for G(Xj),G(Xk):→ℂmG\left( X_{j} \right),G\left( X_{k} \right): \rightarrow \mathbb{C}^{m}, as con­tin­u­ous func­tions, and tr(G(Xj)∗G(Xk))tr\left( G\left( X_{j} \right)^{*}G\left( X_{k} \right) \right) , defin­ing the trace of the prod­uct of the com­plex con­ju­gate trans­form of the ma­tri­ces G(Xj)G\left( X_{j} \right) and G(Xk)G\left( X_{k} \right). The com­ple­tion of 𝒢m\mathcal{G}_{m} con­cern­ing the Hilbert-Schmidt norm

‖G(Xj)‖=(∑j=1m∑k=1n|xj,k|2)12,\begin{array}{r} \left\| G\left( X_{j} \right) \right\| = \left( \sum_{j = 1}^{m}{\sum_{k = 1}^{n}\left| x_{j,k} \right|^{2}} \right)^{\frac{1}{2}}, \end{array}

is the Hilbert Space, de­not­ed as L2(𝕋mℂ)L^{2}\left( \mathbb{T}^{m\mathbb{C}} \right), where 𝕋mℂ\mathbb{T}^{m\mathbb{C}} sig­ni­fies these func­tions are on the sphere (Smℂ)\left( S^{m\mathbb{C}} \right) of di­men­sion mℂm\mathbb{C}.

Lemma (1.1.156) – 𝒢m\mathcal{G}_{m} is a pre-Hilbert space with a de­fined inner prod­uct ⟨G(Xj),G(Xk)⟩=tr⟨G(Xj)∗,G(Xk)⟩\left\langle G\left( X_{j} \right),G\left( X_{k} \right) \right\rangle = tr\left\langle G\left( X_{j} \right)^{*},G\left( X_{k} \right) \right\rangle.

Proof – Ver­i­fi­ca­tion of the inner prod­uct space on 𝒳m\mathcal{X}_{m} holds for all xj∈𝒢mx_{j} \in \mathcal{G}_{m}, and a,b∈ℂa,b\in \mathbb{C}:

  • ⟨x1,x2⟩=⟨x1,x2⟩\left\langle x_{1},x_{2} \right\rangle = \overline{\left\langle x_{1},x_{2} \right\rangle} – Her­mit­ian Sym­met­ric
  • ⟨x1,x1⟩≥0\left\langle x_{1},x_{1} \right\rangle \geq 0 – Non-neg­a­tive
  • ⟨x1,x1⟩=0\left\langle x_{1},x_{1} \right\rangle = 0, iff x1=0x_{1} = 0 – Pos­i­tive-def­i­nite
  • ⟨x1,ax2+bx3⟩=a⟨x1,x2⟩+b⟨x1,x3⟩\left\langle x_{1},ax_{2} + bx_{3} \right\rangle = a\left\langle x_{1},x_{2} \right\rangle + b\left\langle x_{1},x_{3} \right\rangle – Linear in 2nd ar­gu­ment

Taking the inner prod­uct of the gen­er­a­tors G(Xj),G(Xk)G\left( X_{j} \right),G\left( X_{k} \right), we con­firm the inner prod­uct space on 𝒢m\mathcal{G}_{m}:

  • As XjX_{j} are skew-Her­mit­ian ma­tri­ces, we have that eXj=G(Xj)e^{X_{j}} = G\left( X_{j} \right) are uni­tary ma­tri­ces, and as they are uni­tary, G(X1)∗=G(X1)−1.G\left( X_{1} \right)^{*} = G\left( X_{1} \right)^{- 1}. As XjX_{j} are sym­met­ric and uni­tary, we have that G(X1)=G(X1)−1G\overline{\left( X_{1} \right)} = G\left( X_{1} \right)^{- 1} is as well,
∴tr(G(X2)∗G(X1))=tr(G(X2)−1G(X1))−1=tr(G(X1)−1G(X2))=tr(G(X1)∗G(X2)),\therefore\overline{tr\left( G\left( X_{2} \right)^{*}G\left( X_{1} \right) \right)} = tr\left( G\left( X_{2} \right)^{- 1}G\left( X_{1} \right) \right)^{- 1} = tr\left( G\left( X_{1} \right)^{- 1}G\left( X_{2} \right) \right) = tr\left( G\left( X_{1} \right)^{*}G\left( X_{2} \right) \right),

(Her­mit­ian sym­met­ric.)

  • For uni­tary ma­tri­ces, G(X1)∗G(X1)=G(X1)G(X1)∗=Im,G\left( X_{1} \right)^{*}G\left( X_{1} \right) = G\left( X_{1} \right)G\left( X_{1} \right)^{*} = I_{m},
∴tr(G(X1)∗G(X1))=tr(Im)=∑j=1m1=m≠0,\therefore tr\left( G\left( X_{1} \right)^{*}G\left( X_{1} \right) \right) = tr\left( I_{m} \right) = \sum_{j = 1}^{m}{1 = m \neq 0,}

(Non-neg­a­tive and pos­i­tive def­i­nite.)

  • For uni­tary ma­tri­ces and scalars, a,b∈ℂa,b\in \mathbb{C}
⟨G(Xj),aG(Xk)+bG(Xl)⟩=⟨G(Xj),aG(Xk)⟩+⟨G(Xj),bG(Xl)⟩\left\langle G\left( X_{j} \right),aG\left( X_{k} \right) + bG\left( X_{l} \right) \right\rangle = \left\langle G\left( X_{j} \right),aG\left( X_{k} \right) \right\rangle + \left\langle G\left( X_{j} \right),bG\left( X_{l} \right) \right\rangle
=tr(G(Xj)∗aG(Xk))+tr(G(Xj)∗bG(Xl))=tr(aG(Xj)∗G(Xk))+tr(bG(Xj)∗G(Xl))= tr\left( G\left( X_{j} \right)^{*}aG\left( X_{k} \right) \right) + tr\left( G\left( X_{j} \right)^{*}bG\left( X_{l} \right) \right) = tr\left( aG\left( X_{j} \right)^{*}G\left( X_{k} \right) \right) + tr\left( bG\left( X_{j} \right)^{*}G\left( X_{l} \right) \right)
=a∙tr(G(Xj)∗G(Xk))+b∙tr(G(Xj)∗G(Xl))=a⟨G(Xj),G(Xk)⟩+b⟨G(Xj),G(Xl)⟩.= a \bullet tr\left( G\left( X_{j} \right)^{*}G\left( X_{k} \right) \right) + b \bullet tr\left( G\left( X_{j} \right)^{*}G\left( X_{l} \right) \right) = a\left\langle G\left( X_{j} \right),G\left( X_{k} \right) \right\rangle + b\left\langle G\left( X_{j} \right),G\left( X_{l} \right) \right\rangle.

(Lin­ear­i­ty in 2nd ar­gu­ment.)

□\square

Corol­lary (1.1.157) – 𝒳m\mathcal{X}_{m} is a pre-Hilbert space with a de­fined inner prod­uct ⟨Xj,Xk⟩=tr⟨Xj∗,Xk⟩\left\langle X_{j},X_{k} \right\rangle = tr\left\langle X_{j}^{*},X_{k} \right\rangle.

Proof – Ver­i­fi­ca­tion of the inner prod­uct space on 𝒳m\mathcal{X}_{m} holds for all xj∈𝒳mx_{j} \in \mathcal{X}_{m}, and a,b∈ℂa,b\in \mathbb{C}:

  • As XjX_{j} are skew-Her­mit­ian ma­tri­ces and square-sym­met­ric about the main di­ag­o­nal, we have that Xj∗=−Xj,X_{j}^{*} = - X_{j}, Xj∗=(Xj)T=XjT,X_{j}^{*} = \left( \overline{X_{j}} \right)^{T} = \overline{X_{j}^{T}},and XjT=Xj,X_{j}^{T} = X_{j},where­as for the trace we know that tr(XjXk)=tr(XkXj)tr\left( X_{j}X_{k} \right) = tr\left( X_{k}X_{j} \right).
∴(tr(Xj∗Xk)=tr(Xk∗Xj))→(tr(−XjXk)=tr(Xk∗Xj))→(tr(−XjXk)=tr(Xk(−Xj))\therefore\left( tr\left( X_{j}^{*}X_{k} \right) = \overline{tr\left( X_{k}^{*}X_{j} \right)} \right) \rightarrow \left( tr\left( - X_{j}X_{k} \right) = tr\left( \overline{X_{k}^{*}}\overline{X_{j}} \right) \right) \rightarrow \left( tr\left( - X_{j}X_{k} \right) = tr\left( X_{k}( - X_{j} \right) \right)
→(tr(−XjXk)=tr(−XkXj))→(tr(XjXk)=tr(XkXj))→(tr(XjXk)=tr(XjXk)),\rightarrow \left( tr\left( - X_{j}X_{k} \right) = tr\left( {- X}_{k}X_{j} \right) \right) \rightarrow \left( tr\left( X_{j}X_{k} \right) = tr\left( X_{k}X_{j} \right) \right) \rightarrow \left( tr\left( X_{j}X_{k} \right) = tr\left( X_{j}X_{k} \right) \right),

(Her­mit­ian sym­met­ric.)

  • For the skew-Her­mit­ian ma­tri­ces, Xj,X_{j},
tr(Xj∗Xj)=tr(|xj|×Im)=∑j=1m|xj|≠0,tr\left( X_{j}^{*}X_{j} \right) = tr\left( {\left| x_{j} \right| \times I}_{m} \right) = \sum_{j = 1}^{m}{\left| x_{j} \right| \neq 0,}

and tr(Xj∗Xj)=0\ tr\left( X_{j}^{*}X_{j} \right) = 0, only for xj=0.x_{j} = 0.

(Non-neg­a­tive and pos­i­tive def­i­nite.)

  • For skew-Her­mit­ian ma­tri­ces and scalars a,b∈ℂa,b\in \mathbb{C}
⟨Xj,aXk+bXl⟩=⟨Xj,aXk⟩+⟨Xj,bXl⟩\left\langle X_{j},aX_{k} + bX_{l} \right\rangle = \left\langle X_{j},aX_{k} \right\rangle + \left\langle X_{j},bX_{l} \right\rangle
=tr(XjaXk)+tr(XjbXl)=tr(aXj∗Xk)+tr(bXj∗Xl)= tr\left( X_{j}aX_{k} \right) + tr\left( X_{j}bX_{l} \right) = tr\left( aX_{j}^{*}X_{k} \right) + tr\left( {bX}_{j}^{*}X_{l} \right)
=a∙tr(Xj∗Xk)+b∙tr(Xj∗Xl)=a⟨Xj,Xk⟩+b⟨Xj,Xl⟩.= a \bullet tr\left( X_{j}^{*}X_{k} \right) + b \bullet tr\left( X_{j}^{*}X_{l} \right) = a\left\langle X_{j},X_{k} \right\rangle + b\left\langle X_{j}{,X}_{l} \right\rangle.

(Lin­ear­i­ty in 2nd ar­gu­ment.)

□\square

The com­ple­tion of 𝒢m([0,2π])\mathcal{G}_{m}\left( \lbrack 0,2\pi\rbrack \right) as­so­ci­at­ed with the Hilbert-Schmidt norm

‖G(Xj)‖=(∑j=1m∑k=1n|xj,k|2)12=(tr(G(Xj)∗G(Xj)))12,\begin{array}{r} \left\| G\left( X_{j} \right) \right\| = \left( \sum_{j = 1}^{m}{\sum_{k = 1}^{n}\left| x_{j,k} \right|^{2}} \right)^{\frac{1}{2}} = \left( tr\left( G\left( X_{j} \right)^{*}G\left( X_{j} \right) \right) \right)^{\frac{1}{2}}, \end{array}
(1.1.158)

can also be as­so­ci­at­ed with the norm

‖G(Xj)‖=(tr(12π∫02π|G(Xj)|2dxj))12,\begin{array}{r} \left\| G\left( X_{j} \right) \right\| = \left( tr\left( \frac{1}{2\pi}\int_{0}^{2\pi}\left| G\left( X_{j} \right) \right|^{2}dx_{j} \right) \right)^{\frac{1}{2}}, \end{array}
(1.1.159)

for the matrix gen­er­a­tor G(Xj)G\left( X_{j} \right), and for each gen­er­a­tor G(Xj)∈U(3)G\left( X_{j} \right) \in U(3), we find

‖G(Xj)‖=(∑j=1m∑k=1n|xj,k|2)12=(tr(12π∫02π|G(Xj)|2dxj))12=3.\begin{array}{r} \left\| G\left( X_{j} \right) \right\| = \left( \sum_{j = 1}^{m}{\sum_{k = 1}^{n}\left| x_{j,k} \right|^{2}} \right)^{\frac{1}{2}} = \left( tr\left( \frac{1}{2\pi}\int_{0}^{2\pi}\left| G\left( X_{j} \right) \right|^{2}dx_{j} \right) \right)^{\frac{1}{2}} = \sqrt{3}. \end{array}
(1.1.160)

Lemma (1.1.161) – The matrix gen­er­a­tors G(Xj)G\left( X_{j} \right) are a com­plete metric space on the in­ter­val [0,2π]\lbrack 0,2\pi\rbrack, since these matrix func­tions are square-in­te­grable on this in­ter­val

𝒢m([0,2π])⟼ℂ3 square integrable on [0,2π]⟺∫02π|G(Xj)|2dxj<∞.\mathcal{G}_{m}\left( \lbrack 0,2\pi\rbrack \right) \longmapsto \mathbb{C}^{3}\ square\ integrable\ on\ \lbrack 0,2\pi\rbrack \Longleftrightarrow \int_{0}^{2\pi}\left| G\left( X_{j} \right) \right|^{2}dx_{j} < \infty.

Proof – The gen­er­a­tors G(Xj)∈U(3)\ G\left( X_{j} \right) \in U(3) are uni­tary, there­fore |G(Xj)|2=I3\left| G\left( X_{j} \right) \right|^{2} = I_{3}, and when in­te­grat­ed over the in­ter­val [0,2π]\lbrack 0,2\pi\rbrack, con­cern­ing dxjdx_{j}, will be finite for all angles of Xj.X_{j}.

 □ \begin{array}{r} \ \square\ \end{array}

We have now shown the math­e­mat­ics to sup­port De­f­i­n­i­tion (1.1.139), where the com­ple­tion of 𝒢m([0,2π])\mathcal{G}_{m}\left( \lbrack 0,2\pi\rbrack \right) con­cern­ing the Hilbert-Schmidt norm, de­not­ed by L2(𝕋mℂ)L^{2}\left( \mathbb{T}^{m\mathbb{C}} \right), is a Hilbert space. To con­sid­er the in­ter­val of our space 𝒢m([0,2π])\mathcal{G}_{m}\left( \lbrack 0,2\pi\rbrack \right), a de­f­i­n­i­tion of the inner prod­uct of com­plex matrix func­tions on 𝒢m([0,2π])\mathcal{G}_{m}\left( \lbrack 0,2\pi\rbrack \right), that con­sis­tent with the or­thog­o­nal­i­ty re­la­tions of the gen­er­a­tor, is ex­pressed as

⟨G(Xj),G(Xk)⟩=14π2∬02πG(Xj)∗G(Xk) dxk dxj.\begin{array}{r} \left\langle G\left( X_{j} \right),G\left( X_{k} \right) \right\rangle = \frac{1}{4\pi^{2}}\iint_{0}^{2\pi}{G\left( X_{j} \right)^{*}G\left( X_{k} \right)\ dx_{k}\ dx_{j}.} \end{array}
(1.1.162)

We in­te­grate the prod­uct G(Xj)∗G(Xk)G\left( X_{j} \right)^{*}G\left( X_{k} \right) over both xkx_{k} and xj\ x_{j}, on the or­thog­o­nal­ly sit­u­at­ed cir­cles of the Xk/XjX_{k}/X_{j}-planes, having in­ter­vals [0,2π].\lbrack 0,2\pi\rbrack.We would expect ⟨G(Xj),G(Xj)⟩=1\left\langle G\left( X_{j} \right),G\left( X_{j} \right) \right\rangle = 1 for our unit gen­er­a­tors, and j≠kj \neq k ⟨G(Xj),G(Xk)⟩=0\left\langle G\left( X_{j} \right),G\left( X_{k} \right) \right\rangle = 0. The fol­low­ing two ex­am­ples nicely demon­strate this, with the first show­ing the inner prod­uct of G(X1)G\left( X_{1} \right) with itself, yields the iden­ti­ty

⟨G(X1),G(X1)⟩=12π∫02πG(X1)∗G(X1) dx1 \left\langle G\left( X_{1} \right),G\left( X_{1} \right) \right\rangle = \frac{1}{2\pi}\int_{0}^{2\pi}{G\left( X_{1} \right)^{*}G\left( X_{1} \right)\ dx_{1}\ }
=12π∫02π(e−ix1000cos⁡(x1)−isin(x1)0−isin(x1)cos⁡(x1))(eix1000cos⁡(x1)isin(x1)0isin(x1)cos⁡(x1))dx1 = \frac{1}{2\pi}\int_{0}^{2\pi}{\begin{pmatrix} e^{- ix_{1}} & 0 & 0 \\ 0 & \cos\left( x_{1} \right) & - isin\left( x_{1} \right) \\ 0 & - isin\left( x_{1} \right) & \cos\left( x_{1} \right) \end{pmatrix}\begin{pmatrix} e^{ix_{1}} & 0 & 0 \\ 0 & \cos\left( x_{1} \right) & isin\left( x_{1} \right) \\ 0 & isin\left( x_{1} \right) & \cos\left( x_{1} \right) \end{pmatrix}dx_{1}\ }
12π∫02π(100010001)dx1 =12π (x1000x1000x1)|02π=I3,\begin{array}{r} \frac{1}{2\pi}\int_{0}^{2\pi}{\begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix}dx_{1}\ } = {\frac{1}{2\pi}\left. \ \begin{pmatrix} x_{1} & 0 & 0 \\ 0 & x_{1} & 0 \\ 0 & 0 & x_{1} \end{pmatrix} \right|}_{0}^{2\pi} = I_{3}, \end{array}
(1.1.163)

and for the second, we demon­strate the or­thog­o­nal re­la­tion be­tween the func­tions, G(X1),G(X2)G\left( X_{1} \right),G\left( X_{2} \right), where we in­te­grate over both x1x_{1} and x2x_{2}

⟨G(X1),G(X2)⟩=14π2∬02πG(X1)∗G(X2) dx2 dx1\left\langle G\left( X_{1} \right),G\left( X_{2} \right) \right\rangle = \frac{1}{4\pi^{2}}\iint_{0}^{2\pi}{G\left( X_{1} \right)^{*}G\left( X_{2} \right)\ dx_{2}\ dx_{1}}
=14π2∬02π(e−ix1000cos⁡(x1)−isin(x1)0−isin(x1)cos⁡(x1))(cos⁡(x2)0isin(x2)0eix20isin(x2)0cos⁡(x2)) dx2 dx1= \frac{1}{4\pi^{2}}\iint_{0}^{2\pi}{\begin{pmatrix} e^{- ix_{1}} & 0 & 0 \\ 0 & \cos\left( x_{1} \right) & - isin\left( x_{1} \right) \\ 0 & - isin\left( x_{1} \right) & \cos\left( x_{1} \right) \end{pmatrix}\begin{pmatrix} \cos\left( x_{2} \right) & 0 & isin\left( x_{2} \right) \\ 0 & e^{ix_{2}} & 0 \\ isin\left( x_{2} \right) & 0 & \cos\left( x_{2} \right) \end{pmatrix}\ dx_{2}\ dx_{1}}
=14π2∬02π(e−ix1cos⁡(x2)0e−ix1isin(x2)sin⁡(x1)sin⁡(x2)eix2cos⁡(x1)−isin(x1)cos⁡(x2)icos(x1)sin⁡(x2)−isin⁡(x1)eix2cos⁡(x1)cos⁡(x2)) dx2 dx1= \frac{1}{4\pi^{2}}\iint_{0}^{2\pi}{\begin{pmatrix} e^{- ix_{1}}\cos\left( x_{2} \right) & 0 & e^{- ix_{1}}isin\left( x_{2} \right) \\ \sin\left( x_{1} \right)\sin\left( x_{2} \right) & e^{ix_{2}}\cos\left( x_{1} \right) & - isin\left( x_{1} \right)\cos\left( x_{2} \right) \\ icos\left( x_{1} \right)\sin\left( x_{2} \right) & - {i\sin\left( x_{1} \right)e}^{{ix}_{2}} & \cos{\left( x_{1} \right)\cos\left( x_{2} \right)} \end{pmatrix}\ dx_{2}\ dx_{1}}
=14π2  (ie−ix1sin⁡(x2)0e−ix1cos⁡(x2)cos⁡(x1)cos⁡(x2)−ieix2sin⁡(x1)icos(x1)sin⁡(x2)−isin(x1)cos⁡(x2)cos⁡(x1)eix2sin⁡(x1)sin⁡(x2))|02π|02π= \frac{1}{4\pi^{2}}\left. \ \left. \ \begin{pmatrix} ie^{- ix_{1}}\sin\left( x_{2} \right) & 0 & e^{- ix_{1}}\cos\left( x_{2} \right) \\ \cos\left( x_{1} \right)\cos\left( x_{2} \right) & - ie^{ix_{2}}\sin\left( x_{1} \right) & icos\left( x_{1} \right)\sin\left( x_{2} \right) \\ - isin\left( x_{1} \right)\cos\left( x_{2} \right) & \cos\left( x_{1} \right)e^{{ix}_{2}} & \sin{\left( x_{1} \right)\sin\left( x_{2} \right)} \end{pmatrix} \right|_{0}^{2\pi} \right|_{0}^{2\pi}
=14π2(000000000)=0.\begin{array}{r} = \frac{1}{4\pi^{2}}\begin{pmatrix} 0 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix} = 0. \end{array}
(1.1.164)

The inner prod­uct on 𝒢m([0,2π])\mathcal{G}_{m}\left( \lbrack 0,2\pi\rbrack \right) from De­f­i­n­i­tion (1.1.155), can also be ap­plied to vec­tors in the space, seen by set­ting vˆ1,vˆ2{\widehat{v}}_{1},{\widehat{v}}_{2} to be equal to

vˆ1=G(X1)×(1,1,1)T=(eix1,cos⁡(x1)+isin(x1),cos⁡(x1)+isin(x1)),\begin{array}{r} {\widehat{v}}_{1} = G\left( X_{1} \right) \times (1,1,1)^{T} = \left( e^{ix_{1}},\cos\left( x_{1} \right) + isin\left( x_{1} \right),\cos\left( x_{1} \right) + isin\left( x_{1} \right) \right), \end{array}
(1.1.165)
vˆ2=G(X2)×(1,1,1)T=(cos⁡(x2)+isin(x2),eix2,cos⁡(x2)+isin(x2)),\begin{array}{r} {\widehat{v}}_{2} = G\left( X_{2} \right) \times (1,1,1)^{T} = \left( \cos\left( x_{2} \right) + isin\left( x_{2} \right),e^{ix_{2}},\cos\left( x_{2} \right) + isin\left( x_{2} \right) \right), \end{array}
(1.1.166)

we can verify the or­thog­o­nal­i­ty con­di­tion be­tween any vector ∈X1\in X_{1}, with an­oth­er ∈X2\in X_{2}

⟨vˆ1,vˆ2⟩=14π2∬02πvˆ1∗vˆ2 dx2 dx1\left\langle {\widehat{v}}_{1},{\widehat{v}}_{2} \right\rangle = \frac{1}{4\pi^{2}}\iint_{0}^{2\pi}{{\widehat{v}}_{1}^{*}{\widehat{v}}_{2}\ dx_{2}\ dx_{1}}
=94π2∬02πcos⁡(x1−x2)−isin(x1−x2) dx2 dx1= \frac{9}{4\pi^{2}}\iint_{0}^{2\pi}{\cos\left( x_{1} - x_{2} \right) - isin\left( x_{1} - x_{2} \right)}\ dx_{2}\ dx_{1}
=94π2∫02π0 dx1=0.\begin{array}{r} = \frac{9}{4\pi^{2}}\int_{0}^{2\pi}{0\ dx_{1}} = 0. \end{array}
(1.1.167)
  1. Hunter and Nachter­gaele, "Ap­plied Analy­sis", World Sci­en­tif­ic, 2001