Processing math: 100%

数学物理学报, 2022, 42(1): 1-8 doi:

论文

具有已知约化态的无限维量子系统态的秩

杨舒媛,1,3, 贺衎,2,3

1 太原理工大学信息与计算机学院 太原 030024

2 中国科学院量子信息重点实验室 合肥 230026

3 太原理工大学数学学院 太原 030024

Ranks of Quantum States with Prescribed Reduced States in an Infinite Dimensional Quantum System

Yang Shuyuan,1,3, He Kan,2,3

1 College of Information and Computer, Taiyuan University of Technology, Taiyuan 030024

2 Key Laboratory of Quantum Information, University of Science and Technology of China, Chinese Academy of Sciences, Hefei 230026

3 College of Mathematics, Taiyuan University of Technology, Taiyuan 030024

通讯作者: 贺衎, E-mail: hekanquantum@163.com

收稿日期: 2021-01-27  

基金资助: 国家自然科学基金.  11771011
国家自然科学基金.  12071336
中国科学院量子信息重点实验室开放课题

Received: 2021-01-27  

Fund supported: the NSFC.  11771011
the NSFC.  12071336
the Key Laboratory of Quantum Information, University of Science and Technology of China

作者简介 About authors

杨舒媛,E-mail:yangshuyuan2000@163.com , E-mail:yangshuyuan2000@163.com

Abstract

Suppose that H and K are two infinite dimensional quantum systems (i.e. infinite dimensional complex Hilbert space). Let ρ be a quantum state on HK with two reduced states trH(ρ) and trK(ρ). Then all possible ranks of ρ are determined in this paper.

Keywords: Quantum state ; Reduced state ; Partial trace ; Rank

PDF (869KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

杨舒媛, 贺衎. 具有已知约化态的无限维量子系统态的秩. 数学物理学报[J], 2022, 42(1): 1-8 doi:

Yang Shuyuan, He Kan. Ranks of Quantum States with Prescribed Reduced States in an Infinite Dimensional Quantum System. Acta Mathematica Scientia[J], 2022, 42(1): 1-8 doi:

1 介绍

在量子力学的数学框架中, 一个封闭的量子系统用复Hilbert空间H表示, 量子态为H上迹为1的正算子[1]. 令HK分别为两个无限维复Hilbert空间. 用B(H)表示H上的所有有界线性算子构成的代数, S(H)表示复Hilbert空间H上所有量子态. Mm,n表示m×n复矩阵的集合, 其中Mn=Mn,n. 量子态ρHKS(HK)的偏迹映射是一个满足下列运算性质的线性变换

trH(AB)=B,trK(AB)=A.

我们称trH(ρHK)trK(ρHK)是态ρHK的两个约化态[1].

在量子信息理论中, 一个有趣的问题是: 在已知态ρHK的两个约化态的情况下, 如何确定ρHK的秩. 对该问题的研究既受到了量子物理学者们的关注[2-4], 也被许多数学工作者探讨[5-7]. 最近, Li等在文献[3] 和[4]中, 定义了以下集合: 设HK为有限维的两个复Hilbert空间, 对于任意给定的AS(H)BS(K),集合S(A,B)={ρS(HK):trH(ρ)=B,trK(ρ)=A},并获得了计算该集合中任意元ρ秩的所有可能值的方法. 本文研究上述问题在无限维系统中的解决.

在阐述主要结果前, 我们介绍一些术语和标记. 在本文中, S(HK)中的两体态往往需要进行矩阵表示. 令{en}dHn=1{um}dKm=1分别为空间HK中的标准基(dH,dK分别是H,K的维数), 则{enum}dH,dKn=1,m=1成为HK上的标准基, 且eiutj=Gij. 量子态ρS(HK)可用矩阵表示为

ρ=dH,dKi,j=1;k,l=1tij,klEijFkl=dHi,j=1Eij(dKk,l=1tij,klFkl)=dHi,j=1Eijρij.

因此

ρ=(ρij)1i,jdH.

ρ的约化态形式如下

trK(ρ)=(tr(ρij))1i,jdHS(H);

trH(ρ)=dKn=1ρnnS(K).

我们用XtX分别表示一个算子或者向量X的转置和共轭转置

H={AB(H):A=A}  K={BB(K):B=B}.

当空间维数为n时, 将H,K分别记为Hn,Kn.UV分别表示B(H)B(K)上酉算子的集合. Rn表示Rnn元组按递减顺序排列后的集合. 若σH,那么λ(σ)=(λ1(σ),,λn(σ))Rn表示按递减顺序排列的σ的特征值. 为了明确偏迹与算子结构间的关系, 我们首先需要考虑下面的两个线性映射[4]

[]:CmnMm,n,vec:Mm,nCmn.

这两个映射在有限维的情况下已经被广泛的讨论. 这里, 对于w=(w1,,wmn)tCmn,W=[w]m×n矩阵, 其中第j列的元素为(w(j1)n+1,,wjn)t,vec是把m×n矩阵转化为w=vec(W)Cmn的逆映射, 使得W=[w]. 值得注意的是

trH(ww)=Wt(Wt),trK(ww)=WW.

我们可以将上述有限维技巧推广到无限维情形. 假设dimHK=+, TC2(HK), ZC2(KK,HH)(C2(H,K)表示由B(H,K)中的Hilbert-Schmidt算子组成的Hilbert空间, 当H=K时, 简记为C2(H)). {|m}{|μ}分别为HK的标准基. 则TZ有矩阵表示

T=(tmμ,nν),tmμ,nν=m|μ|T|n|ν

Z=(zmn,μν),zmn,μν=m|n|Z|μ|ν.

如果

zmn,μν=tmμ,nν,

那么ZT关于基|m|μ的重排算子, 记为TR=Z. 最后, 我们观察下面的酉相似不变性是成立的, 对于任意酉算子UV, 有

S(UAU,VBV)=(UV)S(A,B)(UV)={(UV)ρ(UV):ρS(A,B)}.

2 S(A,B)中算子的秩

在这一节中, 考虑H,K为无限维Hilbert空间的情况. 对于任意给定的AS(H)BS(K),S(A,B)={ρS(HK):trH(ρ)=B,trK(ρ)=A},我们确定了该集合S(A,B)中量子态的秩的范围.

定理2.1   设(A,B)S(H)×S(K)rank(A)=r1r2=rank(B)<+.S(A,B)中存在一个元素达到最小秩rrr1. 那么, S(A,B)中存在秩为k的量子态的充要条件是rkr1r2.

  由集合S(A,B)的酉相似不变性可知

A=diag(λ1,,λr1,0,,0,0,),B=diag(ˆλ1,,ˆλr2,0,,0,0,).

如果ρ=kj=1zjzjS(A,B),那么trH(ρ)=BtrK(ρ)=A.

由于量子态的迹为1, 那么约化态的秩已知的量子态可以利用这个特征来表示, 具体如下

1ir11jr2时, 令vij=λiˆλjeiuj. 对于l=1,,r1,

zl={1ir2vi+l1,i,1lr1r2+1.1ir1+1lvi+l1,i+r1+1l<ir2vi+l1r1,i,r1r2+1<lr1.

那么, 利用偏迹的计算可以得到

trH(zlzl)={1ir2λi+l1ˆλiGii,1lr1r2+1.1ir1+1lλi+l1ˆλiGii+r1+1l<ir2λi+l1r1ˆλiGii,r1r2+1<lr1.
(2.1)

同理, 可以得到

trK(zlzl)={1ir2λi+l1ˆλiGi+l1i+l1,1lr1r2+1.1ir1+1lλi+l1ˆλiGi+l1i+l1+r1+1l<ir2λi+l1r1ˆλiGi+l1r1i+l1r1,r1r2+1<lr1.
(2.2)

首先, 证明S(A,B)中包含秩为r1,,r1r2的量子态. 令

ρ=r1l=1zlzl,
(2.3)

那么ρ的秩为r1. 由(2.1) 式和(2.2) 式中偏迹的具体表示形式以及ρ的计算方式(2.3) 式可见ρS(A,B). 在(2.3) 式中, 如果我们用

v11v11++vppvpp+(j>pvjj)(j>pvjj),p=1,,r21

代替z1z1,那么S(A,B)中量子态的秩为p+r1,其中p=1,,r21. 类似地, 我们可以用秩为p+1的算子来代替每一个zjzj,其中p=1,,r21. 通过这样的替换之后生成的量子态还是属于集合S(A,B). 因此, S(A,B)中包含秩为r1,,r1r2的量子态.

其次, 利用归纳法证明S(A,B)中包含秩为r,r+1,r1的量子态. 由于S(A,B)中存在量子态ρ,其秩为rr<r1. 那么存在

ρ=rj=1zjzjS(A,B).
(2.4)

利用Schmidt分解可知, 对于任意j,zj=tjl=1sjlxjlyjl,其中tj[zj]的秩. 令wjl=sjlxjlyjl. 类似于前面的证明方法, 我们可以用

lpwjlwjl+(l>pwjl)(l>pwjl)

代替zjzj. 显然, 替换之后的量子态仍然属于S(A,B). 因此, 集合S(A,B)中包含秩为r,r+1,r1的量子态.

结合上面两部分的证明(实际上具体构造出了集合中满足条件的量子态)可以说明S(A,B)中存在量子态的秩为k的充要条件是rkr1r2.

接下来, 我们通过例子来帮助理解上面的定理.

例2.1  若A=diag(a1,a2,a3,,0,0,)S(H)B=diag(b1,b2,,0,0,)S(K),其中AB的特征值由递减顺序排列, rank(A)=r1=3rank(B)=r2=2. 由定理2.1的结论表明:如果ρS(A,B),那么集合S(A,B)中量子态的秩最大为6. 具体构造方式如下

z1=v11+v22=((a1b1)1/2,0,0,,0,,0,(a2b2)1/2,0,,0,0,,0,)t,z2=v21+v32=(0,0,0,,0,,(a2b1)1/2,0,0,,0,0,(a3b2)1/2,0,,0,,0,,0,)t,z3=v31+v12=(0,(a1b2)1/2,0,0,,0,0,0,,0,,(a3b1)1/2,0,0,,0,,0,,0,)t.

那么, ρ1=z1z1+z2z2+z3z3S(A,B)的秩为3. 利用定理2.1证明过程中的替换方式, 我们可以得到ρ2=v11v11+v22v22+z2z2+z3z3S(A,B), 其秩为4. ρ3=v11v11+v22v22+v21v21+v32v32+z3z3S(A,B), 其秩为5. ρ4=v11v11+v22v22+v21v21+v32v32+v31v31+v12v12S(A,B), 其秩为6.

推论2.1  设(A,B)S(H)×S(K), 其中rank(A)<+,rank(B)<+, 那么S(A,B)中存在一个秩为1的量子态ρ的充分必要条件是AB有相同的非零特征值.

进一步可以说明, 存在秩为k的量子态ρS(A,B)当且仅当1kr2, 其中rAB特征值个数.

  若ρ是一个秩为1的量子态, 那么ρ可以表示为ρ=ννS(A,B), 这里ν有奇异值分解ν=i=1si(ν)figi,其中si(ν)ν的奇异值, figi是秩一投影. 而且约化态A,B的非零特征值的个数有限. 因此, ρ=νν当且仅当A=ri=1s2i(ν)fifiB=ri=1s2i(ν)gigi.

AB有相同的非零特征值, 即它们有相同的秩r,从而结合定理2.1即可得出结论:存在秩为k的量子态ρS(A,B)当且仅当1kr2.

同样地, 我们通过例子来帮助理解上面的定理.

例2.2   若A=diag(λ1,λ2,0,,0,)S(H)B=diag(λ1,λ2,0,,0,)S(K),其中λ1λ2>0. 推论2.1表明S(A,B)中量子态ρ的秩r的范围是1r4.z1=ν11+ν22=(λ1,0,0,0,,λ2,0,0,,0,)t, 那么ρ1=z1z1S(A,B),其秩为1.z2=ν21+ν12=(0,(λ1λ2)1/2,0,0,,0,(λ2λ1)1/2,0,,0,)t,那么ρ2=z1z1+z2z2S(A,B),其秩为2. 同理, 可以得到ρ3=(ν11+ν12)(ν11+ν12)+ν21ν21+ν31ν31S(A,B),其秩为3. ρ4=ν11ν11+ν12ν12+ν21ν21+ν31ν31S(A,B),其秩为4.

ρS(A,B),其秩小于等于r, 那么一定存在一个dHdK×r的算子D, 使得trH(DD)=BtrK(DD)=A.Wi=[vi], 其中viD的第i列的元素, 且1ir. 那么可以得到A=W1(W1)++Wr(Wr)B=W1t(W1t)++Wrt(Wrt).

为了确定集合S(A,B)中秩最小的量子态, 我们需要研究特征值与算子和之间的关系. 幸运的是已经有学者研究了这两者之间的关系[5-7], 其中Bercovici等[7] 证明了若一个正紧算子可以表示为一列正紧算子的和, 那么其特征值具有Horn不等式关系. 因此, 首先指出下面的引理, 这一引理为无限维的推广奠定了基础.

引理2.1   令α,β(1),,β(m)是一些极限为0的递减序列, 那么下面两者等价

(1) 存在正紧算子A,B(1),,B(k), 使得A=mk=1B(k), 其中Λ+(A)=α, Λ+(B(k))=β(k), k=1,2,,m.

(2) 对于任意r<N, q1,,qm0, 使得q=q1++qmr, 并且(I,J1,,Jm)TNr(m+1), 有Horn不等式

iIαimk=1jJkβ(k)j
(2.5)

和不等式

iIcqαimk=1jJ(k)cqkβ(k)j,
(2.6)

其中, TNr(m+1)是指[N]=1,2,,NN个正整数中任意取r个数, 每次取完后构成新的集合, 这样重复选取m+1次之后构成的m+1个元组I,J1,,Jm的集合, 且0rN.N是特征值的个数, qk表示正特征值的个数, Ic=NI表示I的补, Icq表示Icq个最小元素组成的集合. Λ+(A),Λ+(B(k))表示各自的特征值序列.

下面给出本文的主要定理, 引理2.1在该定理的证明过程中起着重要的作用.

定理2.2  若AS(H)BS(K). 则下列三者等价

(1) 存在ρS(A,B), 其秩小于等于r.

(2) 存在C1,,CrH˜C1,,˜CrK, 使得

(i)λ(˜Cj)=λ(CjO),(ii)A=rj=1Cj,(iii)B=rj=1˜Cj.

(3) 存在CH, 使得λ(Cj)=λ(C), λ(˜Cj)=λ(COdKdH), 其中Cj˜Cj满足(2), 1jr.

  (1)(2)S(A,B)中存在秩为r的量子态ρ, 那么存在W1,,Wr, 使得

A=W1(W1)++Wr(Wr),B=Wt1(Wt1)++Wtr(Wtr).

Cj=WjWj˜Cj=Wjt(Wjt), 显然Cj˜Cj满足条件(2).

(2)(1)利用(1)(2)证明过程中产生的Wi, 取

ρ=(vec(W1),vec(W2),,vec(Wr))(vec(W1),vec(W2),,vec(Wr)).

容易证明该ρS(A,B), 且其秩小于等于r.

(2)(3)假设λ(A)=(a1,a2,,am,0,0,0,), λ(Cj)=(c1j,c2j,,cmj,0,0,0,), j=1,,r.(I,J1,,Jm)TNr(m+1), 且π=(1,2,,r)1,2,,r的循环重排, 具体地, 有π(j)=j+1, 1j<rπ(r)=1. 那么(I,Jπk(1),Jπk(2),,Jπk(r))TNr(m+1). 因为A=C1++Cr, 由引理2.1, 我们有

iIairj=1iJjciπk(j),

其中0kr.

因此, 上面的式子具体化之后可以得到

iIaijJ1ci2+jJ2ci3++jJr1cir+jJrci1,k=1.

iIaijJ1ci3+jJ2ci4++jJr1ci1+jJrci2,k=2.

iIaijJ1ci4+jJ2ci5++jJr1ci2+jJrci3,k=3.

iIaijJ1ci1+jJ2ci2++jJr1cir1+jJrcir,k=r.

将这r个式子相加后可以得到

riIairj=1iJjci1+ci2++cir.

iIairj=1iJjci,ci=(ci1+ci2++cir)/r.

类似地, 可以说明逆向不等式也是成立的.

如果对于1jr, 取λ(ˆCj)=(c1,c2,,cm), 其中ˆC1,ˆC2,,ˆCrH. 我们可以得到

A=ˆC1+ˆC2++ˆCr.

类似地, 我们可以选择˜Cj, 使得λ(˜Cj)=(c1,c2,,cm,0,0,,0,0,). 此时, 则可以说明确实存在CH, 使得(2)成立, 取C=(c1,c2,,cm)即可.

显然, (3)(2)成立.

注2.1  对于任意ρS(A,B), 其中AS(H)BS(K), 存在一个交换算子F, 使得F(AB)=(BA)F(ρ)S(B,A). 这里F的具体作用结果如下: 若对于任意的量子态ρ, 其表示形式为

ρ=dH,dKi,j=1;k,l=1tij,klEijFkl,

那么

F(ρ)=dH,dKi,j=1;k,l=1tij,klFklEij.

推论2.2  设AS(H)BS(K),存在酉算子UV使得UAU=(˜AˇO), VBV=(˜BˆO), tr(A)=tr(˜A)=1tr(B)=tr(˜B)=1.˜A,˜B的秩分别为n,m, 其中n<, m<, 且n=rm. 那么下列三者等价

(1)存在ρS(A,B),其秩小于等于r.

(2)存在C=(Cij)Hn, 其中CijMm. 使得λ(CˆO)=λ(B)λ(C11ˇO+C22ˇO++CrrˇO)=λ(A).

(3)λ(C11)=λ(C22)==λ(Crr)(2)成立.

  根据注2.1交换算子的存在性, 我们只需要考虑rank(A)>rank(B)的情况即可.

(1)(2)(1)成立, 且˜A,˜B的秩满足关系n=rm, 可以参考有限维情形的证明[3]. 利用文献[3] 中的定理2.5可得, 存在C=(Cij)Hn, C1,,CrHm, 且˜C1,,˜CrHn, 使得˜Cj=Uj(CjO(r1)m)Uj, 这里Uj是一个n×n的酉算子. 取Ujm列元素构成Vj, 此时VjMn,m.Rj=VjC1/2jR=[R1||Rr], 则C=RR. 进而可以得出

(i)λ(˜Cj)=λ(CjOnm),(ii)˜A=rj=1Cj,(iii)˜B=rj=1˜Cj

以及

λ(C)=λ(˜B)=λ(RR)=λ(RR)=λ(rj=1˜Cj).

λ(C11++Crr)=λ(˜A)=λ(rj=1RjRj)=λ(rj=1Cj).

因此

λ(A)=λ(UAU)=λ(˜AˇO)=λ(C11ˇO+C22ˇO++CrrˇO)

λ(B)=λ(VBV)=λ(˜BˆO)=λ(CˆO)

C=RR满足条件(2).

(2)(1)(2)成立. 令C=RR,其中R=[R1||Rr],RjMn,m.Cii=RiRi, 显然λ(C)=λ(RR). 因此, 我们有

˜A=R1R1++RrRr,˜B=R1R1++RrRr,

λ(˜A)=λ(R1R1++RrRr)=λ(C11+C22++Crr),

λ(C)=λ(RR)=λ(R1R1++RrRr)=λ(˜B).

也就是说, λ(A)=λ(C11ˇO+C22ˇO++CrrˇO)λ(B)=λ(CˆO). 由于λ(RiRi)=λ(RiRiO(r1)m),因此利用定理2.2中前两者的等价性即可证明结论成立.

(2)(3)利用定理2.2的结论即可证明.

参考文献

Von Neumann J . Mathematical Foundations of Quantum Mechanics. United States: Princeton University Press, 1932

[本文引用: 2]

Klyachko A .

Quantum marginal problem and N-representability

J Phys Conf Ser, 2006, 36: 72- 86

DOI:10.1088/1742-6596/36/1/014      [本文引用: 1]

Li C K , Poon Y T , Wang X .

Ranks of quantum states with prescribed reduced states

Electron J Linear Algebra, 2017, 34: 331- 342

URL     [本文引用: 3]

Li C K , Poon Y T , Wang X .

Ranks and eigenvalues of states with prescribed reduced states

Electron J Linear Algebra, 2014, 27: 935- 950

URL     [本文引用: 3]

Fulton W .

Eigenvalues, invariant factors, highest weights, and Schubert calculus

Bull Amer Math Soc, 2000, 37: 209- 249

DOI:10.1090/S0273-0979-00-00865-X      [本文引用: 2]

Horn A .

Eigenvalues of sums of Hermitian matrices

Pac J Math, 1962, 12 (1): 225- 241

DOI:10.2140/pjm.1962.12.225     

Bercovici H , Li W S , Timotin D .

The Horn conjecture for sums of compact selfadjoint operators

Am J Math, 2009, 131: 1543- 1567

URL     [本文引用: 3]

/