Source author record

Chi-Kwong Li

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

56works
15topics
4close collaborators

Published work

56 publications

preprint2026arXiv

Cliques and independent subgroups of the Birkhoff polytope graph

The Birkhoff polytope ΩnΩ_n is the polytope of doubly stochastic matrices of order nn. The Birkhoff polytope graph G(Ωn)G(Ω_n) is the skeleton of ΩnΩ_n; it is the Cayley graph whose vertex set consists of the elements of the symmetric group Sym(n){\rm Sym}(n) of degree nn, 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 Sym(n){\rm Sym}(n) whose elements are pairwise nonadjacent in the graph). We obtain maximal subgroups of G(Ωn)G(Ω_n) and establish both a lower bound and an upper bound for its clique number. Especially, we prove that if KK is a subset of Sym(n){\rm Sym}(n) consisting of 3-cycle permutations such that δ1−1δ2δ_1^{-1}δ_2 is a single cycle for all δ1,δ2∈Kδ_1,δ_2\in K, then the maximum size of KK is ⌊(n−1)2/4⌋\lfloor (n-1)^2/4\rfloor, which can be viewed as an Erdős-Ko-Rado-type theorem for Sym(n){\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×nn\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 H{\mathcal H} be a complex Hilbert space and let B(H){\mathcal B}({\mathcal H}) be the algebra of all bounded linear operators on H{\mathcal H}. For a positive integer kk less than the dimension of H{\mathcal H} and A=(A1,…,Am)∈B(H)m{\mathbf A} = (A_1, \dots, A_m)\in {\mathcal B}({\mathcal H})^m, the joint kk-numerical range Wk(A)W_k({\mathbf A}) is the set of (α1,…,αm)∈Cm(α_1, \dots, α_m) \in{\mathbb C}^m such that αi=∑j=1k⟨Aixj,xj⟩α_i = \sum_{j = 1}^k \langle A_ix_j, x_j\rangle for an orthonormal set {x1,…,xk}\{x_1, \ldots, x_k\} in H{\mathcal H}. Relations between the geometric properties of Wk(A)W_k({\mathbf A}) and the algebraic and analytic properties of A1,…,AmA_1, \dots, A_m are studied. It is shown that there is k∈Nk\in {\mathbb N} such that Wk(A)W_k({\mathbf A}) is a polyhedral set, i.e., the convex hull of a finite set, if and only if A1,…,AkA_1, \dots, A_k have a common reducing subspace V{\mathbf V} of finite dimension such that the compression of A1,…,AmA_1, \dots, A_m on the subspace V{\mathbf V} are diagonal operators D1,…,DmD_1, \dots, D_m and Wk(A)=Wk(D1,…,Dm)W_k({\mathbf A}) = W_k(D_1, \dots, D_m). Characterization is also given to A{\bf A} such that the closure of Wk(A)W_k({\mathbf A}) is polyhedral. The conditions are related to the joint essential numerical range of A{\mathbf A}. These results are used to study A{\bf A} such that (a) {A1,…,Am}\{A_1, \dots, A_m\} is a commuting family of normal operators, or (b) Wk(A1,…,Am)W_k(A_1, \dots, A_m) is polyhedral for every positive integer kk. 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 A1,…,AmA_1, \dots, A_m satisfying (a) and (b), respectively. Results are also obtained for general non-compact operators.

preprint2022arXiv

The joint kk-numerical range of operators

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

preprint2022arXiv

Unitarily invariant Norms on Operators

Let ff be a symmetric norm on Rn{\mathbb R}^n and let B(H){\mathcal B}({\mathcal H}) be the set of all bounded linear operators on a Hilbert space H{\mathcal H} of dimension at least nn. Define a norm on B(H){\mathcal B}({\mathcal H}) by ∥A∥f=f(s1(A),…,sn(A))\|A\|_f = f(s_1(A), \dots, s_n(A)), where sk(A)=inf⁡{∥A−X∥:X∈B(H) has rank less than k}s_k(A) = \inf\{\|A-X\|: X\in {\mathcal B}({\mathcal H}) \hbox{ has rank less than } k\} is the kkth singular value of AA. Basic properties of the norm ∥⋅∥f\|\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 LL satisfying ∥L(A)−L(B)∥f=∥A−B∥f\|L(A)-L(B)\|_f=\|A - B\|_f for any A,B∈B(H)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α⊗Nx=0A_α^{\otimes N} {\bf x} = 0 arising from the distinguishability of two quantum operations by NN uses in parallel, where the coefficient matrix Aα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 RNR_N depending on NN. We affirm the necessity part of the conjecture and establish the sufficiency of the conjecture for N≤10N\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 nn-qubit system with error operators that assume the form σx⊗nσ_x^{\otimes n}, σy⊗nσ_y^{\otimes n}, σz⊗nσ_z^{\otimes n}. Previous schemes are improved to facilitate implementation. In particular, when nn is odd and equals 2k+12k+1, we describe a quantum error correction scheme using one arbitrary qubit σσ to protect the data state ρρ in a 2k2k-qubit system. The encoding operation σ⊗ρ↦Φ(σ⊗ρ)σ\otimes ρ\mapsto Φ(σ\otimes ρ) only requires 3k3k 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Φ^{-1} to produce σ~⊗ρ\tilde σ\otimes ρ so that a partial trace operation can recover ρρ. When nn is even and equals 2k+22k+2, we describe a hybrid quantum error correction scheme using any one of the two classical bits σ∈{∣ij⟩⟨ij∣:i,j∈{0,1}}σ\in \{|ij\rangle \langle ij|: i, j \in \{0,1\}\} to protect a 2k2k-qubit state ρρ and 2 classical bits. The encoding operation σ⊗ρ↦Φ(σ⊗ρ)σ\otimes ρ\mapsto Φ(σ\otimes ρ) can be done by 3k+23k+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Φ^{-1} to produce σ⊗ρσ\otimes ρ so that a perfect protection of the two classical bits σσ and the 2k2k-qubit state is achieved. If one uses an arbitrary 22-qubit state σσ, the same scheme will protect 2k2k-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×nn\times n matrices and the generalized joint numerical ranges is studied. For instance, it is shown that F{\cal F} is a family of mutually commuting normal matrices if and only if the joint numerical range Wk(A1,…,Am)W_k(A_1, \dots, A_m) is a polyhedral set for some kk satisfying ∣n/2−k∣≤1|n/2-k|\le 1, where {A1,…,Am}\{A_1, \dots, A_m\} is a basis for the linear span of the family; equivalently, Wk(X,Y)W_k(X,Y) is polyhedral for any two X,Y∈FX, Y \in {\cal F}. More generally, characterization is given for the cc-numerical range Wc(A1,…,Am)W_c(A_1, \dots, A_m) to be polyhedral for any n×nn\times n matrices A1,…,AmA_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 nn, are those linear maps that can be expressed as a convex combination of conjugations by n×nn\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~NN and the Choi rank~rr 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≤r2−r+1N\leq r^2-r+1 is satisfied for every mixed-unitary channel (as is the equality N=2N=2 when r=2r=2), and we exhibit the first known examples of mixed-unitary channels for which N>rN>r. Specifically, we prove that there exist mixed-unitary channels having Choi rank d+1d+1 and mixed-unitary rank 2d2d for infinitely many positive integers dd, including every prime power dd. 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∈{1}∪[2,∞)r \in \{1\} \cup [2, \infty), positive semi-definite matrices Ai, Bi, Ci∈MniA_i,\ B_i,\ C_i\in M_{n_i}, i=1,2i=1,2, and generalized matrix functions dχ,dξ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 mm positive semi-definite matrices for m≥3m \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 AA and BB such that W(A)W(A) and W(B)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 ABAB, and conditions on AA and BB for the set equality holding are obtained. The results cover the case when AA and BB 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ρ_1, ρ_2 be quantum states and (ρ1,ρ2)↦D(ρ1,ρ2)(ρ_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))D(ρ_1, Φ(ρ_2)) for Φ∈SΦ\in \mathcal{S} for different class of functions D(⋅,⋅)D(\cdot, \cdot), where S\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 ℓ1\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 E(d,n)={U⊗n:U∈SU(d)}.{\mathcal E}(d,n) = \{ U^{\otimes n}: U \in {\mathrm{SU}}(d)\}. This is done by analyzing the decomposition of the algebra A(d,n){\mathcal A}(d,n) generated by E(d,n){\mathcal E}(d,n). We summarize the results for the channels on qubits (d=2d=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 dd. 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=3d = 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 K1,K2K_1, K_2 be two compact convex sets in C\mathit{C}. Their Minkowski product is the set K1K2={ab:a∈K1,b∈K2}K_1K_2 = \{ab: a \in K_1, b\in K_2\}. We show that the set K1K2K_1K_2 is star-shaped if K1K_1 is a line segment or a circular disk. Examples for K1K_1 and K2K_2 are given so that K1K_1 and K2K_2 are triangles (including interior) and K1K2K_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 kk positive semi-definite matrices but not fewer, for k=1,2,3,4,5k = 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 TT be a quadratic operator on a complex Hilbert space HH. We show that TT can be written as a product of two positive contractions if and only if TT is of the form aI⊕bI⊕(aIP0bI)onH1⊕H2⊕(H3⊕H3)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∈[0,1]a, b\in [0,1] and strictly positive operator PP with ∥P∥≤∣a−b∣(1−a)(1−b).\|P\| \le |\sqrt{a} - \sqrt{b}|\sqrt{(1-a)(1-b)}. Also, we give a necessary condition for a bounded linear operator TT with operator matrix (T1T30T2)\begin{pmatrix} T_1 & T_3\\ 0 & T_2\cr\end{pmatrix} on H⊕KH\oplus K that can be written as a product of two positive contractions.

preprint2014arXiv

Perturbing eigenvalues of non-negative matrices

Let AA be an irreducible (entrywise) nonnegative n×nn\times n matrix with eigenvalues ρ,b+ic,b−ic,λ4,⋯ ,λn,ρ, b+ic,b-ic, λ_4,\cdots,λ_n, where ρρ is the Perron eigenvalue. It is shown that for any t∈[0,∞)t \in [0, \infty) there is a nonnegative matrix with eigenvalues ρ+t~,λ2+t,λ3+t,λ4⋯ ,λn,ρ+ \tilde t,λ_2+t,λ_3+t, λ_4 \cdots,λ_n, whenever t~≥γnt\tilde t \ge γ_n t with γ3=1,γ4=2,γ5=5γ_3=1, γ_4 = 2, γ_5=\sqrt 5 and γn=2.25γ_n = 2.25 for n≥6n \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 nn-sided convex polygon is bounded by γnγ_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×33\times 3 matrix of the form (∥Ri∗Rj∥tr)1≤i,j≤3(\|R_i^*R_j\|_{\rm tr})_{1 \le i, j \le 3} for any rectangular complex (or real) matrices R1,R2,R3R_1, R_2, R_3 so that the multiplication Ri∗RjR_i^*R_j is compatible for all i,ji, j, where ∥⋅∥tr\|\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 MnM_n the set of n×nn\times n complex matrices. Let f:Mn→[0,∞)f: M_n \rightarrow [0,\infty) be a continuous map such that f(μUAU∗)=f(A)f(μUAU^*)= f(A) for any complex unit μμ, A∈MnA \in M_n and unitary U∈MnU \in M_n, f(X)=0f(X)=0 if and only if X=0X=0 and the induced map t↦f(tX)t \mapsto f(tX) is monotonic increasing on [0,∞)[0,\infty) for any rank 1 nilpotent X∈MnX \in M_n. Characterizations are given for surjective maps ϕϕ on MnM_n satisfying f(AB−BA)=f(ϕ(A)ϕ(B)−ϕ(B)ϕ(A))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,…,ρk}\{ρ_1, \dots, ρ_k\} to another such set {ρ^1,…,ρ^k}\{\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×nn\times n density matrix σ∈Mnσ\in M_n, let S(σ){\cal S}(σ) be the compact convex set of quantum states ρ=(ρij)∈Mm⋅nρ= (ρ_{ij}) \in M_{m\cdot n} with the first partial trace equal to σσ, i.e., tr1(ρ)=ρ11+⋯+ρmm=σ{\rm tr}_1(ρ) = ρ_{11} + \cdots + ρ_{mm} = σ. It is known that if m≥nm \ge n then there is a rank one matrix ρ∈S(σ)ρ\in {\cal S}(σ) satisfying tr1(ρ)=σ{\rm tr}_1(ρ) = σ. If m<nm < n, there may not be rank one matrix in S(σ){\cal S}(σ). In this paper, we determine the ranks of the elements and ranks of the extreme points of the set S{\cal S} We also determine ρ∗∈S(σ)ρ^* \in {\cal S}(σ) with rank bounded by kk such that ∥tr1ρ∗)−σ∥\|{\rm tr}_1ρ^*) - σ\| is minimum for a given unitary similarity invariant norm ∥⋅∥\|\cdot\|. Furthermore, the relation between the eigenvalues of σσ and those of ρ∈S(σ)ρ\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(dd) errors. We find that under this recursive scheme, the asymptotic encoding rate is 1/d1/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 AA is a multiple of a positive semi-definite operator if and only if σ(AB)⊆W(A)W(B)‾,for all (rank one) operators B. σ(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 AA so that equivalence conditions hold.

preprint2013arXiv

Decomposition of quantum gates

A recurrence scheme is presented to decompose an nn-qubit unitary gate to the product of no more than N(N−1)/2N(N-1)/2 single qubit gates with small number of controls, where N=2nN = 2^n. Detailed description of the recurrence steps and formulas for the number of kk-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 dd-by-dd 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−1d-1 classes, where each class is isomorphic to the group of 2×22\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 nn-qubit registers can be expressed as no more than 2n−1(2n−1)2^{n-1}(2^n-1) fully controlled single-qubit gates chosen from 2n−12^n-1 classes, where the quantum gates in each class share the same n−1n-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 A11∈MmA_{11} \in M_m with m≤n/2m \le n/2, be such that the numerical range of AA lies in the set $\{e^{iφ} z \in \IC: |\Im z| \le (\Re z) \tan α\}$, for some φ∈[0,2π)φ\in [0, 2π) and α∈[0,π/2)α\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 Im−A11−1A12A22−1A21I_m - A_{11}^{-1}A_{12} A_{22}^{-1}A_{21} in case A11A_{11} and A22A_{22} are invertible. From this result, one can show ∣det⁡(A)∣≤sec⁡2m(α)∣det⁡(A11)det⁡(A22)∣|\det(A)| \le \sec^{2m}(α) |\det(A_{11})\det(A_{22})|. In particular, if AA is a accretive-dissipative matrix, then ∣det⁡(A)∣≤2m∣det⁡(A11)det⁡(A22)∣|\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 nn, let MnM_n be the set of n×nn\times n complex matrices. Suppose ∥⋅∥\|\cdot\| is the Ky Fan kk-norm with 1≤k≤mn1 \le k \le mn or the Schatten pp-norm with 1≤p≤∞1 \le p \le \infty (p≠2p\ne 2) on MmnM_{mn}, where m,n≥2m,n\ge 2 are positive integers. It is shown that a linear map ϕ:Mmn→Mmnϕ: M_{mn} \rightarrow M_{mn} satisfying ∥A⊗B∥=∥ϕ(A⊗B)∥for allA∈MmandB∈Mn\|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∈MmnU, V \in M_{mn} such that ϕϕ has the form A⊗B↦U(φ1(A)⊗φ2(B))VA\otimes B \mapsto U(φ_1(A) \otimes φ_2(B))V, where φs(X)φ_s(X) is either the identity map X↦XX \mapsto X or the transposition map X↦XtX \mapsto X^t. The results are extended to tensor space Mn1⊗⋯⊗MnmM_{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 S(H){\mathcal S}(H) the convex set of all states (i.e., trace-1 positive operators) on Hilbert space HH with dimH≤∞H\leq \infty, and [ρ1,ρ2][ρ_1, ρ_2] the line segment joining two elements ρ1,ρ2ρ_1, ρ_2 in S(H){\mathcal S}(H). It is shown that a bijective map ϕ:S(H)→S(H)ϕ:{\mathcal S}(H) \rightarrow {\mathcal S}(H) satisfies ϕ([ρ1,ρ2])⊆[ϕ(ρ1),ϕ(ρ2)]ϕ([ρ_1, ρ_2]) \subseteq [ϕ(ρ_1),ϕ(ρ_2)] for any ρ1,ρ2∈Sρ_1, ρ_2 \in {\mathcal S} if and only if ϕϕ has one of the following forms ρ↦MρM∗tr(MρM∗)orρ↦MρTM∗tr(MρTM∗),ρ\mapsto \frac{MρM^*}{{\rm tr}(MρM^*)}\quad \hbox{or} \quad ρ\mapsto \frac{Mρ^T M^*}{{\rm tr}(Mρ^T M^*)}, where MM is an invertible bounded linear operator and ρTρ^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≥2m,n\ge 2 be positive integers, MmM_m the set of m×mm\times m complex matrices and MnM_n the set of n×nn\times n complex matrices. Regard MmnM_{mn} as the tensor space Mm⊗MnM_m\otimes M_n. Suppose ∣⋅∣|\cdot| is the Ky Fan kk-norm with 1≤k≤mn1 \le k \le mn, or the Schatten pp-norm with 1≤p≤∞1 \le p \le \infty (p≠2p\ne 2) on MmnM_{mn}. It is shown that a linear map ϕ:Mmn→Mmnϕ: M_{mn} \rightarrow M_{mn} satisfying ∣A⊗B∣=∣ϕ(A⊗B)∣|A\otimes B| = |ϕ(A\otimes B)| for all A∈MmA \in M_m and B∈MnB \in M_n if and only if there are unitary U,V∈MmnU, V \in M_{mn} such that ϕϕ has the form A⊗B↦U(φ1(A)⊗φ2(B))VA\otimes B \mapsto U(φ_1(A) \otimes φ_2(B))V, where φi(X)φ_i(X) is either the identity map X↦XX \mapsto X or the transposition map X↦XtX \mapsto X^t. The results are extended to tensor space Mn1⊗...⊗MnmM_{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⊗nσ_x^{\otimes n}, σy⊗nσ_y^{\otimes n} and σz⊗nσ_z^{\otimes n}, where n>2n>2 is the number of qubits encoding the codeword. It is proved that (i) nn qubits codeword encodes (n−1)(n-1) data qubits when nn is odd and (ii) nn qubits codeword implements a noiseless subsystem encoding (n−2)(n-2) data qubits when nn 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 nn of physical qubits. A single logical qubit is encoded with n=3n=3 and n=4n=4, while two logical qubits are encoded with n=5n=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 mm logical qubits increases linearly in mm.

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 A1,A2A_1, A_2 on complex Banach spaces X1,X2X_1, X_2, respectively. This includes the usual Jordan product A1∘A2=A1A2+A2A1A_1 \circ A_2 = A_1 A_2 + A_2 A_1, and the triple {A1,A2,A3}=A1A2A3+A3A2A1\{A_1,A_2,A_3\} = A_1 A_2 A_3 + A_3 A_2 A_1. Let a map Φ:A1→A2Φ: A_1 \to A_2 prserving the spectra of the products σ(Φ(A1)∘...∘Φ(Ak))=σ(A1∘...∘Ak) σ(Φ(A_1) \circ ... \circ Φ(A_k)) = σ(A_1\circ ... \circ A_k) whenever any one of A1,...,AkA_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 mmth 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∈MnA \in M_n has eigenvalues a1,.˙.,ana_1, \..., a_n, then its higher rank numerical range Λk(A)Λ_k(A) is the intersection of convex polygons with vertices aj1,.˙.,ajn−k+1a_{j_1}, \..., a_{j_{n-k+1}}, where 1≤j1<.˙.<jn−k+1≤n1 \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 mm distinct eigenvalues can be written as the intersection of no more than max⁡{m,4}\max\{m,4\} closed half planes. In addition, given a convex polygon P{\mathcal P} a construction is given for a normal matrix A∈MnA \in M_n with minimum nn such that Λk(A)=PΛ_k(A) = {\mathcal P}. In particular, if P{\mathcal P} has pp vertices, with p≥3p \ge 3, there is a normal matrix A∈MnA \in M_n with n≤max⁡{p+k−1,2k+2}n \le \max\left\{p+k-1, 2k+2 \right\} such that Λk(A)=PΛ_k(A) = {\mathcal P}.

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 C\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 kk, the rank-kk numerical range Λk(A)Λ_k(A) of an operator AA acting on a Hilbert space $\cH$ of dimension at least kk is the set of scalars λλ such that PAP=λPPAP = λP for some rank kk orthogonal projection PP. In this paper, a close connection between low rank perturbation of an operator AA and Λk(A)Λ_k(A) is established. In particular, for 1≤r<k1 \le r < k it is shown that Λk(A)⊆Λk−r(A+F)Λ_k(A) \subseteq Λ_{k-r}(A+F) for any operator FF with $\rank (F) \le r$. In quantum computing, this result implies that a quantum channel with a kk-dimensional error correcting code under a perturbation of rank ≤r\le r will still have a (k−r)(k-r)-dimensional error correcting code. Moreover, it is shown that if AA is normal or if the dimension of AA is finite, then Λk(A)Λ_k(A) can be obtained as the intersection of Λk−r(A+F)Λ_{k-r}(A+F) for a collection of rank rr operators FF. Examples are given to show that the result fails if AA is a general operator. The closure and the interior of the convex set Λk(A)Λ_k(A) are completely determined. Analogous results are obtained for Λ∞(A)Λ_\infty(A) defined as the set of scalars λλ such that PAP=λPPAP = λP for an infinite rank orthogonal projection PP. It is shown that Λ∞(A)Λ_\infty(A) is the intersection of all Λk(A)Λ_k(A) for k=1,2,>...k = 1, 2, >.... If A−μIA - μI is not compact for any $μ\in \IC$, then the closure and the interior of Λ∞(A)Λ_\infty(A) coincide with those of the essential numerical range of AA. The situation for the special case when A−μIA-μ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-kk numerical range of every nn-by-nn complex matrix is non-empty if n≥3k−2n \ge 3k - 2. The proof is based on a recent characterization of the rank-kk 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)Λ_2(A) is non-empty if n≥4n \ge 4. This confirms a conjecture of Choi et al. If 3k−2>n>03k-2>n>0, an nn-by-nn complex matrix is given for which the rank-kk numerical range is empty. Extension of the result to bounded linear operators acting on an infinite dimensional Hilbert space is also discussed.

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.