Source author record

Chi-Kwong Li

Chi-Kwong Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

56works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

56 published item(s)

preprint2026arXiv

Cliques and independent subgroups of the Birkhoff polytope graph

The Birkhoff polytope $Ω_n$ is the polytope of doubly stochastic matrices of order $n$. The Birkhoff polytope graph $G(Ω_n)$ is the skeleton of $Ω_n$; it is the Cayley graph whose vertex set consists of the elements of the symmetric group ${\rm Sym}(n)$ of degree $n$, where two permutations are adjacent if one equals the product of the other with a cycle. We study the combinatorial structure of this graph, focusing on its maximal and maximum cliques and on its independent subgroups (subgroups of ${\rm Sym}(n)$ whose elements are pairwise nonadjacent in the graph). We obtain maximal subgroups of $G(Ω_n)$ and establish both a lower bound and an upper bound for its clique number. Especially, we prove that if $K$ is a subset of ${\rm Sym}(n)$ consisting of 3-cycle permutations such that $δ_1^{-1}δ_2$ is a single cycle for all $δ_1,δ_2\in K$, then the maximum size of $K$ is $\lfloor (n-1)^2/4\rfloor$, which can be viewed as an Erdős-Ko-Rado-type theorem for ${\rm Sym}(n)$.

preprint2026arXiv

Semigroup automorphisms of total positivity

Totally positive (TP) and totally nonnegative (TN) matrices connect to analysis, mechanics, and to dual canonical bases in reductive groups, by well-known works of Schoenberg, Gantmacher-Krein, Lusztig, and others. TP matrices form a multiplicatively closed semigroup, contained in the larger monoid of invertible totally nonnegative (ITN) matrices. Whitney and Berenstein-Fomin-Zelevinsky found bidiagonal factorizations of all $n\times n$ ITN and TP matrices into multiplicative generators; a natural question now is to classify the multiplicative automorphisms of these semigroups. In this article, we classify all automorphisms of these semigroups of ITN and TP matrices. In particular, we show that the automorphisms are the same, and they respect the multiplicative generators.

preprint2022arXiv

Commuting normal operators and joint numerical range

Let ${\mathcal H}$ be a complex Hilbert space and let ${\mathcal B}({\mathcal H})$ be the algebra of all bounded linear operators on ${\mathcal H}$. For a positive integer $k$ less than the dimension of ${\mathcal H}$ and ${\mathbf A} = (A_1, \dots, A_m)\in {\mathcal B}({\mathcal H})^m$, the joint $k$-numerical range $W_k({\mathbf A})$ is the set of $(α_1, \dots, α_m) \in{\mathbb C}^m$ such that $α_i = \sum_{j = 1}^k \langle A_ix_j, x_j\rangle$ for an orthonormal set $\{x_1, \ldots, x_k\}$ in ${\mathcal H}$. Relations between the geometric properties of $W_k({\mathbf A})$ and the algebraic and analytic properties of $A_1, \dots, A_m$ are studied. It is shown that there is $k\in {\mathbb N}$ such that $W_k({\mathbf A})$ is a polyhedral set, i.e., the convex hull of a finite set, if and only if $A_1, \dots, A_k$ have a common reducing subspace ${\mathbf V}$ of finite dimension such that the compression of $A_1, \dots, A_m$ on the subspace ${\mathbf V}$ are diagonal operators $D_1, \dots, D_m$ and $W_k({\mathbf A}) = W_k(D_1, \dots, D_m)$. Characterization is also given to ${\bf A}$ such that the closure of $W_k({\mathbf A})$ is polyhedral. The conditions are related to the joint essential numerical range of ${\mathbf A}$. These results are used to study ${\bf A}$ such that (a) $\{A_1, \dots, A_m\}$ is a commuting family of normal operators, or (b) $W_k(A_1, \dots, A_m)$ is polyhedral for every positive integer $k$. It is shown that conditions (a) and (b) are equivalent for finite rank operators but it is no longer true for compact operators. Characterizations are given for compact operators $A_1, \dots, A_m$ satisfying (a) and (b), respectively. Results are also obtained for general non-compact operators.

preprint2022arXiv

The joint $k$-numerical range of operators

Let ${\mathcal B}({\mathcal H})$ be the algebra of all bounded linear operators on the Hilbert space ${\mathcal H}$. For a positive integer $k$ less than the dimension of ${\mathcal H}$ and ${\mathbf A} = (A_1, \dots, A_m)\in {\mathcal B}({\mathcal H})^m$, the joint $k$-numerical range $W_k({\mathbf A})$ is the set of vector $(α_1, \dots, α_m) \in{\mathbb C}^m$ such that $α_i = \sum_{j = 1}^k \langle A_ix_j, x_j\rangle$ for an orthonormal set $\{x_1, \ldots, x_k\}$ in ${\mathcal H}$. Geometrical properties of $W_k({\mathbf A})$ and their relations with the algebraic properties of $\{A_1, \dots, A_m\}$ are investigated in this paper. For example, conditions for $W_k({\mathbf A})$ to be convex are studied. Descriptions are given for the closure of $W_k({\mathbf A})$ and the closure of ${\rm conv}\, W_k({\mathbf A})$ in terms of the joint essential numerical range of ${\mathbf A}$ for infinite dimensional operators $A_1, \dots, A_m$. Characterizations are obtained for $W_k({\mathbf A})$ or ${\rm conv}\, W_k({\mathbf A})$ to be closed. It is shown that $W_k({\mathbf A})$ is a polyhedral set if and only if $A_1, \dots, A_k$ have a common reducing subspace ${\mathbf V}$ of finite dimension such that the compression of $A_1, \dots, A_m$ on the subspace ${\mathbf V}$ are diagonal operators $D_1, \dots, D_m$ and $W_k({\mathbf A}) = W_k(D_1, \dots, D_m)$. Similar results are obtained for ${\bf A}$ such that the closure of $W_k({\mathbf A})$ is polyhedral. Classifications are given for operators satisfying (1) $\{A_1, \dots, A_m\}$ is a commuting family of normal operators, or (2) $W_k(A_1, \dots, A_m)$ is polyhedral for every positive integer $k$ less than $\dim {\mathcal H}$.

preprint2022arXiv

Unitarily invariant Norms on Operators

Let $f$ be a symmetric norm on ${\mathbb R}^n$ and let ${\mathcal B}({\mathcal H})$ be the set of all bounded linear operators on a Hilbert space ${\mathcal H}$ of dimension at least $n$. Define a norm on ${\mathcal B}({\mathcal H})$ by $\|A\|_f = f(s_1(A), \dots, s_n(A))$, where $s_k(A) = \inf\{\|A-X\|: X\in {\mathcal B}({\mathcal H}) \hbox{ has rank less than } k\}$ is the $k$th singular value of $A$. Basic properties of the norm $\|\cdot\|_f$ are obtained including some norm inequalities and characterization of the equality case. Geometric properties of the unit ball of the norm are obtained; the results are used to determine the structure of maps $L$ satisfying $\|L(A)-L(B)\|_f=\|A - B\|_f$ for any $A, B \in {\mathcal B}({\mathcal H})$.

preprint2020arXiv

A Note on Parallel Distinguishability of two Quantum Operations

We consider a homogeneous system of linear equations of the form $A_α^{\otimes N} {\bf x} = 0$ arising from the distinguishability of two quantum operations by $N$ uses in parallel, where the coefficient matrix $A_α$ depends on a real parameter $α$. It was conjectured by Duan et al. that the system has a non-trivial nonnegative solution if and only if $α$ lies in a certain interval $R_N$ depending on $N$. We affirm the necessity part of the conjecture and establish the sufficiency of the conjecture for $N\leq 10$ by presenting explicit non-trivial nonnegative solutions for the linear system.

preprint2020arXiv

Construction of quantum states with special properties by projection methods

We use projection methods to construct (global) quantum states with prescribed reduced (marginal) states, and possibly with some special properties such as having specific eigenvalues, having specific rank and extreme von Neumann or Renyi entropy. Using convex analysis, optimization techniques on matrix manifolds, we obtain algorithms to solve the problem. Matlab programs are written based on these algorithms and numerical examples are illustrated. The numerical results reveal new patterns leading to new insights and research problems on the topic.

preprint2020arXiv

Error correction schemes for fully correlated quantum channels protecting both quantum and classical information

We study efficient quantum error correction schemes for the fully correlated channel on an $n$-qubit system with error operators that assume the form $σ_x^{\otimes n}$, $σ_y^{\otimes n}$, $σ_z^{\otimes n}$. Previous schemes are improved to facilitate implementation. In particular, when $n$ is odd and equals $2k+1$, we describe a quantum error correction scheme using one arbitrary qubit $σ$ to protect the data state $ρ$ in a $2k$-qubit system. The encoding operation $σ\otimes ρ\mapsto Φ(σ\otimes ρ)$ only requires $3k$ CNOT gates (each with one control bit and one target bit). After the encoded state $Φ(σ\otimes ρ)$ goes through the channel, we can apply the inverse operation $Φ^{-1}$ to produce $\tilde σ\otimes ρ$ so that a partial trace operation can recover $ρ$. When $n$ is even and equals $2k+2$, we describe a hybrid quantum error correction scheme using any one of the two classical bits $σ\in \{|ij\rangle \langle ij|: i, j \in \{0,1\}\}$ to protect a $2k$-qubit state $ρ$ and 2 classical bits. The encoding operation $σ\otimes ρ\mapsto Φ(σ\otimes ρ)$ can be done by $3k+2$ CNOT gates and a single quibt Hadamard gate. After the encoded state $Φ(σ\otimes ρ)$ goes through the channel, we can apply the inverse operation $Φ^{-1}$ to produce $σ\otimes ρ$ so that a perfect protection of the two classical bits $σ$ and the $2k$-qubit state is achieved. If one uses an arbitrary $2$-qubit state $σ$, the same scheme will protect $2k$-qubit states. The scheme was implemented using Matlab, Mathematica, Python, and the IBM's quantum computing framework qiskit.

preprint2020arXiv

Joint numerical ranges and communtativity of matrices

The connection between the commutativity of a family of $n\times n$ matrices and the generalized joint numerical ranges is studied. For instance, it is shown that ${\cal F}$ is a family of mutually commuting normal matrices if and only if the joint numerical range $W_k(A_1, \dots, A_m)$ is a polyhedral set for some $k$ satisfying $|n/2-k|\le 1$, where $\{A_1, \dots, A_m\}$ is a basis for the linear span of the family; equivalently, $W_k(X,Y)$ is polyhedral for any two $X, Y \in {\cal F}$. More generally, characterization is given for the $c$-numerical range $W_c(A_1, \dots, A_m)$ to be polyhedral for any $n\times n$ matrices $A_1, \dots, A_m$. Other results connecting the geometrical properties of the joint numerical ranges and the algebraic properties of the matrices are obtained. Implications of the results to representation theory, and quantum information science are discussed.

preprint2020arXiv

On the mixed-unitary rank of quantum channels

In the theory of quantum information, the mixed-unitary quantum channels, for any positive integer dimension $n$, are those linear maps that can be expressed as a convex combination of conjugations by $n\times n$ complex unitary matrices. We consider the mixed-unitary rank of any such channel, which is the minimum number of distinct unitary conjugations required for an expression of this form. We identify several new relationships between the mixed-unitary rank~$N$ and the Choi rank~$r$ of mixed-unitary channels, the Choi rank being equal to the minimum number of nonzero terms required for a Kraus representation of that channel. Most notably, we prove that the inequality $N\leq r^2-r+1$ is satisfied for every mixed-unitary channel (as is the equality $N=2$ when $r=2$), and we exhibit the first known examples of mixed-unitary channels for which $N>r$. Specifically, we prove that there exist mixed-unitary channels having Choi rank $d+1$ and mixed-unitary rank $2d$ for infinitely many positive integers $d$, including every prime power $d$. We also examine the mixed-unitary ranks of the mixed-unitary Werner--Holevo channels.

preprint2016arXiv

Bounds on probability of state transfer with respect to readout time and edge weight

We analyse the sensitivity of a spin chain modelled by an undirected weighted connected graph exhibiting perfect state transfer to small perturbations in readout time and edge weight in order to obtain physically relevant bounds on the probability of state transfer. At the heart of our analysis is the concept of the numerical range of a matrix; our analysis of edge weight errors additionally makes use of the spectral and Frobenius norms.

preprint2016arXiv

Inequalities on generalized matrix functions

We prove inequalities on non-integer powers of products of generalized matrices functions on the sum of positive semi-definite matrices. For example, for any real number $r \in \{1\} \cup [2, \infty)$, positive semi-definite matrices $A_i,\ B_i,\ C_i\in M_{n_i}$, $i=1,2$, and generalized matrix functions $d_χ, d_ξ$ such as the determinant and permanent, etc., we have \begin{eqnarray*}&&\left(d_χ(A_1+B_1+C_1)d_ξ(A_2+B_2+C_2)\right)^r \\ &&\hskip 1in + \left(d_χ(A_1)d_ξ(A_2)\right)^r + \left(d_χ(B_1)d_ξ(B_2)\right)^r + \left(d_χ(C_1)d_ξ(C_2)\right)^r \\ & \ge &\left(d_χ(A_1+B_1 )d_ξ(A_2+B_2 )\right)^r + \left(d_χ(A_1+ C_1)d_ξ(A_2+ C_2)\right)^r + \left(d_χ( B_1+C_1)d_ξ( B_2+C_2)\right)^r\,.\end{eqnarray*} A general scheme is introduced to prove more general inequalities involving $m$ positive semi-definite matrices for $m \ge 3$ that extend the results of other authors.

preprint2016arXiv

Numerical Ranges of the product of Operators

We study containment regions of the numerical range of the product of operators $A$ and $B$ such that $W(A)$ and $W(B)$ are line segments. It is shown that the containment region is equal to the convex hull of elliptical disks determined by the spectrum of $AB$, and conditions on $A$ and $B$ for the set equality holding are obtained. The results cover the case when $A$ and $B$ are self-adjoint operators extending the previous results on the numerical range of the product of two orthogonal projections.

preprint2016arXiv

Optimal Bounds on Functions of Quantum States under Quantum Channels

Let $ρ_1, ρ_2$ be quantum states and $(ρ_1,ρ_2) \mapsto D(ρ_1, ρ_2)$ be a scalar function such as the trace norm, the fidelity, and the relative entropy, etc. We determine optimal bounds for $D(ρ_1, Φ(ρ_2))$ for $Φ\in \mathcal{S}$ for different class of functions $D(\cdot, \cdot)$, where $\mathcal{S}$ is the set of unitary quantum channels, the set of mixed unitary channels, the set of unital quantum channels, and the set of all quantum channels.

preprint2016arXiv

Quantifying the coherence of pure quantum states

In recent years, several measures have been proposed for characterizing the coherence of a given quantum state. We derive several results that illuminate how these measures behave when restricted to pure states. Notably, we present an explicit characterization of the closest incoherent state to a given pure state under the trace distance measure of coherence. We then use this result to show that the states maximizing the trace distance of coherence are exactly the maximally coherent states. We define the trace distance of entanglement and show that it coincides with the trace distance of coherence for pure states. Finally, we give an alternate proof to a recent result that the $\ell_1$ measure of coherence of a pure state is never smaller than its relative entropy of coherence.

preprint2015arXiv

Discontinuity of Maximum Entropy Inference and Quantum Phase Transitions

In this paper, we discuss the connection between two genuinely quantum phenomena --- the discontinuity of quantum maximum entropy inference and quantum phase transitions at zero temperature. It is shown that the discontinuity of the maximum entropy inference of local observable measurements signals the non-local type of transitions, where local density matrices of the ground state change smoothly at the transition point. We then propose to use the quantum conditional mutual information of the ground state as an indicator to detect the discontinuity and the non-local type of quantum phase transitions in the thermodynamic limit.

preprint2015arXiv

Maximal noiseless code rates for collective rotation channels on qudits

We study noiseless subsystems on collective rotation channels of qudits, i.e., quantum channels with operators in the set ${\mathcal E}(d,n) = \{ U^{\otimes n}: U \in {\mathrm{SU}}(d)\}.$ This is done by analyzing the decomposition of the algebra ${\mathcal A}(d,n)$ generated by ${\mathcal E}(d,n)$. We summarize the results for the channels on qubits ($d=2$), and obtain the maximum dimension of the noiseless subsystem that can be used as the quantum error correction code for the channel. Then we extend our results to general $d$. In particular, it is shown that the code rate, i.e., the number of protected qudits over the number of physical qudits, always approaches 1 for a suitable noiseless subsystem. Moreover, one can determine the maximum dimension of the noiseless subsystem by solving a non-trivial discrete optimization problem. The maximum dimension of the noiseless subsystem for $d = 3$ (qutrits) is explicitly determined by a combination of mathematical analysis and the symbolic software Mathematica.

preprint2015arXiv

Minkowski product of convex sets and product numerical range

Let $K_1, K_2$ be two compact convex sets in $\mathit{C}$. Their Minkowski product is the set $K_1K_2 = \{ab: a \in K_1, b\in K_2\}$. We show that the set $K_1K_2$ is star-shaped if $K_1$ is a line segment or a circular disk. Examples for $K_1$ and $K_2$ are given so that $K_1$ and $K_2$ are triangles (including interior) and $K_1K_2$ is not star-shaped. This gives a negative answer to a conjecture by Puchala et. al concerning the product numerical range in the study of quantum information science. Additional results and open problems are presented.

preprint2015arXiv

Product of positive semi-definite matrices

It is known that every complex square matrix with nonnegative determinant is the product of positive semi-definite matrices. There are characterizations of matrices that require two or five positive semi-definite matrices in the product. However, the characterizations of matrices that require three or four positive semi-definite matrices in the product are lacking. In this paper, we give a complete characterization of these two types of matrices. With these results, we give an algorithm to determine whether a square matrix can be expressed as the product of $k$ positive semi-definite matrices but not fewer, for $k = 1,2,3,4,5$.

preprint2014arXiv

Conditions for degradability of tripartite quantum states

Alice, Bob, and Eve share a pure quantum state. We introduce the notion of state degradability by asking whether the joint density of Alice and Eve can be transformed to the joint density of Alice and Bob by processing Eve's part through a quantum channel, in order words, degrading Eve. We prove necessary and sufficient conditions for state degradability and provide an efficient method to quickly rule out degradability for a given state. The problem of determining degradability of states is different from that of quantum channels, although the notion is similar. One application of state degradability is that it can be used to test channel degradability. In particular, the degradability of the output state of a channel obtained from the maximally entangled input state gives information about the degradability of the channel.

preprint2014arXiv

Factoring a quadratic operator as a product of two positive contractions

Let $T$ be a quadratic operator on a complex Hilbert space $H$. We show that $T$ can be written as a product of two positive contractions if and only if $T$ is of the form $$aI \oplus bI \oplus\begin{pmatrix} aI & P \cr 0 & bI \cr \end{pmatrix} \quad \text{on} \quad H_1\oplus H_2\oplus (H_3\oplus H_3)$$ for some $a, b\in [0,1]$ and strictly positive operator $P$ with $\|P\| \le |\sqrt{a} - \sqrt{b}|\sqrt{(1-a)(1-b)}.$ Also, we give a necessary condition for a bounded linear operator $T$ with operator matrix $\begin{pmatrix} T_1 & T_3\\ 0 & T_2\cr\end{pmatrix}$ on $H\oplus K$ that can be written as a product of two positive contractions.

preprint2014arXiv

Perturbing eigenvalues of non-negative matrices

Let $A$ be an irreducible (entrywise) nonnegative $n\times n$ matrix with eigenvalues $$ρ, b+ic,b-ic, λ_4,\cdots,λ_n,$$ where $ρ$ is the Perron eigenvalue. It is shown that for any $t \in [0, \infty)$ there is a nonnegative matrix with eigenvalues $$ρ+ \tilde t,λ_2+t,λ_3+t, λ_4 \cdots,λ_n,$$ whenever $\tilde t \ge γ_n t$ with $γ_3=1, γ_4 = 2, γ_5=\sqrt 5$ and $γ_n = 2.25$ for $n \ge 6$. The result improves that of Guo et al. Our proof depends on an auxiliary result in geometry asserting that the area of an $n$-sided convex polygon is bounded by $γ_n$ times the maximum area of the triangle lying inside the polygon.

preprint2014arXiv

Positivity of Partitioned Hermitian Matrices with Unitarily Invariant Norms

We give a short proof of a recent result of Drury on the positivity of a $3\times 3$ matrix of the form $(\|R_i^*R_j\|_{\rm tr})_{1 \le i, j \le 3}$ for any rectangular complex (or real) matrices $R_1, R_2, R_3$ so that the multiplication $R_i^*R_j$ is compatible for all $i, j$, where $\|\cdot\|_{\rm tr}$ denotes the trace norm. We then give a complete analysis of the problem when the trace norm is replaced by other unitarily invariant norms.

preprint2014arXiv

Preservers of Unitary Similarity Functions on Lie Products of Matrices

Denote by $M_n$ the set of $n\times n$ complex matrices. Let $f: M_n \rightarrow [0,\infty)$ be a continuous map such that $f(μUAU^*)= f(A)$ for any complex unit $μ$, $A \in M_n$ and unitary $U \in M_n$, $f(X)=0$ if and only if $X=0$ and the induced map $t \mapsto f(tX)$ is monotonic increasing on $[0,\infty)$ for any rank 1 nilpotent $X \in M_n$. Characterizations are given for surjective maps $ϕ$ on $M_n$ satisfying $f(AB-BA) = f(ϕ(A)ϕ(B)-ϕ(B)ϕ(A))$. The general theorem are then used to deduce results on special cases when the function is the pseudo spectrum and the pseudo spectral radius, that answers a question of Molnar raised at the 2014 CMS summer meeting.

preprint2014arXiv

Projection methods in quantum information science

We consider the problem of constructing quantum operations or channels, if they exist, that transform a given set of quantum states $\{ρ_1, \dots, ρ_k\}$ to another such set $\{\hatρ_1, \dots, \hatρ_k\}$. In other words, we must find a {\em completely positive linear map}, if it exists, that maps a given set of density matrices to another given set of density matrices. This problem, in turn, is an instance of a positive semi-definite feasibility problem, but with highly structured constraints. The nature of the constraints makes projection based algorithms very appealing when the number of variables is huge and standard interior point-methods for semi-definite programming are not applicable. We provide emperical evidence to this effect. We moreover present heuristics for finding both high rank and low rank solutions. Our experiments are based on the \emph{method of alternating projections} and the \emph{Douglas-Rachford} reflection method.

preprint2014arXiv

Ranks and eigenvalues of states with prescribed reduced states

For a quantum state represented as an $n\times n$ density matrix $σ\in M_n$, let ${\cal S}(σ)$ be the compact convex set of quantum states $ρ= (ρ_{ij}) \in M_{m\cdot n}$ with the first partial trace equal to $σ$, i.e., ${\rm tr}_1(ρ) = ρ_{11} + \cdots + ρ_{mm} = σ$. It is known that if $m \ge n$ then there is a rank one matrix $ρ\in {\cal S}(σ)$ satisfying ${\rm tr}_1(ρ) = σ$. If $m < n$, there may not be rank one matrix in ${\cal S}(σ)$. In this paper, we determine the ranks of the elements and ranks of the extreme points of the set ${\cal S}$ We also determine $ρ^* \in {\cal S}(σ)$ with rank bounded by $k$ such that $\|{\rm tr}_1ρ^*) - σ\|$ is minimum for a given unitary similarity invariant norm $\|\cdot\|$. Furthermore, the relation between the eigenvalues of $σ$ and those of $ρ\in {\cal S}(σ)$ is analyzed. Extension of our results and open problems will be mentioned.

preprint2014arXiv

Recursive encoding and decoding of the noiseless subsystem for qudits

We give a full explanation of the noiseless subsystem that protects a single-qubit against collective errors and the corresponding recursive scheme described by C.-K. Li et. al. [Phys. Rev. A 84, 044301 (2011)] from a representation theory point of view. Furthermore, we extend the construction to qudits under the influence of collective SU($d$) errors. We find that under this recursive scheme, the asymptotic encoding rate is $1/d$.

preprint2014arXiv

Solution to time-energy costs of quantum channels

We derive a formula for the time-energy costs of general quantum channels proposed in [Phys. Rev. A 88, 012307 (2013)]. This formula allows us to numerically find the time-energy cost of any quantum channel using positive semidefinite programming. We also derive a lower bound to the time-energy cost for any channels and the exact the time-energy cost for a class of channels which includes the qudit depolarizing channels and projector channels as special cases.

preprint2014arXiv

The spectrum of the product of operators, and the product of their numerical ranges

We show that a compact operator $A$ is a multiple of a positive semi-definite operator if and only if $$ σ(AB) \subseteq \overline{W(A)W(B)}, \quad\text{for all (rank one) operators $B$}. $$ An example of a normal operator is given to show that the equivalence conditions may fail in general. We then obtain conditions to identify other classes of operators $A$ so that equivalence conditions hold.

preprint2013arXiv

Decomposition of quantum gates

A recurrence scheme is presented to decompose an $n$-qubit unitary gate to the product of no more than $N(N-1)/2$ single qubit gates with small number of controls, where $N = 2^n$. Detailed description of the recurrence steps and formulas for the number of $k$-controlled single qubit gates in the decomposition are given. Comparison of the result to a previous scheme is presented, and future research directions are discussed.

preprint2013arXiv

Decomposition of unitary matrices and quantum gates

A general scheme is presented to decompose a $d$-by-$d$ unitary matrix as the product of two-level unitary matrices with additional structure and prescribed determinants. In particular, the decomposition can be done by using two-level matrices in $d-1$ classes, where each class is isomorphic to the group of $2\times 2$ unitary matrices. The proposed scheme is easy to apply, and useful in treating problems with the additional structural restrictions. A Matlab program is written to implement the scheme, and the result is used to deduce the fact that every quantum gate acting on $n$-qubit registers can be expressed as no more than $2^{n-1}(2^n-1)$ fully controlled single-qubit gates chosen from $2^n-1$ classes, where the quantum gates in each class share the same $n-1$ control qubits. Moreover, it is shown that it is easy to adjust the proposed decomposition scheme to take advantage of additional structure evolving in the process.

preprint2013arXiv

Determinantal and eigenvalue inequalities for matrices with numerical ranges in a sector

Let $A = \pmatrix A_{11} & A_{12} \cr A_{21} & A_{22}\cr\pmatrix \in M_n$, where $A_{11} \in M_m$ with $m \le n/2$, be such that the numerical range of $A$ lies in the set $\{e^{iφ} z \in \IC: |\Im z| \le (\Re z) \tan α\}$, for some $φ\in [0, 2π)$ and $α\in [0, π/2)$. We obtain the optimal containment region for the generalized eigenvalue $λ$ satisfying $$λ\pmatrix A_{11} & 0 \cr 0 & A_{22}\cr\pmatrix x = \pmatrix 0 & A_{12} \cr A_{21} & 0\cr\pmatrix x \quad \hbox{for some nonzero} x \in \IC^n,$$ and the optimal eigenvalue containment region of the matrix $I_m - A_{11}^{-1}A_{12} A_{22}^{-1}A_{21}$ in case $A_{11}$ and $A_{22}$ are invertible. From this result, one can show $|\det(A)| \le \sec^{2m}(α) |\det(A_{11})\det(A_{22})|$. In particular, if $A$ is a accretive-dissipative matrix, then $|\det(A)| \le 2^m |\det(A_{11})\det(A_{22})|$. These affirm some conjectures of Drury and Lin.

preprint2013arXiv

Linear maps preserving Ky Fan norms and Schatten norms of tensor products of matrices

For a positive integer $n$, let $M_n$ be the set of $n\times n$ complex matrices. Suppose $\|\cdot\|$ is the Ky Fan $k$-norm with $1 \le k \le mn$ or the Schatten $p$-norm with $1 \le p \le \infty$ ($p\ne 2$) on $M_{mn}$, where $m,n\ge 2$ are positive integers. It is shown that a linear map $ϕ: M_{mn} \rightarrow M_{mn}$ satisfying $$\|A\otimes B\| = \|ϕ(A\otimes B)\| \quad \hbox{for all} A \in M_m \hbox{and} B \in M_n$$ if and only if there are unitary $U, V \in M_{mn}$ such that $ϕ$ has the form $A\otimes B \mapsto U(φ_1(A) \otimes φ_2(B))V$, where $φ_s(X)$ is either the identity map $X \mapsto X$ or the transposition map $X \mapsto X^t$. The results are extended to tensor space $M_{n_1} \otimes \cdots \otimes M_{n_m}$ of higher level. The connection of the problem to quantum information science is mentioned.

preprint2012arXiv

A geometric characterization of invertible quantum measurement maps

A geometric characterization is given for invertible quantum measurement maps. Denote by ${\mathcal S}(H)$ the convex set of all states (i.e., trace-1 positive operators) on Hilbert space $H$ with dim$H\leq \infty$, and $[ρ_1, ρ_2]$ the line segment joining two elements $ρ_1, ρ_2$ in ${\mathcal S}(H)$. It is shown that a bijective map $ϕ:{\mathcal S}(H) \rightarrow {\mathcal S}(H)$ satisfies $ϕ([ρ_1, ρ_2]) \subseteq [ϕ(ρ_1),ϕ(ρ_2)]$ for any $ρ_1, ρ_2 \in {\mathcal S}$ if and only if $ϕ$ has one of the following forms $$ρ\mapsto \frac{MρM^*}{{\rm tr}(MρM^*)}\quad \hbox{or} \quad ρ\mapsto \frac{Mρ^T M^*}{{\rm tr}(Mρ^T M^*)},$$ where $M$ is an invertible bounded linear operator and $ρ^T$ is the transpose of $ρ$ with respect to an arbitrarily fixed orthonormal basis.

preprint2012arXiv

Entanglement transformation between sets of bipartite pure quantum states using local operations

Alice and Bob are given an unknown initial state chosen from a set of pure quantum states. Their task is to transform the initial state to a corresponding final pure state using local operations only. We prove necessary and sufficient conditions on the existence of such a transformation. We also provide efficient algorithms that can quickly rule out the possibility of transforming a set of initial states to a set of final states.

preprint2012arXiv

Linear preservers and quantum information science

Let $m,n\ge 2$ be positive integers, $M_m$ the set of $m\times m$ complex matrices and $M_n$ the set of $n\times n$ complex matrices. Regard $M_{mn}$ as the tensor space $M_m\otimes M_n$. Suppose $|\cdot|$ is the Ky Fan $k$-norm with $1 \le k \le mn$, or the Schatten $p$-norm with $1 \le p \le \infty$ ($p\ne 2$) on $M_{mn}$. It is shown that a linear map $ϕ: M_{mn} \rightarrow M_{mn}$ satisfying $$|A\otimes B| = |ϕ(A\otimes B)|$$ for all $A \in M_m$ and $B \in M_n$ if and only if there are unitary $U, V \in M_{mn}$ such that $ϕ$ has the form $A\otimes B \mapsto U(φ_1(A) \otimes φ_2(B))V$, where $φ_i(X)$ is either the identity map $X \mapsto X$ or the transposition map $X \mapsto X^t$. The results are extended to tensor space $M_{n_1} \otimes ... \otimes M_{n_m}$ of higher level. The connection of the problem to quantum information science is mentioned.

preprint2012arXiv

Physical transformations between quantum states

Given two sets of quantum states {A_1, ..., A_k} and {B_1, ..., B_k}, represented as sets of density matrices, necessary and sufficient conditions are obtained for the existence of a physical transformation T, represented as a trace-preserving completely positive map, such that T(A_i) = B_i for i = 1, ..., k. General completely positive maps without the trace-preserving requirement, and unital completely positive maps transforming the states are also considered.

preprint2011arXiv

A note on the realignment criterion

For a quantum state in a bipartite system represented as a density matrix, researchers used the realignment matrix and functions on its singular values to study the separability of the quantum state. We obtain bounds for elementary symmetric functions of singular values of realignment matrices. This answers some open problems proposed by Lupo, Aniello, and Scardicchio. As a consequence, we show that the proposed scheme by these authors for testing separability would not work if the two subsystems of the bipartite system have the same dimension.

preprint2011arXiv

Efficient Quantum Error Correction for Fully Correlated Noise

We investigate an efficient quantum error correction of a fully correlated noise. Suppose the noise is characterized by a quantum channel whose error operators take fully correlated forms given by $σ_x^{\otimes n}$, $σ_y^{\otimes n}$ and $σ_z^{\otimes n}$, where $n>2$ is the number of qubits encoding the codeword. It is proved that (i) $n$ qubits codeword encodes $(n-1)$ data qubits when $n$ is odd and (ii) $n$ qubits codeword implements a noiseless subsystem encoding $(n-2)$ data qubits when $n$ is even. Quantum circuits implementing these schemes are constructed.

preprint2011arXiv

Recovery in quantum error correction for general noise without measurement

It is known that one can do quantum error correction without syndrome measurement, which is often done in operator quantum error correction (OQEC). However, the physical realization could be challenging, especially when the recovery process involves high-rank projection operators and a superoperator. We use operator theory to improve OQEC so that the implementation can always be done by unitary gates followed by a partial trace operation. Examples are given to show that our error correction scheme outperforms the existing ones in various scenarios.

preprint2011arXiv

Recursive Encoding and Decoding of Noiseless Subsystem and Decoherence Free Subspace

When the environmental disturbace to a quantum system has a wavelength much larger than the system size, all qubits localized within a small area are under action of the same error operators. Noiseless subsystem and decoherence free subspace are known to correct such collective errors. We construct simple quantum circuits, which implement these collective error correction codes, for a small number $n$ of physical qubits. A single logical qubit is encoded with $n=3$ and $n=4$, while two logical qubits are encoded with $n=5$. The recursive relations among the subspaces employed in noiseless subsystem and decoherence free subspace play essential rôles in our implementation. The recursive relations also show that the number of gates required to encode $m$ logical qubits increases linearly in $m$.

preprint2011arXiv

The automorphism group of separable states in quantum information theory

We show that the linear group of automorphism of Hermitian matrices which preserves the set of separable states is generated by \emph{natural} automorphisms: change of an orthonormal basis in each tensor factor, partial transpose in each tensor factor, and interchanging two tensor factors of the same dimension. We apply our results to preservers of the product numerical range.

preprint2010arXiv

"Maps preserving the spectrum of generalized Jordan product of operators", and its "Addendum"

In the paper "Maps preserving the spectrum of generalized Jordan product of operators", we define a generalized Jordan products on standard operator algebras $A_1, A_2$ on complex Banach spaces $X_1, X_2$, respectively. This includes the usual Jordan product $A_1 \circ A_2 = A_1 A_2 + A_2 A_1$, and the triple $\{A_1,A_2,A_3\} = A_1 A_2 A_3 + A_3 A_2 A_1$. Let a map $Φ: A_1 \to A_2$ prserving the spectra of the products $$ σ(Φ(A_1) \circ ... \circ Φ(A_k)) = σ(A_1\circ ... \circ A_k) $$ whenever any one of $A_1, ..., A_k$ has rank at most one. It is shown in this paper that if the range of $Φ$ contains all operators of rank at most three, then $Φ$ must be a Jordan isomorphism multiplied by an $m$th root of unity. Similar results for maps between self-adjoint operators acting on Hilbert spaces are also obtained. After our paper "Maps preserving the spectrum of generalized Jordan product of operators" was published in Linear Algebra Appl. 432 (2010), 1049-1069, Jianlian Cui pointed out that some arguments in the proof of Theorem 3.1 are not entirely clear and accurate. Here we supply some details in the "Addendum".

preprint2010arXiv

Evolution of unconditional dispersal in periodic environments

Organisms modulate their fitness in heterogeneous environments by dispersing. Prior work shows that there is selection against "unconditional" dispersal in spatially heterogeneous environments. "Unconditional" means individuals disperse at a rate independent of their location. We prove that if within-patch fitness varies spatially and between two values temporally, then there is selection for unconditional dispersal: any evolutionarily stable strategy (ESS) or evolutionarily stable coalition (ESC) includes a dispersive phenotype. Moreover, at this ESS or ESC, there is at least one sink patch (i.e. geometric mean of fitness less than one) and no sources patches (i.e. geometric mean of fitness greater than one). These results coupled with simulations suggest that spatial-temporal heterogeneity due to abiotic forcing result in either an ESS with a dispersive phenotype or an ESC with sedentary and dispersive phenotypes. In contrast, spatial-temporal heterogeneity due to biotic interactions can select for higher dispersal rates that ultimately spatially synchronize population dynamics.

preprint2010arXiv

Higher rank numerical ranges of normal matrices

The higher rank numerical range is closely connected to the construction of quantum error correction code for a noisy quantum channel. It is known that if a normal matrix $A \in M_n$ has eigenvalues $a_1, \..., a_n$, then its higher rank numerical range $Λ_k(A)$ is the intersection of convex polygons with vertices $a_{j_1}, \..., a_{j_{n-k+1}}$, where $1 \le j_1 < \... < j_{n-k+1} \le n$. In this paper, it is shown that the higher rank numerical range of a normal matrix with $m$ distinct eigenvalues can be written as the intersection of no more than $\max\{m,4\}$ closed half planes. In addition, given a convex polygon ${\mathcal P}$ a construction is given for a normal matrix $A \in M_n$ with minimum $n$ such that $Λ_k(A) = {\mathcal P}$. In particular, if ${\mathcal P}$ has $p$ vertices, with $p \ge 3$, there is a normal matrix $A \in M_n$ with $n \le \max\left\{p+k-1, 2k+2 \right\}$ such that $Λ_k(A) = {\mathcal P}$.

preprint2010arXiv

Interpolation problems by completely positive maps

Given commuting families of Hermitian matrices {A1, ..., Ak} and {B1, ...., Bk}, conditions for the existence of a completely positive map L, such that L(Aj) = Bj for j = 1, ...,k, are studied. Additional properties such as unital or / and trace preserving on the map ? are also considered. Connections of the study to dilation theory, matrix inequalities, unitary orbits, and quantum information science are mentioned.

preprint2008arXiv

Canonical forms, higher rank numerical range, convexity, totally isotropic subspace, matrix equations

Results on matrix canonical forms are used to give a complete description of the higher rank numerical range of matrices arising from the study of quantum error correction. It is shown that the set can be obtained as the intersection of closed half planes (of complex numbers). As a result, it is always a convex set in $\mathcal C$. Moreover, the higher rank numerical range of a normal matrix is a convex polygon determined by the eigenvalues. These two consequences confirm the conjectures of Choi et al. on the subject. In addition, the results are used to derive a formula for the optimal upper bound for the dimension of a totally isotropic subspace of a square matrix, and verify the solvability of certain matrix equations.

preprint2008arXiv

Higher rank numerical ranges and low rank perturbations of quantum channels

For a positive integer $k$, the rank-$k$ numerical range $Λ_k(A)$ of an operator $A$ acting on a Hilbert space $\cH$ of dimension at least $k$ is the set of scalars $λ$ such that $PAP = λP$ for some rank $k$ orthogonal projection $P$. In this paper, a close connection between low rank perturbation of an operator $A$ and $Λ_k(A)$ is established. In particular, for $1 \le r < k$ it is shown that $Λ_k(A) \subseteq Λ_{k-r}(A+F)$ for any operator $F$ with $\rank (F) \le r$. In quantum computing, this result implies that a quantum channel with a $k$-dimensional error correcting code under a perturbation of rank $\le r$ will still have a $(k-r)$-dimensional error correcting code. Moreover, it is shown that if $A$ is normal or if the dimension of $A$ is finite, then $Λ_k(A)$ can be obtained as the intersection of $Λ_{k-r}(A+F)$ for a collection of rank $r$ operators $F$. Examples are given to show that the result fails if $A$ is a general operator. The closure and the interior of the convex set $Λ_k(A)$ are completely determined. Analogous results are obtained for $Λ_\infty(A)$ defined as the set of scalars $λ$ such that $PAP = λP$ for an infinite rank orthogonal projection $P$. It is shown that $Λ_\infty(A)$ is the intersection of all $Λ_k(A)$ for $k = 1, 2, >...$. If $A - μI$ is not compact for any $μ\in \IC$, then the closure and the interior of $Λ_\infty(A)$ coincide with those of the essential numerical range of $A$. The situation for the special case when $A-μI$ is compact for some $μ\in \IC$ is also studied.

preprint2007arXiv

Condition for the higher rank numerical range to be non-empty

It is shown that the rank-$k$ numerical range of every $n$-by-$n$ complex matrix is non-empty if $n \ge 3k - 2$. The proof is based on a recent characterization of the rank-$k$ numerical range by Li and Sze, the Helly's theorem on compact convex sets, and some eigenvalue inequalities. In particular, the result implies that $Λ_2(A)$ is non-empty if $n \ge 4$. This confirms a conjecture of Choi et al. If $3k-2>n>0$, an $n$-by-$n$ complex matrix is given for which the rank-$k$ numerical range is empty. Extension of the result to bounded linear operators acting on an infinite dimensional Hilbert space is also discussed.