周期驱动的二能级系统中的准宇称-时间对称动力学*

胡洲 曾招云 唐佳 罗小兵

1) (浙江理工大学物理系,杭州 310018)

2) (井冈山大学数理学院,吉安 343009)

本文研究了一个周期驱动的非宇称-时间对称二能级量子系统的非厄米动力学.通过经典相空间分析方法,解出了该非厄米系统的Floquet 态和准能谱,并解析构造了由该非厄米哈密顿量支配下的量子态的非幺正时间演化算符,给出了不同参数区域的量子态演化.本文数值和分析证明,该非宇称-时间对称二能级Floquet系统,类似于宇称-时间对称系统,存在一个准能谱从实数谱到复数谱的相变.本文还揭示了在量子态的动态演化中存在一种准宇称-时间对称动力学,即,该系统的粒子布居概率演化完全满足时间空间对称(宇称-时间对称),但是由于相位演化违反了宇称-时间对称性的要求,因此包含相位信息的量子态演化不满足时间空间对称(宇称-时间对称).这些结果加深了对非厄米物理的理解,拓展和推广了传统的宇称-时间对称概念.

1 引言

众所周知,量子体系的状态演化由哈密顿量确定并服从薛定谔方程.传统量子力学认为,描述一个量子力学系统的哈密顿量必须具有厄米性,才具有实数的能量本征值,从而保证能量的可观测性和体系概率守恒.1998 年,Bender和 Boettcher[1]发现,满足宇称-时间(parity-time,PT)对称性的非厄米哈密顿量的本征值也可以为实数.当非厄米参数强度增大时,PT 对称系统的能量本征值存在一个从实数谱到复数谱的对称破缺相变.尽管PT 对称量子理论源于量子力学,但随后被类比到经典物理系统中,在光学[2-6]、力学[7]、电学[8]、声学[9,10]等多种经典系统中被实验验证.但受制于量子力学的退相干性及量子实验系统的复杂性,真正的量子体系中实现宇称时间对称哈密顿量的演化具有巨大挑战.最近,人们在光子[11-13]、冷原子[14]、腔自旋波混合系统[15,16]、金刚石氮-空位色心[17]、核磁共振[18,19]、离子阱[20,21]等量子体系实验平台中成功构造了等效PT 对称量子哈密顿量.PT 对称物理的发展不但加深了人们对量子力学基本问题的理解,而且为单向隐身[10,22-24]、相干完美吸收[25]、非互易性传输[5,6,26]、超灵敏传感[27,28]、单模激光[29,30]、无线传能[31]、拓扑相的控制[4,32]等一系列应用技术带来了新的前景.

同时,PT 对称的非厄米量子系统相应的数学结构也吸引了人们的关注[33-37].实际上,算符厄米性只是本征值为实数的充分条件,而非必要条件.2002 年,Mostafazadeh[35,36]在数学上拓展了PT对称理论:对于一个任意哈密顿量H,若在线性厄米算符η的作用下,可以将该哈密顿量变成它的厄米共轭,即满足ηHη-1H†,则这个哈密顿量是赝厄米的,其本征值也可以为实数.这种赝厄米哈密顿量更为广泛,它包含了PT 对称哈密顿量.由于非厄米物理带来了丰富的新奇效应,构造和理解具有实本征谱的非PT 对称的非厄米哈密顿量也引起了人们的兴趣.例如,Nixon和Yang[38]最近在数学上构造了一类具有任意增益和耗散分布的复数势函数,证明了这类非厄米哈密顿量不仅具有实数谱,而且和PT 对称体系一样,能量本征谱存在一个从实数谱到复数谱的相变.这类非PT 对称势函数简洁而且容易调控,有望在多能级原子系统[39]等实验平台中实现.

对非厄米哈密顿量量子体系加以周期性调控,可以实现许多平衡态系统中无法达到的新奇的物理学现象,如非厄米Floquet 拓扑物态[12,40]、Floquet调控下的PT 相变[41-44]、Floquet 调控实现增益和双向隐身[45]等.周期调制的宇称-时间对称性物理已经在经典电路[42]、冷原子[14]、光子[12]等体系中实验实现.二能级系统(称为量子比特)是量子信息和量子计算的基本单元,PT 对称二能级非厄米量子系统自然引起了理论和实验的关注和兴趣.人们可以将宇称-时间对称哈密顿量扩展成更大空间中的厄米哈密顿量,从而实现量子比特子系统在非厄米哈密顿量支配下的演化.在前期工作中研究了周期性驱动对PT 对称二能级系统的量子调控[41,46-48],并获得了一些精确解[47,48],有助于对非厄米哈密顿量的精确操控.但是,目前对具有实能谱的非PT 对称二能级非厄米量子系统的Floquet调控的研究较少.

本文研究了一个周期驱动的非PT 对称二能级量子系统的非厄米动力学.通过经典相空间分析方法,首先得到了该非厄米系统的Floquet 态和准能谱,从而解析构造了非幺正时间演化算符,分析了不同参数区域量子态演化行为.通过解析和数值证明,该非PT 对称二能级Floquet 系统,类似于PT 对称系统,存在一个准能谱从实数谱到复数谱的相变.而且,在量子态的动态演化过程中,存在一种准PT 对称量子动力学:该系统的粒子布居概率演化满足时间空间对称(PT 对称),但是量子态演化(由于相位演化违反PT 对称性要求)不满足时间空间对称(PT 对称).这种准PT 对称量子动力学推广了传统的PT 对称概念,加深了对非厄米物理和非传统量子力学的理解.

2 模 型

考虑一个一般的周期驱动的非厄米二能级哈密顿量系统,其动力学演化满足的薛定谔方程可写成:

其中,态矢量|ψ〉(ψ1,ψ2)T(文中上标T表示矩阵转置),哈密顿量H具有如下形式

这里,γ0为两能级的静态失谐,γ1为非厄米参数,V(t)为两能级间的耦合强度,假定两能级的耦合强度具有时间周期性,即V(t+T′)V(t) .

首先介绍PT 对称系统的Floquet 本征模和准能量的基本概念和性质.若哈密顿量H(t) 与PT算符满足对易关系,即[H,PT]0,则该哈密顿量被认为具有PT 不变性,这里的为宇称算符,T是时间反演算符,具有如下效应i→-i,t0+t→t0-t这里,t0为某个时间反演点.可以证明:当γ00,且V(t0+t)V(t0-t) 时,哈密顿量(2)系统变成了标准的具有PT 对称性的二能级Floquet 系统,此时[H,PT]0 .与厄米Floquet 系统一样,周期驱动的非厄米(不限于PT)系统(1)同样具有Floquet 态解,即|ψ〉=|φ〉e-iεt,|φ(t+T′)〉|φ(t)〉在把时间作为另外一个维度的广义希尔伯特空间中,态矢量|φ〉(φ1,φ2)T满足如下本征值方程:

这里的F称为Floquet 哈密顿量,其定义为F:H -i∂/∂t,ε称为准能量,|φ〉称为Floquet 态本征模(或简称Floquet态).若[H,PT]0,则[F,PT]0也同样成立,该Floquet 系统具有PT 对称性.对于具有PT 对称性的Floquet 系统,如果Floquet态|φ〉满足本征值方程(3),那么PT |φ〉也是Floquet哈密顿量的本征模,具有复共轭的准能量ε*,即FPT |φ〉ε*PT |φ〉.当准能量为实数时,即εε*,PT算符和F算符具有共同的本征态,PT |φ〉eiβ|φ〉,β为任意相位(可以通过一个简单的规范变换去掉),此时系统具有稳定的本征动力学模,处于PT 对称未破缺相.当准能量为复数时,PT算符和F算符不再拥有共同的本征态,系统处于PT 对称破缺相.在PT 对称系统中,准能谱从实能谱到复能谱的转变,称为PT 对称破缺相变.

其次,简单讨论具有PT 对称性的二能级系统的一般量子态演化的性质.定义t0±t0±t,如果|ψ(t0±)〉[ψ1(t0+),ψ2(t0+)]T是薛定谔方程(1)的 解,那 么PT |ψ(t0+)〉[ψ2*(t0-),ψ1*(t0-)]T和[ψ1(t0+)e-iξ,ψ2(t0+)e-iξ]T(ξ为任意相位)也是薛定谔方程(1)的解.假如在t0时刻,ψ1(t0)ψ2*(t0)eiξ,则PT |ψ(t0+)〉和 e-iξ|ψ(t0+)〉是薛定谔方程(1)的同一个解,即:

方程(4)等价于

其物理含义是,在PT 对称二能级系统中,从t0时刻开始演化的量子动力学,沿着 +t方向演化和沿着-t方向演化在宇称算符(交换两个基矢,空间反演)操作下是完全相同的.

一个有趣的问题自然提出:是否存在非PT 对称的非厄米Floquet 系统,它具有类似PT 对称的性质,其准能谱也存在从实数谱到复数谱的相变? 如果存在,那么其一般量子态的动力学演化是否有可能也和PT 对称系统一样,仍然存在某种对称性呢?

本文以一个具体的Floquet 模型为例尝试去回答上述问题.在哈密顿模型(2)中,选择简单的驱动[49]:

这样的哈密顿量可以描述一个被外磁场B(Bx,By,Bz)驱动的二能级非厄米系统.在实验方面,以双势阱系统为例,这样一个依赖于时间的复数隧穿系数,可以同时调制两个势阱间的势垒高度和两个势阱的量子态的相角来实现.容易证明,对于任意时刻t0,等式V(t0+t)V(t0-t) 都不成立,即这个二能级非厄米Floquet 系统本质上不具有PT 对称性.在下一节,本文将求解该二能级非厄米系统的Floquet 态本征模和任意量子态的非幺正演化.

3 Floquet 态和准能谱

3.1 约化的有效哈密顿量

由于哈密顿量的非厄米性,该系统不再保持概率守恒,即量子态的模(范数)n(t)〈ψ(t)| ψ(t)〉|ψ1|2+|ψ2|2不再为常数,它随时间的演化可由非厄米薛定谔方程(1)得到:

因此,引入归一化的量子态:

该归一化的量子态满足新的薛定谔方程:

约化后的有效哈密顿量的形式为:

其中,σx,σy,σz为泡利矩阵的3 个分量.值得强调的是,虽然有效哈密顿量H′是非厄米的,但此时新的量子态的范数保持守恒,即〈ψ′(t)|ψ′(t)〉0 .

3.2 赝定点

以布居数差Z和相对相位θ′为两个独立变量,薛定谔方程(10)变成了如下的两个耦合方程:

现在我们可以在 (Z,θ′) 的相空间研究系统的动力学演化,在相空间定义一种“赝定点”,

这样一种“赝定点”可以认为是非含时哈密顿量系统的传统定点在Floquet 系统的一种拓展.它描述量子态在两能级的布居数差不随时间变化,但相对相位线性增长(每隔一个驱动周期相对相位增长2π,等价于每隔一个驱动周期T′2π/ω,相对相位不变).实际上,这样一种“赝定点”对应周期驱动量子系统的Floquet 态.为了方便研究这种“赝定点”,我们引入一个新的变量,θθ′-ωt,方程(13)变成:

首先,考虑共振(即两能级的失谐等于外场驱动频率,2γ0ω)情况的定点.为了方便计算,若无特殊说明,本文设定:当＜1时,sinα当＞1时,coshα.在共振情况下,方程(15)的定点可以分成以下三类情况:1)当＜1 时,存在两个定点 (Zf,θf)(0,2π-α),(Zf,θf)(0,α+π),这里 sinα;2)当1 时,只存在一个定点 (Zf,θf);3)当＞1 时,存在两个定点 (Zf,θf)这里 coshα

其次,考虑非共振(2γ0ω)情况.由方程(15)可知,非共振时,Zf0 .令可以得到关于Zf的四次方程,

求解方程(16),得到定点的布居数差:

这些定点解的稳定性可以通过计算雅可比矩阵:

的本征值λ1,2来分析.如果雅可比矩阵两个本征值的实部都小于或等于0,则这个定点是稳定的,否则,定点是不稳定的.在共振(2γ0ω)情况下,当＜1时,两个定点 (Zf,θf)(0,2π-α),(Zf,θf)(0,α+π)对应的雅可比矩阵的本征值为λ1,2雅可比矩阵的本征值为纯虚数,因此,这两个定点是稳定的.当1时,定点(Zf,θf)对应的雅可比矩阵的本征值λ0,对应临界稳定.当＞1时,两个定点(Zf,θf)对应的雅可比矩阵的本征值为λ1-2γ1Zf,λ2其中,定点对应的雅可比矩阵的两个本征值都为负的实数,为稳定的结点;另外一个定点对应的雅可比矩阵的两个本征值都为正的实数,为不稳定的结点.同理,非共振(2γ0ω)情况的雅可比矩阵分析可知,对于任意大小的非厄米参数γ1,一个定点是稳定的,一个定点是非稳定的.

图1 给出了 (Z,θ) 的经典相图.图1(a)—(c)给出了共振(2γ0ω)的3 种情况,图1(d)给出了非共振(2γ0ω)情况.图中的红点表示定点.文中所有参数都已经无量纲化.

图1 不同系统参数下的相空间轨道(ω=2,ν=1) (a)γ0=1,γ1=0.5;(b) γ0=1,γ1=1;(c) γ0=1,γ1=2 ;(d)γ0=2,γ1=0.5Fig.1.Phase-space trajectories with different system parameters(ω=2,ν=1):(a) γ0=1,γ1=0.5;(b)γ0=1,γ1=1;(c) γ0=1,γ1=2;(d) γ0=2,γ1=0.5 .

3.3 Floquet 态和准能量

本节将从相空间的赝定点构造系统的Floquet态和准能量.约化后的薛定谔方程(10)存在归一化的Floquet 态解|ψ′〉,|φ′(t+T′)〉|φ′(t)〉,这里,T′2π/ω,量子态的归一化条件要求ε′为实数.Floquet 态|φ′〉满足本征值方程(H′-i∂t)|φ′〉ε′|φ′〉.由定点 (Zf,θf) 可以构造出归一化的Floquet 态解|ψ′〉:

显然,态向量中两组分的布居数差为Zf,相位差为θ2-θ1θf+ωt,对应经典相空间的赝定点.式中的相位-mωt反映了准能谱的类似布里渊区的周期性结构.本文中如无特殊说明取m0,把准能量限制在一个布里渊区内.

由(9)式可以得到原始哈密顿量的Floquet 态解:

此处用到关系式n〈φ| φ〉e2Im(ε)t,εRe(ε)+i Im(ε).(20)式和原始Floquet 态解|ψ〉|φ〉e-iεt,相比较,可得到具有时间周期性的Floquet 态本征模以及ε′Re(ε),也就是说,归一化的Floquet 态解的准能量是原始Floquet 态解的准能量的实部,其值求解如下:

由(8)式可知,原始哈密顿量Floquet 态解的模(范数)随时间的演化方程满足,

式中的n(0)表示初始时刻t0 的原始哈密顿量Floquet 态解的模(范数).(22)式与n〈φ| φ〉e2Im(ε)t比较,可以得到原始Floquet 态解的准能量的虚部,

以及,

假如设定n(0)1,则原始哈密顿量的Floquet态解(20)可以简化为:

准能量为:

由(26)式可知,如果Floquet 态在两能级上的占有概率(布居数)相等,Zf0,那么准能谱是全实数谱,如果两能级上的占有概率(布居数)不相等,Zf0,准能谱就会变成复数谱.

本文所考虑的周期驱动非厄米系统,在归一化的相空间中具有特殊的“赝定点”,它对应不依赖时间的布居数差和依赖于时间的相对相位.但是对于更一般的周期驱动非厄米系统,在归一化的经典相空间中,Floquet 态本征模应该对应庞加莱截面(每隔驱动周期间隔采样)的定点.即,布居数差一般是时间的周期函数,ZZd+Za(t),这里,Zd是常数代表庞加莱截面的定点位置,Za(t) 是平均值为0 的时间周期函数.这时,Floquet态解的模(范数)随时间的演化方程满足n(t)与n〈φ|φ〉e2Im(ε)t比较,可以得到,准能量的虚部为,Im(ε)γ1Zd,以及,这样,原始哈密顿量的概率不守恒的Floquet 态解可以构造出来,即:|ψ〉.可以知道,原始哈密顿量Floquet 态解的准能量的虚部由Zd决定.

结合(19),(25),(26)式,可以确定薛定谔方程(1)的Floquet 态和准能量.考虑非共振(2γ0ω)情况,其定点解(17)Zf0,结合给出准能量(m0)为:

显然,在非共振情况下,由于Zf0,对于任意非零的非厄米强度γ1,其准能量都是复数.

接下来,本文侧重于研究共振(2γ0ω)情况的本征模和量子动力学演化,将分别给出共振情况三类定点对应的Floquet 态和准能量.＜1 .这时,方程(15)有两个布居数对称的定点.定 点 (Zf,θf)(0,2π-α) 对应的Floquet 态解为:

准能量为:

另外一个定点 (Zf,θf)(0,α+π) 对应的Floquet态解为:

准能量为:

准能量为ε-γ0.这时,两个Floquet 态本征模塌缩成一个,这个参数点被称为奇异点(exceptional point).系统能谱在奇异点附近有独特的拓扑结构,将带来许多新奇的物理现象和应用,如增强传感灵敏度[27,28,50],非对称模式转化和手性传输[4]等.

准能量为:

准能量为:

图2 给出了非厄米系统(1)在周期驱动(7)作用下共振(2γ0ω)情况时的准能量的实部(见图2(a))和虚部(见图2(b))随参数γ1/ν的关系.红线和蓝线代表直接对角化一个驱动周期的时间演化算符的数值结果,圆圈代表经典相图中赝定点对应的准能量解析结果.如图所示,数值结果和解析结果完全符合.尽管在周期驱动(7)作用下,V(t0+t)V(t0-t),从而非厄米系统(1)不具有PT 对称性,[H,PT]0,但可发现,它和PT 对称系统一样,存在一个准能谱从实数谱到复数谱的相变.当非厄米系统(1)的准能谱为实数时,将PT算符作用在实数准能谱ε-(γ0-νcosα) 对应的Floquet 态本征模上,可以验证,对于任意的时刻点t0,

图2 非厄米系统(1)在周期驱动(7)作用下共振(2γ0=ω)情况时的准能谱的实部(a)和虚部(b)随参数γ1/ν 的关系.红线和蓝线代表直接对角化一个驱动周期的时间演化算符的数值结果,圆圈代表经典相图中赝定点对应的准能量解析结果.参数取为γ0=1,ω=2Fig.2.Real (a) and imaginary (b) parts of the quasienergies as a function of γ1/ν for the non-Hermitian system (1)subject to a periodic modulation (7) in the resonant (2γ0=ω) case.The red and blue lines denote the numerical results of quasienergies computed through direct diagonalization of the time-evolution operator over one period of the driving,while the circles denote exact analytical results of quasienergies corresponding to the pseudo fixed points in phase space.The system parameters are set as γ0=1,ω=2.

式中,C为任意常数.这说明,该实数准能谱对应的Floquet 态|φ(t0+t)〉不是PT算符的本征模.同理可证,另外一支实数准能谱ε-(γ0+νcosα)对应的Floquet 态也不是PT算符的本征模.这进一步说明,在周期驱动(7)作用下的非厄米Floquet系统(1)具有实数准能谱,但该系统不具有PT 对称性.

4 量子态的非幺正动力学演化

在第3 节中,已经通过相空间中的定点方法解析构造了系统的Floquet 态本征模和准能量.由这些Floquet 态和准能量,我们可以求解任意量子态的时间演化.在共振情况下,即使本文所研究的Floquet 系统不具有PT 对称性,也存在能谱从实数谱到复数谱的相变.那么,一般量子态的动力学演化又会具有什么独特的量子性质呢? 分3 种情况求解在共振(2γ0ω)情况下的一般量子态的动力学演化.

式中叠加系数C1,C2由初态决定.

假设给定t0时刻的态矢量,|ψ(0)〉[ψ1(0),ψ2(0)]T,时间演化算符作用在该态矢量上,可得:

由初态|ψ(0)〉决定(38)式的叠加系数C1,C2后,比较(38)和(39)式,可以写出从t0 开始的非幺正时间演化算符的具体的矩阵形式:

以ψ1(0)1,ψ2(0)0 为初态直接数值模拟原始薛定谔方程(1),图3 给出了在共振情况下,当＜1时的系统动力学.从图3 可以看出,两能级上的占有概率|ψ1|2和|ψ2|2随时间周期振荡,当(k为整数)时,在两能级上的占有概率相等,|ψ1|2|ψ2|2

.图中以t0-0.605为例说明了这个特性.图3 显示,从t0时刻开始,沿着正向演化(+t方向)和沿着反向演化(-t方向)在宇称算符(交换两个基矢|1〉和|2〉,空间反演)操作下,概率是完全相同的,即|ψ1(t0+t)|2|ψ2(t0-t)|2.但相位之和,θ1(t0+t)+θ2(t0-t),却不等于t0时刻的值ξθ1(t0)+θ2(t0),反而随时间t线性变化.注意,在数值中,相位每隔2π 有个跳跃.图中的数值结果完全符合理论预言.

图3 在共振情况下的系统动力学(＜1)(a)两能级 上的占有概率|ψn|2(n=1,2)随时间演化;(b)相位θn (n=1,2)随时间演化;(c)以 t0=-0.605 为时间反演点的θ1(t0+t)+θ2(t0-t)演化.系统参数取为γ0=1,γ1=0.5,ω=2,ν=1初态为ψ1(0)=1,ψ2(0)=0Fig.3.System dynamics for the resonance case with the non-Hermitian parameters ＜1,starting the system with the state ψ1(0)=1,ψ2(0)=0 :(a) Time evolutions of the occupation probabilities|ψ1|2and|ψ2|2 ;(b) time evolutions of phases θ1(t)and θ2(t) ;(c) time evolution of the sum of phases,θ1(t0+t)+θ2(t0-t) .Here we choose the time-inversion point t0=-0.605 .The system parameters are γ0=1,γ1=0.5,ω=2,ν=1 .

2)γ1/ν1 .这时,系统处于奇异点,只有一个Floquet 态解(32),但可以构造出另外一个线性无关的基本解:

因此,由这两个线性无关的基本解,可以写出薛定谔方程(1)的通解:

式中叠加系数C1,C2由初态决定.由(44)和(45)式,可以得到,从t0 开始的非幺正时间演化算符为

当t0时刻的态制备为ψ1(0)1,ψ2(0)0 时,把时间演化算符(46)作用在t0 的态矢量上,可得:

可以证明,当t0-时,|ψ1(t0+t)||ψ2(t0-t)|,θ1(t0+t)+θ2(t0-t)π/2-ωt.

取ψ1(0)1,ψ2(0)0 为初态,通过直接数值积分模拟原始薛定谔方程(1),图4 给出了在共振情况下,当γ1/ν1 时的系统动力学演化.从图4可以看出,两能级上的占有概率从t0-0.5 开始沿着正向演化(+t方向),呈时间t的平方增长.如理论预期一样,数值结果显示,|ψ1(t0+t)|2|ψ2(t0-t)|2,相位之和,θ1(t0+t)+θ2(t0-t),却不恒等于t0时刻的相位之和,而是随时间t线性减小.注意,在数值中,相位每隔 2π 有个跳跃.

图4 在共振情况下的系统动力学(γ1/ν=1) (a) 两能级上的占有概率|ψn|2 (n=1,2)随时间演化;(b)相位 θn (n=1,2)随时间演化;(c)以 t0=-0.5 为时间反演点的θ1(t0+t)+θ2(t0-t)演化.系统参数取为γ0=1,γ1=1,ω=2,ν=1初态为ψ1(0)=1,ψ2(0)=0Fig.4.System dynamics for the resonance case with the non-Hermitian parameters γ1/ν=1,starting the system with the state ψ1(0)=1,ψ2(0)=0 :(a) Time evolutions of the occupation probabilities|ψ1|2and|ψ2|2 ;(b) time evolutions of phases θ1(t)and θ2(t) ;(c) time evolution of the sum of phases,θ1(t0+t)+θ2(t0-t) .Here the timeinversion point is given by t0=-0.5 .The system parameters are γ0=1,γ1=1,ω=2,ν=1 .

3)γ1/ν ＞1 .由系统的两个基本的Floquet 态解(33)和(35),可以构造在这种情况下的通解

式中叠加系数C1,C2由初态决定.同理,由(49)和(50)式,可以得到,从t0 开始的非幺正时间演化算符为

当t0时刻的态制备为ψ1(0)1,ψ2(0)0时,把时间演化算符(51)作用在该态矢量上,可得

图5 给出了在共振情况下,当γ1/ν ＞1 时的系统动力学演化.同样地,这里也是以ψ1(0)1,ψ2(0)0为初态直接数值模拟了原始薛定谔方程(1).如图5 所示,两能级上的占有概率从t0-0.380 开始沿着正向演化(+t方向),呈指数形式增长.类似地,如理论预期的一样,数值结果显示,|ψ1(t0+t)|2|ψ2(t0-t)|2,相位之和,θ1(t0+t)+θ2(t0-t),随时间t线性减小.

以上分析和数值证明,具有周期驱动形式(7)的系统(1)虽然不具有PT 对称性,但在共振情况下,其量子态演化具有类似PT 对称的性质.当非厄米参数强度小于某个临界值时,两能级系统布居概率呈周期振荡;当非厄米参数强度大于该临界值时,布居概率呈指数增强.而且,在共振情况下,从t0时刻开始演化的量子态,沿着正方向(+t方向)演化和沿着反方向(-t方向)演化在宇称算符(交换两个基矢,空间反演)操作下,概率是完全相同的.也就是说,从系统粒子布居概率来看,该非厄米系统的动力学和PT 对称系统完全一样.但实际上,两能级系统的概率幅的相位不满足PT 对称性的要求.从数学上讲,该系统只满足PT 对称的其中一个条件:|ψ1(t0+t)||ψ2(t0-t)|,但不满足另外一个条件θ1(t0+t)+θ2(t0-t)ξ,这里ξ是个常数,由t0时刻的概率幅的相位决定.我们把这种类似PT 对称性的量子动力学称为准PT对称量子动力学.

5 结论

采用经典相空间分析方法研究了一个周期驱动的非PT 对称二能级量子系统的非厄米物理特性.由于哈密顿量的非厄米性,体系概率不守恒,引入一个归一化的量子态,并把薛定谔方程转化成了一个以两能级布居数差和相对相位为两个独立变量的经典动力学方程.研究发现,该量子系统的Floquet 态对应经典相空间的“赝定点”,即不随时间变化的布居数差和依赖于时间的相对相位.在经典相空间“赝定点”的基础上,我们解析构造了由该非厄米哈密顿量支配下的量子态的非幺正时间演化算符,给出了不同参数区域的量子态演化.研究结果表明,该非PT 对称二能级量子系统也存在类似PT 对称系统一样的从对称未破缺区到破缺区的过渡,而且在该系统量子态的动力学演化中,存在一种准PT 对称量子动力学,即系统粒子布居概率演化完全满足时间空间对称(PT 对称),但由于相位演化不满足PT 对称性的要求,因此,该量子态演化不满足时间空间对称(PT 对称).这些研究结果有利于进一步加深对非厄米量子系统,特别是对PT 对称概念的理解.

对于更一般的周期驱动非厄米系统,在归一化的经典相空间中,Floquet 态本征模对应庞加莱截面(每隔驱动周期间隔采样)的定点.本文建立的经典相空间定点同Floquet 态和准能量的数学联系,可以推广到非线性的非厄米量子系统中,有助于研究开放量子多体系统的经典和量子对应,有助于构造新的超越标准量子极限的快速演化非幺正量子门[18,51],这些都是未来进一步的研究工作的重点.