Source author record

Yiu-Tung Poon

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

26works
7topics
4close collaborators

Published work

26 publications

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

NISQ: Error Correction, Mitigation, and Noise Simulation

Error-correcting codes were invented to correct errors on noisy communication channels. Quantum error correction (QEC), however, may have a wider range of uses, including information transmission, quantum simulation/computation, and fault-tolerance. These invite us to rethink QEC, in particular, about the role that quantum physics plays in terms of encoding and decoding. The fact that many quantum algorithms, especially near-term hybrid quantum-classical algorithms, only use limited types of local measurements on quantum states, leads to various new techniques called Quantum Error Mitigation (QEM). This work examines the task of QEM from several perspectives. Using some intuitions built upon classical and quantum communication scenarios, we clarify some fundamental distinctions between QEC and QEM. We then discuss the implications of noise invertibility for QEM, and give an explicit construction called Drazin-inverse for non-invertible noise, which is trace preserving while the commonly-used Moore-Penrose pseudoinverse may not be. Finally, we study the consequences of having an imperfect knowledge about the noise, and derive conditions when noise can be reduced using QEM.

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}.

preprint2021arXiv

Observing geometry of quantum states in a three-level system

In quantum mechanics, geometry has been demonstrated as a useful tool for inferring non-classical behaviors and exotic properties of quantum systems. One standard approach to illustrate the geometry of quantum systems is to project the quantum state space to the Euclidean space via measurements of observables on the system. Despite the great success of this method in studying two-level quantum systems (qubits) with the celebrated Bloch sphere representation, there is always the difficulty to reveal the geometry of multi-dimensional quantum systems. Here we report the first experiment measuring the geometry of such projections beyond the qubit. Specifically, we observe the joint numerical ranges (JNRs) of a triple of observables in a three-level photonic system, providing complete classification of the JNRs. We further show that the geometry of different classes reveal ground-state degeneracies of a Hamiltonian as a linear combination of the observables, which is related to quantum phases in the thermodynamic limit. Our results offer a versatile geometric approach for exploring the properties of higher-dimensional quantum systems.

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

Experimental self-characterization of quantum measurements

The accurate and reliable description of measurement devices is a central problem in both observing uniquely non-classical behaviors and realizing quantum technologies from powerful computing to precision metrology. To date quantum tomography is the prevalent tool to characterize quantum detectors. However, such a characterization relies on accurately characterized probe states, rendering reliability of the characterization lost in circular argument. Here we report a self-characterization method of quantum measurements based on reconstructing the response range, the entirety of attainable measurement outcomes, eliminating the reliance on known states. We characterize two representative measurements implemented with photonic setups and obtain fidelities above 99.99% with the conventional tomographic reconstructions. This initiates range-based techniques in characterizing quantum systems and foreshadows novel device-independent protocols of quantum information applications.

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.

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.

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.

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

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

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.

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

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

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.