Physlib

PhyslibAlpha.ProbabilisticTheory.CStarAlgebra.Stinespring.Kernel

Completely positive maps as positive kernels

i. Overview

A completely positive map `J` into the bounded operators on `H` is a positive operator-valued kernel: applied to the Gram matrix of finitely many elements of `A`, it gives a positive block operator on `Hⁿ`. A Stinespring witness of `J` is a representation `π` on a Hilbert space `K` and an operator `V : H → K` with `J a = V⋆ π(a) V`.

ii. Key results

  • `blockMatrixMap` : a matrix of operators acting on `Hⁿ`.
  • `cpKernel_inner_nonneg_natural'` : completely positive maps give positive kernels.
  • `StinespringWitness` : a Stinespring dilation of `J`.
  • `completelyPositiveMap_map_star_general` : completely positive maps preserve the adjoint.

iii. Table of contents

  • A. `blockMatrixMap`: the finite block-operator representation
  • B. The CP-map positivity kernel
  • C. Stinespring witnesses

A. `blockMatrixMap`: the finite block-operator representation

This part uses only the codomain Hilbert space `H`, no algebra `A` at all.

B. The CP-map positivity kernel

From here on, `A` is a unital C⋆-algebra with the compatible Mathlib order (`[CStarAlgebra A] [PartialOrder A] [StarOrderedRing A]`) — the field content of `OperatorAlgebra A`.

C. Stinespring witnesses

21 declarations

definition

Bounded linear operator on HnH^n induced by a block matrix MM of operators on HH

Given a Hilbert space HH and a natural number nn, let MM be an n×nn \times n matrix whose entries MijM_{ij} are bounded linear operators on HH (i.e., Mij∈B(H)M_{ij} \in \mathcal{B}(H)). The function `blockMatrixMap` defines a bounded linear operator on the direct sum space Hn≅⨁k=1nHH^n \cong \bigoplus_{k=1}^n H such that for any vector x=(x1,…,xn)∈Hn\mathbf{x} = (x_1, \dots, x_n) \in H^n, the ii-th component of the output is given by ∑j=1nMijxj\sum_{j=1}^n M_{ij} x_j.

theorem

(blockMatrixMap M x)i=∑jMijxj(\text{blockMatrixMap} \, M \, \mathbf{x})_i = \sum_j M_{ij} x_j

Let HH be a Hilbert space and nn be a natural number. Let MM be an n×nn \times n matrix of bounded linear operators on HH, where Mij∈B(H)M_{ij} \in \mathcal{B}(H) denotes the entry in the ii-th row and jj-th column. For any vector x=(x1,…,xn)\mathbf{x} = (x_1, \dots, x_n) in the direct sum space Hn≅⨁k=1nHH^n \cong \bigoplus_{k=1}^n H, the ii-th component of the vector produced by the operator blockMatrixMap(M)\text{blockMatrixMap}(M) acting on x\mathbf{x} is given by the sum: ((blockMatrixMap M)x)i=∑j=1nMijxj ((\text{blockMatrixMap} \, M) \mathbf{x})_i = \sum_{j=1}^n M_{ij} x_j where (⋅)i(\cdot)_i denotes the projection onto the ii-th coordinate.

theorem

blockMatrixMap(M∗)=(blockMatrixMap(M))∗\text{blockMatrixMap}(M^*) = (\text{blockMatrixMap}(M))^*

Let HH be a Hilbert space and nn be a natural number. Let MM be an n×nn \times n matrix whose entries MijM_{ij} are bounded linear operators on HH (i.e., Mij∈B(H)M_{ij} \in \mathcal{B}(H)). Let M∗M^* denote the adjoint matrix, whose entries are given by (M∗)ij=(Mji)∗(M^*)_{ij} = (M_{ji})^*, where (Mji)∗(M_{ji})^* is the adjoint operator of MjiM_{ji}. Then the bounded linear operator on the direct sum space Hn≅⨁k=1nHH^n \cong \bigoplus_{k=1}^n H induced by the matrix M∗M^* is equal to the adjoint of the bounded linear operator induced by MM. That is, blockMatrixMap(M∗)=(blockMatrixMap(M))∗\text{blockMatrixMap}(M^*) = (\text{blockMatrixMap}(M))^*

theorem

blockMatrixMap(M⋅N)=blockMatrixMap(M)∘blockMatrixMap(N)\text{blockMatrixMap}(M \cdot N) = \text{blockMatrixMap}(M) \circ \text{blockMatrixMap}(N)

Let HH be a Hilbert space and nn be a natural number. For any two n×nn \times n matrices MM and NN whose entries are bounded linear operators on HH (i.e., Mij,Nij∈B(H)M_{ij}, N_{ij} \in \mathcal{B}(H)), the bounded linear operator on the direct sum space HnH^n induced by the matrix product M⋅NM \cdot N is equal to the composition of the bounded linear operators induced by MM and NN individually. That is, blockMatrixMap(M⋅N)=blockMatrixMap(M)∘blockMatrixMap(N)\text{blockMatrixMap}(M \cdot N) = \text{blockMatrixMap}(M) \circ \text{blockMatrixMap}(N) where blockMatrixMap\text{blockMatrixMap} is the map from operator-valued matrices to operators on HnH^n.

theorem

blockMatrixMap(1)=id\text{blockMatrixMap}(1) = \text{id}

Let HH be a Hilbert space and nn be a natural number. Let 11 denote the n×nn \times n identity matrix whose entries are bounded linear operators on HH (where the diagonal entries are the identity operator idH\text{id}_H and the off-diagonal entries are the zero operator). The bounded linear operator on the direct sum space Hn≅⨁k=1nHH^n \cong \bigoplus_{k=1}^n H induced by this matrix is the identity operator on HnH^n. That is, blockMatrixMap(1)=idHn\text{blockMatrixMap}(1) = \text{id}_{H^n} where blockMatrixMap\text{blockMatrixMap} is the map from operator-valued matrices to operators on HnH^n.

theorem

blockMatrixMap(M+N)=blockMatrixMap(M)+blockMatrixMap(N)\text{blockMatrixMap}(M + N) = \text{blockMatrixMap}(M) + \text{blockMatrixMap}(N)

Let HH be a Hilbert space and nn be a natural number. For any two n×nn \times n matrices MM and NN whose entries are bounded linear operators on HH (i.e., Mij,Nij∈B(H)M_{ij}, N_{ij} \in \mathcal{B}(H)), the bounded linear operator on the direct sum space Hn≅⨁k=1nHH^n \cong \bigoplus_{k=1}^n H induced by the matrix sum M+NM + N is equal to the sum of the bounded linear operators induced by MM and NN individually. That is, blockMatrixMap(M+N)=blockMatrixMap(M)+blockMatrixMap(N)\text{blockMatrixMap}(M + N) = \text{blockMatrixMap}(M) + \text{blockMatrixMap}(N) where blockMatrixMap\text{blockMatrixMap} is the map from operator-valued matrices to operators on HnH^n.

theorem

`blockMatrixMap` of the zero matrix is the zero operator

For any natural number nn and Hilbert space HH, the bounded linear operator on Hn≅⨁k=1nHH^n \cong \bigoplus_{k=1}^n H induced by the n×nn \times n zero matrix of operators (where each entry is the zero operator on HH) is the zero operator on HnH^n.

definition

Star-algebra representation of operator-valued matrices on HnH^n

For a natural number nn and a complex Hilbert space HH, `blockMatrixRepresentation` is the unital C\mathbb{C}-star-algebra homomorphism from the algebra of n×nn \times n matrices with entries in B(H)\mathcal{B}(H) (the bounded linear operators on HH) to the algebra of bounded linear operators on the direct sum Hilbert space Hn≅⨁k=1nHH^n \cong \bigoplus_{k=1}^n H. This homomorphism maps an operator-valued matrix MM to the operator that acts on a vector x=(x1,…,xn)∈Hn\mathbf{x} = (x_1, \dots, x_n) \in H^n such that its ii-th component is ∑j=1nMijxj\sum_{j=1}^n M_{ij} x_j. It preserves the additive, multiplicative, scalar multiplication, identity, and adjoint (star) structures of the matrix algebra.

theorem

M≥0M \ge 0 implies `blockMatrixMap` MM is a positive operator

Let HH be a complex Hilbert space and nn be a natural number. Let MM be an n×nn \times n matrix where each entry MijM_{ij} is a bounded linear operator on HH. If MM is a positive operator-valued matrix (M≥0M \ge 0), then the bounded linear operator on the direct sum Hilbert space Hn≅⨁k=1nHH^n \cong \bigoplus_{k=1}^n H induced by MM (via `blockMatrixMap`) is a positive operator.

theorem

∑i,j⟨xi,Mijxj⟩≥0\sum_{i,j} \langle x_i, M_{ij} x_j \rangle \ge 0 for Positive Operator Matrices

Let HH be a complex Hilbert space and let MM be an n×nn \times n matrix where each entry MijM_{ij} is a bounded linear operator on HH. If MM is a positive operator-valued matrix (M≥0M \ge 0), then for any vector x=(x1,…,xn)x = (x_1, \dots, x_n) in the direct sum Hilbert space Hn=⨁i=1nHH^n = \bigoplus_{i=1}^n H, the following inequality holds: ∑i=1n∑j=1n⟨xi,Mijxj⟩≥0\sum_{i=1}^n \sum_{j=1}^n \langle x_i, M_{ij} x_j \rangle \ge 0 where ⟨⋅,⋅⟩\langle \cdot, \cdot \rangle denotes the inner product on HH.

theorem

Completely positive maps are star-preserving: J(a∗)=J(a)∗J(a^*) = J(a)^*

Let AA and BB be C∗C^*-algebras. For any completely positive map J:A→BJ : A \to B and any element a∈Aa \in A, the map JJ preserves the star operation, satisfying J(a∗)=J(a)∗J(a^*) = J(a)^*.

definition

Row matrix of aa with non-zero zeroth row

Given a sequence of nn elements a=(a0,a1,…,an−1)a = (a_0, a_1, \dots, a_{n-1}) in an algebra AA (typically a C∗C^*-algebra), where nn is a positive integer, the function defines an n×nn \times n matrix MM whose entries MijM_{ij} are given by aja_j if i=0i = 0 and 00 otherwise. In other words, MM is a square matrix where the zeroth row is the vector aa and all other rows consist of zero elements.

definition

Gram matrix Gij=ai∗ajG_{ij} = a_i^* a_j for a∈Ana \in A^n

Given a positive integer nn and a sequence of nn elements a=(a0,a1,…,an−1)a = (a_0, a_1, \dots, a_{n-1}) in a C∗C^*-algebra AA, the Gram matrix of aa is the n×nn \times n matrix over AA whose entry at row ii and column jj is defined as ai∗aja_i^* a_j, where ai∗a_i^* denotes the adjoint of aia_i.

theorem

The (i,j)(i, j)-th entry of the Gram matrix is ai∗aja_i^* a_j

For a positive integer nn and a sequence of elements a=(a0,a1,…,an−1)a = (a_0, a_1, \dots, a_{n-1}) in a C∗C^*-algebra AA, the entry of the Gram matrix of aa at row ii and column jj is given by ai∗aja_i^* a_j, where ai∗a_i^* denotes the adjoint of aia_i.

theorem

The Gram matrix Gij=ai∗ajG_{ij} = a_i^* a_j is positive

Let AA be a unital C∗C^*-algebra. For any positive integer nn and any sequence of elements a0,a1,…,an−1∈Aa_0, a_1, \dots, a_{n-1} \in A, the Gram matrix G∈Mn(A)G \in M_n(A) with entries Gij=ai∗ajG_{ij} = a_i^* a_j (where ai∗a_i^* denotes the adjoint of aia_i) is positive, i.e., G≥0G \ge 0.

theorem

A Completely Positive Map applied to a Gram Matrix is Positive

Let AA be a unital C∗C^*-algebra and HH be a complex Hilbert space. Let J:A→B(H)J: A \to \mathcal{B}(H) be a completely positive map, where B(H)\mathcal{B}(H) denotes the space of bounded linear operators on HH. For any positive integer nn and any sequence of elements a0,a1,…,an−1∈Aa_0, a_1, \dots, a_{n-1} \in A, let GG be the n×nn \times n Gram matrix with entries Gij=ai∗ajG_{ij} = a_i^* a_j. Then the matrix J(G)J(G), defined by applying JJ to each entry of GG such that (J(G))ij=J(ai∗aj)(J(G))_{ij} = J(a_i^* a_j), is a positive operator on the direct sum Hilbert space HnH^n.

theorem

Positivity of the kernel of a completely positive map

Let AA be a unital C∗C^*-algebra and HH be a complex Hilbert space. For any completely positive map J:A→B(H)J: A \to \mathcal{B}(H), any finite sequence of elements a1,…,an∈Aa_1, \dots, a_n \in A (where n>0n > 0), and any sequence of vectors x1,…,xn∈Hx_1, \dots, x_n \in H, the following inequality holds: ∑i=1n∑j=1n⟨J(ai∗aj)xj,xi⟩≥0\sum_{i=1}^n \sum_{j=1}^n \langle J(a_i^* a_j) x_j, x_i \rangle \ge 0 where B(H)\mathcal{B}(H) denotes the space of bounded linear operators on HH, ai∗a_i^* is the adjoint of aia_i, and ⟨⋅,⋅⟩\langle \cdot, \cdot \rangle is the inner product on HH.

theorem

Positivity of the Kernel for Completely Positive Maps

Let AA be a unital C∗C^*-algebra and HH be a complex Hilbert space. Let J:A→B(H)J: A \to \mathcal{B}(H) be a completely positive map, where B(H)\mathcal{B}(H) denotes the space of bounded linear operators on HH. For any n∈Nn \in \mathbb{N}, any finite collection of elements a1,…,an∈Aa_1, \dots, a_n \in A, and any vectors x1,…,xn∈Hx_1, \dots, x_n \in H, the following inequality holds: ∑i,j=1n⟨J(ai∗aj)xj,xi⟩≥0\sum_{i, j=1}^n \langle J(a_i^* a_j) x_j, x_i \rangle \ge 0 where ⟨⋅,⋅⟩\langle \cdot, \cdot \rangle denotes the inner product on HH and ai∗a_i^* is the adjoint of aia_i in AA.

theorem

Positivity of ∑i,j⟨xi,J(ai∗aj)xj⟩\sum_{i,j} \langle x_i, J(a_i^* a_j) x_j \rangle for a completely positive map JJ

Let AA be a unital C∗C^*-algebra and HH be a complex Hilbert space. Let J:A→B(H)J : A \to \mathcal{B}(H) be a completely positive map, where B(H)\mathcal{B}(H) denotes the algebra of bounded linear operators on HH. For any positive integer nn, any sequence of elements a0,…,an−1∈Aa_0, \dots, a_{n-1} \in A, and any sequence of vectors x0,…,xn−1∈Hx_0, \dots, x_{n-1} \in H, the following inequality holds: 0≤∑i,j=0n−1⟨xi,J(ai∗aj)xj⟩ 0 \leq \sum_{i,j=0}^{n-1} \langle x_i, J(a_i^* a_j) x_j \rangle where ⟨⋅,⋅⟩\langle \cdot, \cdot \rangle denotes the inner product on HH, which is linear in the second argument.

theorem

Positivity of the kernel ∑i,j⟨xi,J(ai∗aj)xj⟩≥0\sum_{i,j} \langle x_i, J(a_i^* a_j) x_j \rangle \ge 0 for completely positive maps JJ

Let AA be a unital C∗C^*-algebra and HH be a complex Hilbert space. Let B(H)\mathcal{B}(H) denote the space of bounded linear operators on HH. If J:A→B(H)J: A \to \mathcal{B}(H) is a completely positive map, then for any natural number nn, any collection of elements a0,…,an−1∈Aa_0, \dots, a_{n-1} \in A, and any collection of vectors x0,…,xn−1∈Hx_0, \dots, x_{n-1} \in H, the following inequality holds: ∑i,j=0n−1⟨xi,J(ai∗aj)xj⟩≥0\sum_{i,j=0}^{n-1} \langle x_i, J(a_i^* a_j) x_j \rangle \ge 0 where ⟨⋅,⋅⟩\langle \cdot, \cdot \rangle denotes the inner product on HH.

theorem

Stinespring Dilation Identity J(a)=V∗π(a)VJ(a) = V^* \pi(a) V for Completely Positive Maps

Let AA be a unital C∗C^*-algebra and HH be a complex Hilbert space. Let J:A→B(H)J : A \to \mathcal{B}(H) be a completely positive map, where B(H)\mathcal{B}(H) denotes the algebra of bounded linear operators on HH. Given a Stinespring witness WW for JJ—consisting of a representation π:A→B(K)\pi : A \to \mathcal{B}(K) on a Hilbert space KK and a bounded linear operator V:H→KV : H \to K—then for any a∈Aa \in A, the following equality holds: J(a)=V∗π(a)V J(a) = V^* \pi(a) V where V∗V^* is the adjoint operator of VV.