Preconditioner: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Manoguru
→‎Jacobi (or diagonal) preconditioner: 'diag' in roman, not italic
Line 1: Line 1:
In [[mathematics]], '''Choi's theorem on completely positive maps''' (after [[Man-Duen Choi]]) is a result that classifies completely positive maps between finite-dimensional (matrix) [[C*-algebra]]s. An infinite-dimensional algebraic generalization of Choi's theorem is known as [[Viacheslav Belavkin|Belavkin]]'s "[[Radon–Nikodym theorem|Radon–Nikodym]]" theorem for completely positive maps.
Her name is Felicidad Ahmad. Managing individuals is what I do in my working day job. To play croquet is some thing that I've done for many years. Her family members lives in Idaho.<br><br>Look at my web-site :: auto warranty [[http://Dev.holesofsa.com/UserProfile/tabid/57/userId/14239/Default.aspx dig this]]
 
== Some preliminary notions ==
 
Before stating Choi's result, we give the definition of a completely positive map and fix some notation.  '''C'''<sup>''n''&nbsp;&times;&nbsp;''n''</sup> will denote the C*-algebra of ''n'' &times; ''n'' complex matrices. We will call ''A'' &isin; '''C'''<sup>''n'' &times; ''n''</sup> '''positive''', or symbolically, ''A''&nbsp;&ge;&nbsp;0, if ''A'' is Hermitian and the [[spectrum of a matrix|spectrum]] of ''A'' is nonnegative. (This condition is also called '''positive semidefinite'''.)
 
A linear map &Phi; : '''C'''<sup>''n'' &times; ''n''</sup> &rarr; '''C'''<sup>''m'' &times; ''m''</sup> is said to be a '''positive map''' if &Phi;(''A'') &ge; 0 for all ''A'' &ge; 0. In other words, a map Φ is positive if it preserves Hermiticity and the cone of positive elements.
 
Any linear map Φ induces another map
 
:<math>I_k \otimes \Phi : \mathbb{C} ^{k \times k} \otimes \mathbb{C} ^{n \times n} \rightarrow \mathbb{C} ^{k \times k} \otimes \mathbb{C} ^{m \times m}</math>
 
in a natural way: define
 
:<math>
( I_k \otimes \Phi ) (M \otimes A) = M \otimes \Phi (A)
</math>
 
and extend by linearity. In matrix notation, a general element in
 
:<math>\mathbb{C} ^{k \times k} \otimes \mathbb{C} ^{n \times n}</math>
 
can be expressed as a ''k'' &times; ''k'' operator matrix:
 
:<math>
\begin{bmatrix}
A_{11} & \cdots & A_{1k} \\
\vdots & \ddots & \vdots \\
A_{k1} & \cdots & A_{kk}
\end{bmatrix},
</math>
 
and its image under the induced map is
 
:<math>
(I_k \otimes \Phi)
(\begin{bmatrix} A_{11} & \cdots & A_{1k} \\ \vdots & \ddots & \vdots \\A_{k1} & \cdots & A_{kk} \end{bmatrix})
=
\begin{bmatrix}
\Phi (A_{11}) & \cdots & \Phi( A_{1k} ) \\
\vdots & \ddots & \vdots \\
\Phi (A_{k1}) & \cdots & \Phi( A_{kk} )
\end{bmatrix}.
</math>
 
Writing out the individual elements in the above matrix-of-matrices amounts to the natural identification of algebras
:<math>
\mathbb{C}^{k\times k}\otimes\mathbb{C}^{m\times m}\cong\mathbb{C}^{km\times km}.
</math>
We say that Φ is '''k-positive''' if <math>I_k \otimes \Phi</math>, considered as an element of '''C'''<sup>''km''&times;''km''</sup>, is a positive map, and Φ is called '''completely positive'''
if Φ is k-positive for all k.
 
The [[Transpose|transposition map]] is a standard example of a positive map that fails to be 2-positive. Let T denote this map on '''C''' <sup>2 &times; 2</sup>. The following is a positive matrix in <math> \mathbb{C} ^{2 \times 2} \otimes \mathbb{C}^{2 \times 2} </math>:
 
:<math>
\begin{bmatrix}
\begin{pmatrix}1&0\\0&0\end{pmatrix}&
\begin{pmatrix}0&1\\0&0\end{pmatrix}\\
\begin{pmatrix}0&0\\1&0\end{pmatrix}&
\begin{pmatrix}0&0\\0&1\end{pmatrix}
\end{bmatrix}
=
\begin{bmatrix}
1 & 0 & 0 & 1 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 \\
1 & 0 & 0 & 1 \\
\end{bmatrix} .
</math>
 
The image of this matrix under <math>I_2 \otimes T</math> is
 
:<math>
\begin{bmatrix}
\begin{pmatrix}1&0\\0&0\end{pmatrix}^T&
\begin{pmatrix}0&1\\0&0\end{pmatrix}^T\\
\begin{pmatrix}0&0\\1&0\end{pmatrix}^T&
\begin{pmatrix}0&0\\0&1\end{pmatrix}^T
\end{bmatrix}
=
\begin{bmatrix}
1 & 0 & 0 & 0 \\
0 & 0 & 1 & 0 \\
0 & 1 & 0 & 0 \\
0 & 0 & 0 & 1 \\
\end{bmatrix} ,
</math>
 
which is clearly not positive, having determinant -1.
 
Incidentally, a map Φ is said to be '''co-positive''' if the composition &Phi;  <math>\circ</math> ''T'' is positive. The transposition map itself is a co-positive map.
 
The above notions concerning positive maps extend naturally to maps between C*-algebras.
 
== Choi's result ==
 
=== Statement of theorem ===
Choi's theorem reads as follows:
 
Let
 
:<math>\Phi : \mathbb{C} ^{n \times n} \rightarrow \mathbb{C} ^{m \times m}</math>
 
be a positive map. The following are equivalent:
 
i) <math>\Phi</math> is ''n''-positive.
 
ii) The matrix with operator entries
 
:<math>\; C_\Phi=(I_n\otimes\Phi)(\sum_{ij}E_{ij}\otimes E_{ij}) = \sum_{ij}E_{ij}\otimes\Phi(E_{ij}) \in \mathbb{C} ^{nm \times nm}</math>
 
is positive, where <math>E_{ij}\in\mathbb{C}^{n\times n}</math> is the matrix with 1 in the <math>ij</math>-th entry and 0s elsewhere. (The matrix <math>C_\Phi</math> is sometimes called the ''Choi matrix'' of <math>\Phi</math>.)
 
iii) <math>\Phi</math> is completely positive.
 
=== Proof ===
 
To show i) implies ii), we observe that if
:<math>
\; E=\sum_{ij}E_{ij}\otimes E_{ij},
</math>
then ''E''=''E''<sup>*</sup> and ''E''<sup>2</sup>=''nE'', so ''E''=''n''<sup>−1</sup>''EE''<sup>*</sup> which is positive and
''C''<sub>&Phi;</sub>=(''I<sub>n</sub>''&otimes;&Phi;)(''E'') is positive by the ''n''-positivity of &Phi;.
 
If iii) holds, then so does i) trivially.
 
We now turn to the argument for ii) &rArr; iii). This mainly involves chasing the different ways of looking at '''C'''<sup>''nm''&times;''nm''</sup>:
:<math>
    \mathbb{C}^{nm\times nm}
\cong\mathbb{C}^{nm}\otimes(\mathbb{C}^{nm})^*
\cong\mathbb{C}^n\otimes\mathbb{C}^m\otimes(\mathbb{C}^n\otimes\mathbb{C}^m)^*
\cong\mathbb{C}^n\otimes(\mathbb{C}^n)^*\otimes\mathbb{C}^m\otimes(\mathbb{C}^m)^*
\cong\mathbb{C}^{n\times n}\otimes\mathbb{C}^{m\times m}.
</math>
 
Let the eigenvector decomposition of ''C''<sub>&Phi;</sub> be
 
:<math>
\; C_\Phi = \sum _{i = 1} ^{nm} \lambda_i v_i v_i ^* ,
</math>
 
where the vectors <math>v_i</math> lie in '''C'''<sup>''nm''</sup> . By assumption, each eigenvalue <math>\lambda_i</math> is non-negative so we can absorb the eigenvalues in the eigenvectors and redefine <math>v_i</math> so that
 
:<math>
\; C_\Phi = \sum _{i = 1} ^{nm} v_i v_i ^* .
</math>
 
The vector space '''C'''<sup>''nm''</sup> can be viewed as the direct sum <math>\textstyle \oplus_{i=1}^n \mathbb{C}^m</math>
compatibly with the above identification <math>\textstyle\mathbb{C}^{nm}\cong\mathbb{C}^n\otimes\mathbb{C}^m</math>
and the standard basis of '''C'''<sup>''n''</sup>.
 
If ''P<sub>k</sub>'' &isin; '''C'''<sup>''m'' &times; ''nm''</sup> is projection onto the ''k''-th copy of '''C'''<sup>''m''</sup>, then ''P<sub>k</sub>''<sup>*</sup> &isin; '''C'''<sup>''nm''&times;''m''</sup>
is the inclusion of '''C'''<sup>''m''</sup> as the ''k''-th summand of the direct sum and
 
:<math>
\; \Phi (E_{kl}) = P_k \cdot C_\Phi \cdot P_l^* = \sum _{i = 1} ^{nm} P_k v_i ( P_l v_i )^*.
</math>
 
Now if the operators ''V<sub>i</sub>'' &isin; '''C'''<sup>''m''&times;''n''</sup> are defined on the ''k''-th standard
basis vector ''e<sub>k</sub>'' of '''C'''<sup>''n''</sup> by
 
:<math>\; V_i e_k = P_k v_i,</math>
 
then
 
:<math>
\; \Phi (E_{kl}) = \sum _{i = 1} ^{nm} P_k v_i ( P_l v_i )^* = \sum _{i = 1} ^{nm} V_i e_k e_l ^*  V_i ^*
= \sum _{i = 1} ^{nm} V_i E_{kl}  V_i ^*.
</math>
 
Extending by linearity gives us
 
:<math>
\; \Phi(A) = \sum_{i=1}^{nm} V_i A V_i^*
</math>
for any ''A'' &isin; '''C'''<sup>''n'' &times; ''n''</sup>. Since any map of this form is manifestly completely positive, we have the desired result.
 
The above is essentially Choi's original proof. Alternative proofs have also been known.
 
== Consequences ==
 
=== Kraus operators ===
 
In the context of [[quantum information theory]], the operators {''V<sub>i</sub>''} are called the ''[[Kraus operator]]s'' (after [[Karl Kraus (physicist)|Karl Kraus]]) of Φ. Notice, given a completely positive Φ, its Kraus operators need not be unique. For example, any "square root" factorization of the Choi matrix
:<math>
\; C_\Phi = B^* B.
</math>
 
gives a set of Kraus operators. (Notice ''B'' need not be the unique positive [[square root of a matrix|square root]] of the Choi matrix.)
 
Let
 
:<math>B^* = [b_1, \ldots, b_{nm}] ,</math>
 
where ''b''<sub>i</sub>*'s are the row vectors of ''B'', then
 
:<math>
\; C_\Phi = \sum _{i = 1} ^{nm} b_i b_i ^*.  
</math>
 
The corresponding Kraus operators can be obtained by exactly the same argument from the proof.
 
When the Kraus operators are obtained from the eigenvector decomposition of the Choi matrix, because the eigenvectors form an orthogonal set, the corresponding Kraus operators are also orthogonal in the [[Hilbert–Schmidt operator|Hilbert–Schmidt]] [[inner product]]. This is not true in general for Kraus operators obtained from square root factorizations. (Positive semidefinite matrices do not generally have a unique square-root factorizations.)
 
If two sets of Kraus operators {''A<sub>i</sub>''}<sub>1</sub><sup>''nm''</sup> and {''B<sub>i</sub>''}<sub>1</sub><sup>''nm''</sup> represent the same completely positive map Φ, then there exists a unitary ''operator'' matrix
 
:<math>\{U_{ij}\}_{ij} \in \mathbb{C}^{nm^2 \times nm^2} \quad \text{such that}
\quad A_i = \sum _{i = 1} U_{ij} B_j.
</math>
 
This can be viewed as a special case of the result relating two [[Stinespring factorization theorem|minimal Stinespring representations]].
 
Alternatively, there is an isometry ''scalar'' matrix {''u<sub>ij</sub>''}<sub>''ij''</sub> &isin; '''C'''<sup>''nm'' &times; ''nm''</sup> such that
 
:<math>
\; A_i = \sum _{i = 1} u_{ij} B_j.
</math>
 
This follows from the fact that for two square matrices ''M'' and ''N'', ''M M*'' = ''N N*'' if and only if ''M = N U'' for some unitary ''U''.
 
=== Completely copositive maps ===
 
It follows immediately from Choi's theorem that Φ is completely copositive if and only if it is of the form
 
:<math>\Phi(A) = \sum _i V_i A^T V_i ^* .</math>
 
=== Hermitian-preserving maps ===
 
Choi's technique can be used to obtain a similar result for a more general class of maps. &Phi; is said to be Hermitian-preserving if ''A'' is Hermitian implies &Phi;(''A'') is also Hermitian. One can show &Phi; is Hermitian-preserving if and only if it is of the form
 
:<math>\Phi (A) = \sum_{i=1} ^{nm} \lambda_i V_i A V_i ^*</math>
 
where &lambda;<sub>''i''</sub> are real numbers, the eigenvalues of ''C''<sub>&Phi;</sub>,
and each ''V''<sub>''i''</sub> corresponds to an eigenvector of ''C''<sub>&Phi;</sub>. Unlike the completely positive case, ''C''<sub>&Phi;</sub> may fail to be positive. Since Hermitian matrices do not admit factorizations of the form ''B*B'' in general, the Kraus representation is no longer possible for a given Φ.
 
== See also ==
 
*[[Stinespring factorization theorem]]
*[[Quantum operation]]
*[[Holevo's theorem]]
 
== References ==
 
* M. Choi, ''Completely Positive Linear Maps on Complex matrices'', Linear Algebra and Its Applications, 285–290, 1975
 
* V. P. Belavkin, P. Staszewski, ''Radon-Nikodym Theorem for Completely Positive Maps,'' Reports on Mathematical Physics, v.24, No 1, 49–55, 1986.
* J. de Pillis, ''Linear Transformations Which Preserve Hermitian and Positive Semidefinite Operators'', Pacific Journal of Mathematics, 129–137, 1967.
 
[[Category:Linear algebra]]
[[Category:Operator theory]]
[[Category:Articles containing proofs]]
[[Category:Theorems in functional analysis]]

Revision as of 17:33, 24 February 2014

Her name is Felicidad Ahmad. Managing individuals is what I do in my working day job. To play croquet is some thing that I've done for many years. Her family members lives in Idaho.

Look at my web-site :: auto warranty [dig this]