摘要:本文以具有非对称性自旋相互作用的三体自旋系统为研究对象, 重点研究了三体量子相干含时演化规律. 采用精确量子对角化和基于量子主方程的数值模拟方法, 讨论了三体量子系统中多种量子相干组分及其退相干. 研究发现, 量子相干组分的含时演化与整个系统的初态量子特性紧密相关. 当初态为可分离纯态时, 在较短时间内, 非对称相互作用有利于增加多体量子相干度. 这些量子相干度因受噪声影响而逐渐衰减. 当初态为类Werner态时, 量子相干度的分布满足加和性, 即三体量子相干度等于所有两体量子相干度之和. 自旋之间非对称相互作用和环境噪声都会引起三体量子相干度大于所有两体量子相干度之和. 这些结论有助于多体量子资源的制备.
关键词: 多体量子相干/
非对称自旋系统/
量子退相干/
量子主方程

English Abstract

全文HTML

--> --> -->

作为重要的量子资源, 量子相干^[1-3]已被广泛应用在量子信息处理^[4,5]、量子测量^[6]和生物物理^[7]等领域. 根据量子态叠加原理, 量子相干可以表征物理系统的非经典特性. 一方面, 由于量子相干可以是单个系统中不同能态的叠加, 人们通过实验可以观测到不同能态的振荡现象. 另一方面, 量子相干也包含在多个系统之间, 表现出一种非经典关联, 与量子纠缠和量子失协都存在紧密联系. 在量子光学领域中, 人们通常用光场相空间分布和多点关联函数来描述光场的量子相干性. 利用光子偏振态的相干叠加, 人们可以观测到量子干涉条纹, 为实验测量量子相干性提供了一种有效方法. 随着量子调控技术的发展, 人们对量子相干开展了丰富的理论和实验研究, 研究对象从量子光学系统扩展到固态自旋系统和光学混合量子系统^[8,9]等. 这些量子相干系统已经从单体二维封闭系统^[10]发展到多体高维开放系统^[11]. 同时, 人们将量子相干资源与量子热力学^[12]、量子相变^[13]、量子退相干^[14]和量子耗散^[15]等理论紧密相联, 利用量子相干度阐述了一些重要量子现象的产生机理. 为了定量评价物理系统的量子相干性, 2014年Baumgratz等^[16]基于量子资源理论提出了一些有效的量子相干度量方法, 和作为量子相干度量应满足的必要条件. 一些常用判据被用于量子相干度量, 如基于${l_1}$范数^[17,18]的量子相干度量、基于相对熵^[19,20]的量子相干度量和基于其他量子态距离判据^[21]的度量方法.
目前, 人们非常关注多体系统的量子相干分布特性^[22]和动力学特性^[23]. 一方面, 这些量子相干特性可以用来解释多体系统的集合特性(如量子相变); 另一方面, 不同量子相干组分的含时演化可以帮助人们更好地理解量子资源在时空的分布规律, 为发展量子网络通信提供必要的量子资源. 于是, 一些典型量子多体系统, 如量子自旋系统^[24-26]、光学混合量子系统、甚至一些生物系统^[27]都成为了多体量子相干性的研究对象. 其中, 海森堡自旋系统^[28]的量子相干就是一种具有代表性的多体量子相干. 最近, 人们讨论了量子相干对基于海森堡XXX自旋工作物质的量子热机效率的影响^[29]. 在海森堡自旋系统中, 人们分析了不确定性关系及量子相干性的联系^[30], 定量评价了海森堡XY自旋系统的量子相干非局域优势及量子相变表征^[31]. 在这些研究工作中, 人们都考虑了由于自旋轨道耦合引起的非对称Dzyaloshinsky-Moriya (DM)相互作用^[32,33]. 然而, 任何系统都会不可避免地与周围环境发生相互作用, 这种作用将引发量子退相干和量子耗散. 于是, 我们以具有非对称自旋DM相互作用的多体量子系统为研究对象, 分析多体量子相干组分的含时演化规律, 讨论环境噪声对多体量子相干组分的影响, 揭示量子相干的时空分布特性.
本文借助精确对角化^[34]和量子主方程数值计算方法, 定量研究了具有非对称自旋-轨道耦合作用的三体系统量子相干. 研究发现, 存在一类多体量子相干态, 其整体量子相干度与任意两组分量子相干度之和具有一种紧密联系. 首先, 介绍了具有非对称相互作用的量子自旋系统模型, 通过对三体系统的精确对角化, 解析得到了量子相干态的含时演化规律; 基于${l_1}$范数定量分析了多种组分与整体量子相干度之间的关系, 证实了量子相干加和性的存在条件. 其次, 采用量子主方程数值计算方法, 研究了退振幅型量子噪声对多体量子相干组分的影响. 最后, 得到了一些有效结论, 这些研究结果有助于规模化量子相干的制备, 为量子通信网络提供必要的量子资源.

本文考虑一种具有自旋-轨道耦合相互作用的多体非对称量子自旋系统. 在半导体量子点和分子磁体中, 自旋-轨道耦合产生的相互作用是一种典型的非对称作用, 通常描述导带电子的相互作用形式^[35-37]. 这种模型可以用来研究固态量子系统的磁性特征和量子临界现象等^[38]. 人们也可以利用微腔耦合系统模拟实现该模型, 用于量子态远程传输和量子通信网络^[39,40]. 在外磁场中, 系统模型的哈密顿量表示为

$\begin{split}H =\;& J\sum\limits_{i = 1}^{N - 1} \frac{{1 + {{( - 1)}^{i + 1}}\gamma }}{2} \cdot \left( {\sigma _i^ + \sigma _{i + 1}^ - + \sigma _i^ - \sigma _{i + 1}^ + } \right) \\&+ {{D}} \cdot \left( {{{{\sigma}} _i} \times {{{\sigma}} _{i + 1}}} \right) + h\sum\limits_{i = 1}^N {\sigma _i^z},\end{split}$

其中J表示自旋之间量子隧穿的相互作用强度, $\gamma $为相邻自旋之间各项异性耦合参数, D为自旋与轨道耦合的非对称参量, h为沿z方向的外磁场强度, $\sigma _i^ \pm, {{{\sigma}} _i}$为第i个自旋的上升(下降)算子和Pauli自旋算符. 其中, 每个自旋态满足${\sigma ^z}$ $\left| {1(0)} \right\rangle = \pm \left| {1(0)} \right\rangle $. 考虑一种简单情况$N = 3$, ${{D}} = D{{{e}}_z}$. 通过对三体自旋系统的研究, 可以得到一些多体量子相干的空间分布特性.
首先, 通过精确对角化, 在自旋表象${\big\{ }\left| {111} \right\rangle, \left| {110} \right\rangle, \cdots \left| {000} \right\rangle {\big\}}$中, 得到(1)式的能量本征值及其相应本征态:

$\begin{split}&{E_1} = 3h,\\&\left| {{\varphi _1}} \right\rangle = \left| {111} \right\rangle ,\\&{E_{2,3}} = h \pm J\sqrt {8{D^2} + 2(1 + {\gamma ^2})} ,\\&\left| {{\phi _{2,3}}} \right\rangle = A\left| {110} \right\rangle \pm ({R_1} \!+\! {\rm{i}}{I_1})\left| {101} \right\rangle \!+\! ({R_2} \!+\! {\rm{i}}{I_2})\left| {011} \right\rangle ,\\&{E_{4,5}} = \pm h,\\&|{\varphi }_{4}\rangle =B(|110\rangle \!+\!K|011\rangle ), \left| {{\varphi _5}} \right\rangle \!=\! A(\left| {110} \right\rangle \!+\! {K^ * }\left| {011} \right\rangle ),\\&{E_{6,7}} = - h \pm J\sqrt {8{D^2} + 2(1 + {\gamma ^2})},\\&\left| {{\phi _{6,7}}} \right\rangle = A\left| {001} \right\rangle \!\pm \!({R_1} \!-\! {\rm{i}}{I_1})\left| {010} \right\rangle \!+\! ({R_2} \!-\! {\rm{i}}{I_2})\left| {100} \right\rangle ,\\&{E_8} = - 3h,\\&\left| {{\varphi _8}} \right\rangle = \left| {000} \right\rangle ,\\[-12pt]\end{split}$

其中, 参量

$\begin{split} &A = - \dfrac{1}{2}\sqrt {\dfrac{{4{D^2} + {{(1 - \gamma )}^2}}}{{4{D^2} + (1 + {\gamma ^2})}}},\\ &{I_1} = \dfrac{{ \pm \sqrt 2 D}}{{\sqrt {4{D^2} + {{(1 - \gamma )}^2}} }},~~ {R_1} = - \dfrac{{(1 - \gamma )}}{{2 D}}{I_1}\\ & B = \dfrac{{4{D^2} + {{(1 + \gamma )}^2}}}{{\sqrt {{{\left[ {4{D^2} + {{(1 + \gamma )}^2}} \right]}^2} + \left[ {{{(4{D^2} + {\gamma ^2} - 1)}^2} + 16{D^2}} \right]} }}\\ &{I_2} = \dfrac{{2 D}}{{\sqrt {\left[ {4{D^2} + \left( {1 + {\gamma ^2}} \right)} \right]\left[ {4{D^2} + {{\left( {1 - \gamma } \right)}^2}} \right]} }},\\ &{R_2} = \dfrac{{4{D^2} - {{(1 - \gamma )}^2}}}{{8 D}}{I_2},\\ &K = \dfrac{{4{D^2} + {\gamma ^2} - 1 + 4{\rm{i}}D}}{{4{D^2} + {{(1 + \gamma )}^2}}}.\end{split}$

经过含时演化算符$U(t) = \sum\nolimits_j^{} {\exp ( - {\rm{i}}{E_j}t)} \left| {{\phi _j}} \right\rangle \left\langle {{\phi _j}} \right|$作用, 这个多体系统在任意t时刻的量子态密度矩阵表示为

${{\rho}}(t) = {{U}}(t){{\rho}} (0){{{U}}^ + }(t).$

这里, ${{\rho}} (0)$表示系统初态.
选取两种初态: 类Werner态和类GHZ态. 其中, 整个系统的基态符合类Werner态特征, 是一种典型的非经典态. 首先, 以类Werner态为系统初态, 即$\rho (0) = \left| \varphi \right\rangle \left\langle \varphi \right|, \left| \varphi \right\rangle = a\left| {001} \right\rangle \!+\! b\left| {010} \right\rangle \!+\! c\left| {100} \right\rangle$(${\left| a \right|^2} + {\left| b \right|^2} + {\left| c \right|^2} = 1$). 由于系统哈密顿量满足对易关系$\left[H, \displaystyle\sum\nolimits_i {\sigma _i^z} \right] = 0$, 所以此时系统态的演化发生在希尔伯特子空间$\{ \left| {001} \right\rangle, \left| {010} \right\rangle, \left| {100} \right\rangle \} $里, 任意t时刻量子态满足$\rho (t) = \left| {\psi (t)} \right\rangle \left\langle {\psi (t)} \right|$. 利用(2)式和(3)式得到,

$ \begin{split}&\left| {\psi (t)} \right\rangle = \left( {\begin{aligned} {{y_1}} \\ {{y_2}} \\ {{y_3}} \end{aligned}} \right) = {{S}}\left( {\begin{aligned} a \\ b \\ c \end{aligned}} \right), \\&{{S}} = \left( {\begin{array}{*{20}{c}} {{S_{11}}}&{{S_{12}}}&{{S_{13}}} \\ {S_{12}^ * }&{{S_{22}}}&{{S_{23}}} \\ {S_{13}^ * }&{S_{23}^ * }&{{S_{33}}} \end{array}} \right), \end{split}$

这里的幺正操作S,

$\begin{split}& {S_{11}} = F{A^2} + \exp ( - {\rm{i}}{E_4}t) {B^2}, \\& {S_{12}} = FA({R_1} - {\rm{i}}{I_1}), \\& {S_{13}} = FA({R_2} - {\rm{i}}{I_2}) + \exp ( - {\rm{i}}{E_4}t){A^2}{K^*}, \\& {S_{22}} = F(R_1^2 + I_1^2), \\& {S_{23}} = F({R_1} + {\rm{i}}{I_1}) ({R_2} - {\rm{i}}{I_2}), \\& {S_{33}} = F(R_2^2 + I_2^2) + \exp ( - {\rm{i}}{E_4}t){B^2}{| K |^2}, \\& F = \exp ( - {\rm{i}}{E_2}t) + \exp ( - {\rm{i}}{E_3}t). \end{split}$

为了研究多体量子相干的空间分布特性, 分别计算了任意两体系统量子态的约化密度矩阵及其量子相干度. 采用${l_1}$相干性范数来计算量子相干,

$C(\rho ) = \sum\limits_{ i,j ~~ i \ne j} {\left| {{\rho _{ij}}} \right|} ,$

即密度矩阵的非对角元素绝对值之和. Baumgratz等^[16]基于量子资源理论, 提出了一种便于计算度量的量子${l_1}$相干性范数判据. 对于d维希尔伯特空间, 选取固定基矢$\{ \left| i \right\rangle \} $来表示量子态, 人们将$\rho =\displaystyle \sum\nolimits_i {{p_i}\left| i \right\rangle \left\langle i \right|}$定义为非相干态. 密度矩阵的非对角元素就包含了系统的量子相干性. 相比其他判据, 这种相干性范数不仅更加直观, 而且满足量子相干判据的基本条件, 如非负性、单调性和凸函数性等性质. 利用(5)式, 得到任意两组分的量子相干, ${C_{12}} = 2\left| {{y_2}y_3^ * } \right|$, ${C_{23}} = 2\left| {{y_1}y_2^ * } \right|$, ${C_{13}} = 2\left| {{y_1}y_3^ * } \right|$. 同时, 当系统初态为类Werner态时, 单个自旋的量子相干为零, 并且三体量子相干表示为

${C_{123}} = {C_{12}} + {C_{23}} + {C_{13}}.$

(6)式表明: 对于类Werner态, 系统在非对称性相互作用下, 三体量子相干度与任意两体量子相干度之间满足一种加和性. 对于一般情况$N > 3$, 根据(1)式描述的哈密顿量特性$\left[H, \displaystyle\sum\nolimits_i {\sigma _i^z} \right] \!=\! 0$, 可以得到某个N体类Werner态, $\left| {{\varphi _N}} \right\rangle \!=\! \displaystyle\sum\nolimits_{i = 1}^N {{a_i}| {{1_i}} \rangle }$, 其中$\left| {{1_i}} \right\rangle = \left| {00 \cdots 1 \cdots 0} \right\rangle $表示第i个自旋格点为$\left| 1 \right\rangle $且其余格点为$\left| 0 \right\rangle $. 利用(5)式的相干性范数判据, 可以得到整个多体系统的量子相干度${C_{{\rm{tot}}}} = \displaystyle \sum\nolimits_{i \ne j} {2\left| {{a_i}a_j^ * } \right|}$. 任意两体量子相干度与整体系统量子相干度也满足这种加和性, 即${C_{{\rm{tot}}}} = \displaystyle\sum\nolimits_{i \ne j} {{C_{ij}}}$. 以$N = 4$为例, ${C_{12}} \!=\! 2\left| {{a_3}a_4^ * } \right|, {C_{13}} \!=\! 2\left| {{a_2}a_4^ * } \right|,$ ${C_{14}} = 2 \left| {{a_2}a_3^ * } \right|, {C_{23}} = 2\left| {{a_1}a_4^ * } \right|$, ${C_{34}} = 2\left| {{a_1}a_2^*} \right|$. 然而, 当初态选取类GHZ态, $\left| \varphi \right\rangle = \sqrt q \left| {000} \right\rangle + \sqrt {1 - q} \left| {111} \right\rangle, (0 < q < 1)$三体量子相干度$C_{123}^{{\rm{GHZ}}} = 2\sqrt {q(1 - q)} \cos (6 ht)$, 任意两组分的量子相干度总为零. 这个结论进一步表明: 类GHZ态是一种不可区分的真实多体量子纠缠态, 任意两组分之间量子关联为零. 多体量子相干度的空间分布依赖于不同初态的性质.

在实际情况下, 系统总是与周围环境存在相互作用, 其量子相干特性往往在环境影响下发生量子退相干. 接下来考虑一种典型的量子退振幅噪声, 采用量子主方程和量子Kraus算符, 讨论噪声因素对多体量子相干度的影响.
关于非对称自旋系统, 合理考虑每个量子自旋i与其定域性环境${E_i}$存在弱耦合作用, 这种相互作用可以用量子退振幅噪声模型来描述^[41]. 为了有效描述这种量子噪声模型, 人们通常采用量子主方程和量子信道来研究噪声对量子系统含时演化的影响. 自旋系统与周围环境初态为${\rho _{SE}} = \rho (0) \otimes ({\left| {000} \right\rangle _E}\left\langle {000} \right|)$, 其中${\left| {000} \right\rangle _E} = {\left| 0 \right\rangle _{{E_1}}} \otimes {\left| 0 \right\rangle _{{E_2}}} \otimes {\left| 0 \right\rangle _{{E_3}}}$表示三个定域性环境初态为真空态, $\rho (0)$表示三个自旋系统初态. 在弱耦合近似条件下, 整个自旋系统状态演化可以通过量子主方程描述:

$\dot \rho = - {\rm{i}}\left[ {H,\rho } \right] + L(\rho )$

其中, 超算符Lindblad算符表示为

$\begin{split} L(\rho ) =\;& \displaystyle\sum\limits_i {{\varGamma _i}(n + 1)(2\sigma _i^ - \rho \sigma _i^ + - \rho \sigma _i^ + \sigma _i^ - - \sigma _i^ + \sigma _i^ - \rho )} \\&+{ \varGamma _i}n(2\sigma _i^ + \rho \sigma _i^ - - \rho \sigma _i^ - \sigma _i^ + - \sigma _i^ - \sigma _i^ + \rho ).\\[-10pt] \end{split} $

这里, n表示有限温度环境的平均光子数, 当温度越高时, n的数值越大. ${\varGamma _i} = \varGamma $表示每个自旋自发衰减的阻尼系数. 这些噪声因素会影响自旋系统的演化.
为了说明噪声对多体量子相干度的影响, 首先考虑一种最简单情况, $J = 0, \begin{array}{*{20}{c}} {h = 0} \end{array}$, 即自旋间相互作用为零. 整个量子态演化可以由量子Kraus算符表示为

$\rho (t) = \sum\limits_{i = 1}^3 {\sum\limits_{j = 0}^3 {{K_{ij}} \cdot \rho (0)} } \cdot K_{ij}^ + ,$

其中, 关于第i个自旋演化的量子Kraus算符分别是:

$\begin{split}&{K_{i0}} = \sqrt {\frac{{n + 1}}{{2n + 1}}} \left( {{{\left| 0 \right\rangle }_i}\left\langle 0 \right| + \sqrt {1 - p} {{\left| 1 \right\rangle }_i}\left\langle 1 \right|} \right),\\&{K_{i1}} = \sqrt {\frac{{n + 1}}{{2n + 1}}p} {\left| 0 \right\rangle _i}\left\langle 1 \right|,\end{split}$

$\begin{split}&{K_{i2}} = \sqrt {\frac{n}{{2n + 1}}} \left( {\sqrt {1 - p} {{\left| 0 \right\rangle }_i}\left\langle 0 \right| + {{\left| 1 \right\rangle }_i}\left\langle 1 \right|} \right), \\&{K_{i3}} = \sqrt {\frac{{np}}{{2n + 1}}} {\left| 1 \right\rangle _i}\left\langle 0 \right|.\end{split}$

利用Kraus算符, 得到任意自旋演化过程,

$ \begin{split}&\displaystyle\sum\limits_{j = 0}^3 {K_{ij}}\left( {{{\left| 0 \right\rangle }_i}\left\langle 0 \right|} \right)K_{ij}^ + = \left( {1 - \frac{{np}}{{2n + 1}}} \right){{\left| 0 \right\rangle }_i}\left\langle 0 \right| \\&\qquad\qquad\qquad\qquad\qquad+ \frac{{np}}{{2n + 1}}{{\left| 1 \right\rangle }_i}\left\langle 1 \right|, \\&\sum\limits_{j = 0}^3 {K_{ij}}\left( {{{\left| 1 \right\rangle }_i}\left\langle 1 \right|} \right)K_{ij}^ + = \left( {1 - \frac{{n + 1}}{{2n + 1}}p} \right){{\left| 1 \right\rangle }_i}\left\langle 1 \right| \\&\qquad\qquad\qquad\qquad\qquad+ \frac{{n + 1}}{{2n + 1}}p{{\left| 0 \right\rangle }_i}\left\langle 0 \right|, \\&{\sum\limits_{j = 0}^3 {{K_{ij}}\left( {{{\left| 0 \right\rangle }_i}\left\langle 1 \right|} \right)K_{ij}^ + = \sqrt {1 - p} {{\left| 0 \right\rangle }_i}\left\langle 1 \right|,} }\\&{\sum\limits_{j = 0}^3 {{K_{ij}}\left( {{{\left| 1 \right\rangle }_i}\left\langle 0 \right|} \right)K_{ij}^ + = \sqrt {1 - p} {{\left| 1 \right\rangle }_i}\left\langle 0 \right|.} }\end{split}$

这里, 参量$p(t) = 1 - {{\rm e}^{\tfrac{{\varGamma (2 n + 1)}}{2}t}}$与自发衰减系数相关. 对于三体量子态$\rho (0) = \displaystyle\sum\limits_{kls, mnr} {a_{mnr, kls}}({{\left| k \right\rangle }_1} \left\langle m \right| \otimes {{\left| l \right\rangle }_2}\left\langle n \right| \otimes {{\left| s \right\rangle }_3}\left\langle r \right|)$, 任意t时刻的量子态为

$\begin{split}\rho (t) =\;& \sum\limits_{kls,mnr} {{a_{mnr,kls}}\sum\limits_j {K{}_j({{\left| k \right\rangle }_1}\left\langle m \right|)K_j^ + } } \\&\times\sum\limits_{j'} {K{}_{j'}({{\left| l \right\rangle }_2}\left\langle n \right|)K_{j'}^ + } \\&\times\sum\limits_{j''} {K{}_{j''}({{\left| s \right\rangle }_3}\left\langle r \right|)K_{j''}^ + } \end{split}$

根据量子相干度计算公式(5)式, 选取类Werner态为初态, 得到三体量子相干度和任意两组分量子相干度为:

$\begin{split}&{C_{123}} = 2(\left| {a{b^ * }} \right| + \left| {a{c^ * }} \right| + \left| {b{c^ * }} \right|)(1 - p),\\&{C_{12}} = 2\left| {b{c^ * }} \right|(1 - p), {{C_{23}} = } 2\left| {a{b^ * }} \right|(1 - p), \\&{{C_{13}} = 2\left| {a{c^ * }} \right|(1 - p)} .\end{split}$

如果系统初态选择类GHZ态, 那么在环境噪声影响下, 任意t时刻的量子相干度为$C_{123}^{{\rm{GHZ}}} = 2{(1 - p)^{ \textstyle\frac{3}{2}}}\sqrt {q(1 - q)}$, 其他两体量子相干仍然为零. 当自旋相互作用为零时, 初态处于类Werner态的系统, 其量子相干空间分布仍然满足(6)式所描述的加和性. 而且, 量子相干度随着时间呈现指数衰减规律. 噪声参数p越大, 量子相干数值衰减得越快.
然而, 当自旋之间存在非对称性相互作用时, 利用量子主方程的数值计算, 发现三体量子相干度与任意两组分量子相干之和存在差异. 当初态为 $\left| {001} \right\rangle $时, 图1描述了三体量子相干度${C_{123}}$、近邻两组分量子相干${C_{12}}$和${C_{23}}$、以及次近邻两体量子相干${C_{13}}$的含时演化. 如图1所示, 当时间较短时, 状态演化处于初始阶段, 一些能级组分态$\left| {{\psi _{6, 7}}} \right\rangle $的混合, 使得类Werner态组分在量子态中占有一定比重, 所以量子相干度的数值从零逐渐增大. 其中, 对于相互作用较小的两个自旋, 其量子相干度${C_{13}}$和${C_{23}}$在演化初期增长较快. 但是, 随着演化进一步发展, 环境噪声又会导致量子相干度的衰减, 量子退相干现象明显.
图 1 当初态为$ \left| {001} \right\rangle $时, 量子相干组分演化, 参量选取为$ D = 0.2, $ $ \gamma = 0.2, $ $ \varGamma = 0.5, $ $ n = 0.2$, 黑色实线为$ {C_{123}}$, 红色虚线和绿色点划线分别为$ {C_{12}}, {C_{13}}$, 蓝色点线为$ {C_{23}}$
Figure1. The dynamics of all fractions of quantum coherence for the initial state $ \left| {001} \right\rangle $. The parameters are chosen to be$ D = 0.2, $ $ \gamma = 0.2, $ $ \varGamma = 0.5, $ $ n = 0.2$. The black solid line denotes $ {C_{123}}$, the red dashed and green dot-dashed line are $ {C_{12}}, {C_{13}}$ respectively, and the blue dotted line represents $ {C_{23}}$.

为了进一步研究多体量子相干与两体量子相干的关系, 定义一种量子相干组分差值${C_{\rm{r}}} = {C_{123}} - ({C_{12}} + {C_{23}} + {C_{13}})$. 如图2所示, 分析了非均匀自旋相互作用参量D和$\gamma $对量子相干性质的影响. 图2(a)的曲线变化表明, 随着自旋与轨道耦合作用D的增加, 量子相干组分差值${C_{\rm{r}}}$也会增大, 但是环境噪声又会大大抑制量子相干, 导致量子退相干. 同样, 图2(b)的曲线变化表明, 各项异性参量$\gamma $的增加也会引起量子相干组分差值${C_{\rm{r}}}$的增加. 这个结论说明, 非均匀相互作用更有利于三体量子相干度的增加.
图 2 当初态为$ \left| {001} \right\rangle $时, 量子相干组分差值的演化, 参量选取为$ \varGamma = 0.5, $ $ n = 0.2$　(a) 当$ \gamma = 0.2$时, 黑色实线对应参数$ D = 0.3, $ 红色虚线对应$ D = 0.5;$ (b)当$ D = 0.2$时, 黑色实线对应参数$ \gamma = 0.3, $ 红色虚线对应$ \gamma = 0.7$
Figure2. The dynamics of the difference of quantum coherence for the initial state $ \left| {001} \right\rangle $. The parameters are chosen to be $ \varGamma = 0.5, $ $ n = 0.2$: (a) When $ \gamma = 0.2$, the black solid line denotes $ D = 0.3$ the red dashed line is $ D = 0.5;$ (b) When $ D = 0.2$, the black solid line denotes$ \gamma = 0.3, $ the red dashed line is $ \gamma = 0.7$.

当初态选择W态$\left(a = b = c = {1}/{{\sqrt 3 }}\right)$时, 量子相干的各个组分随时间演化行为如图3所示. 所有量子相干组分都呈现出单调衰减现象. 这与(12)式的解析结果相对应.
图 3 当初态为W态时, 量子相干组分演化, 参量$ D = 0.2, $ $ \gamma = 0.2, $ $ \varGamma = 0.5, $ $ n = 0.2$, 　黑色实线为$ {C_{123}}$, 红色虚线和绿色点划线分别为$ {C_{12}}, {C_{13}}$, 蓝色点线为$ {C_{23}}$
Figure3. The dynamics of all fractions of quantum coherence for the initial W state. The parameters are chosen to be $ D = 0.2, $ $ \gamma = 0.2, $ $ \varGamma = 0.5, $ $ n = 0.2$. The black solid line denotes $ {C_{123}}$, the red dashed and green dot-dashed line are $ {C_{12}}, {C_{13}}$ respectively, and the blue dotted line represents $ {C_{23}}$.

但是, 由于自旋相互作用的非对称性, 所以量子相干组分之间存在差异, 相应差值${C_{\rm{r}}}$也会不为零. 图4(a)和图4(b)分别表现出, 自旋与轨道耦合作用D和各项异性参量$\gamma $的增加都会引起三体量子相干度的增大, 使得相应量子相干组分差值在演化初期呈现增长现象. 但是, 由于环境噪声的影响, 无论是三体量子相干度, 还是任意两体量子相干度, 都在较长时间里发生衰减, 直至完全消失.
图 4 当初态为W态时, 量子相干组分差值的演化, 参量选取为$ \varGamma = 0.5, $ $ n = 0.2$　(a) 当$ \gamma = 0.2$时, 黑色实线对应参数$ D = 0.3, $ 红色虚线对应$ D = 0.5;$ (b)当$ D = 0.2$时, 黑色实线对应参数$ \gamma = 0.3, $ 红色虚线对应$ \gamma = 0.7$
Figure4. The dynamics of the difference of quantum coherence for the initial W state. The parameters are chosen to be $ \varGamma = 0.5, $ $ n = 0.2$; (a) When $ \gamma = 0.2$, the black solid line denotes $ D = 0.3$ the red dashed line is $ D = 0.5;$ (b) When $ D = 0.2$, the black solid line denotes$ \gamma = 0.3, $ the red dashed line is $ \gamma = 0.7$.

本文以非对称三体自旋系统为研究对象, 重点分析了量子相干的空间分布特性, 及其在环境噪声背景中的量子退相干现象. 发现类Werner态的量子相干度组分之间存在一种特殊加和性, 即三体量子相干度数值等于所有两体量子相干度之和. 在忽略自旋相互作用的极限条件下, 尽管环境噪声会导致系统所有量子相干组分发生衰减, 但是这种量子加和性仍然存在. 对于可分离初态, 自旋之间非均匀相互作用更有利于三体量子相干度的增加. 这些研究在多体系统可能的量子相干特性规律方面提供了一种思路. 同时, 本文采用的精确对角化和量子主方程, 可以推广到任意多体量子自旋系统, 为进一步研究多体量子相干特性提供了有效方法.