删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

纳米机械谐振器耦合量子比特非厄米哈密顿量诱导的声子阻塞

本站小编 Free考研考试/2021-12-29
摘要:研究了纳米机械谐振器耦合量子比特系统中的声子阻塞现象. 发现在弱驱动的条件下, 非厄米哈密顿量的所有本征值都等于零时, 出现了强声子反聚束现象. 于是提出了零本征值方法来诱导声子阻塞. 本文详细分析了此方法对声子阻塞的影响, 给出了最佳条件并解释了背后的物理机理. 不同于传统的声子阻塞和非传统的声子阻塞, 不管是强非线性还是弱非线性, 都可以实现声子的反聚束效应.
关键词: 声子阻塞/
纳米机械谐振器/
量子比特

English Abstract


全文HTML

--> --> -->
1.引　言
超导量子电路作为人工设计的量子力学系统, 显示了宏观的量子相干性, 可以用来研究量子光学和原子物理的现象[1]. 在超导量子电路中, 超导传输线谐振器代替传统的光学腔, 超导量子比特扮演自然原子的角色. 超导量子比特作为宏观的人工原子, 其他超导量子器件, 例如超导传输线、机械纳米振子等的强耦合, 可以实现一种崭新的全固态的腔量子电动力学(cavity quantum electrodynamics, C-QED)结构, 而且这些人工的原子和光场, 具有以前微腔QED没有的可控性和可集成化的优点. 包含超导量子比特和纳米机械振子的混合量子电路是另一种新兴的混合量子系统[2]. 这种混合量子电路的有趣的理论研究和有实质性进展的实验结果, 证明了该系统可以作为未来量子器件, 并拥有广泛的应用[2-4].
2010年, 加州大学O'Connell等[5]发现机械谐振器在超低温(超导)条件下的振动可以激发量子效应, 产生单个量子(声子). 随后, 清华大学Liu等[6]根据这个实验, 研究了纳米机械谐振器(nanomechanical resonator, NAMR)中的声子阻塞效应. 声子阻塞, 顾名思义, 类似于光子阻塞[7], 是一种纯量子现象. 这种现象是在非线性机械振荡器[8]中, 第一个声子被外部驱动场激发, 第二个声子被第一个声子阻挡导致不能传输. 光子阻塞效应的一项重要应用是可以用于实现理想的单光子源[9,10], 而声子阻塞效应同样可以作为实现单声子源的机制之一[11,12]. 类比于产生光子阻塞的两种机制, 声子阻塞同样地对应两种发生机制. 声子阻塞的实现有两种主要的物理机理: 基于失谐机制的声子反聚束需要较大的非线性来改变系统的能级结构, 这被称为传统声子阻塞(conventional phonon blockade, CPNB)[13-16]. 除了第一种物理机理外, 还可以获得强声子反聚束, 即利用双声子激发的不同路径产生的破坏性干涉, 即使在非线性远小于腔模的衰减率(弱非线性)的条件下也能实现, 称为非传统声子阻塞(unconventional phonon blockade, UPNB)[16-19].
自发现声子阻塞效应起, 国内外人员对这种物理现象进行了大量的研究. 例如双光力系统中, 光子和声子不同耦合路径的干涉产生的光谱差别[20]、光学模式与机械模式之间产生的二阶非线性相互作用导致的强声子反聚束现象[21]、机械谐振器与两能级系统耦合的光力系统的双色二阶边带的生成[22]. 多模光力系统[23]、有限温度下耦合非线性微米/纳米机电系统的谐振器的声子反聚束效应并与热噪声的关系[24]、机械振动位移的平方耦合的声子阻塞现象研究[25]、光声阻塞现象[26]、俘获多个离子的库仑晶体中声子阻塞研究[27]等.
最近, 关于非厄米哈密顿量的系统特性研究得到大量关注[28-44]. 其中包含了量子纠缠[28]、兰姆位移[31]、光子统计特性[44]等. 本文研究NAMR-量子比特系统中的声子统计性质, 其中的系统是由一个纳米机械谐振器耦合一个超导量子比特(以下简称量子比特)组成. 这里利用了非厄米哈密顿量的零本征值方法, 并假设$\kappa < 0$, 研究了NAMR-量子比特系统中的声子统计特性. 同时给出探测声子阻塞或声子统计特性的方法.
2.NAMR-量子比特系统的物理模型
图1研究的是由一个纳米机械谐振器(右边)耦合一个量子比特(左边)组成的物理系统. 纳米机械谐振器受到一个强度为$\varepsilon $的外加场驱动, 量子比特受到一个强度为$\varOmega $的外加场驱动. 整个系统的哈密顿量为[16]
图 1 NAMR-量子比特系统的物理模型图
Figure1. Schematic diagram of the nanomechanical resonator coupled with a qubit.

$\begin{split}{H_1} =\;& {\omega _0}{\sigma _ + }{\sigma _ - } + {\omega _{\rm{m}}}{a^ + }a + g\left( {{\sigma _ + }a + {a^ + }{\sigma _ - }} \right) \\ &+ \left( {\varOmega {{\rm{e}}^{ - {\rm{i}}\phi }}{{\rm{e}}^{ - {\rm{i}}{\omega _{\rm{q}}}t}}{\sigma _ + } + \varepsilon {a^ + }{{\rm{e}}^{ - {\rm{i}}{\omega _{\rm{a}}}t}} + {\rm{H}}.{\rm{c}}.} \right), \end{split}$
其中$\hbar = 1$, ${\sigma _ + }$和${\sigma _ - }$是量子比特的上升和下降算符, 其跃迁频率是${\omega _0}$; a和${a^ + }$表示NAMR中声子的湮灭和产生算符, 其频率是${\omega _{\rm{m}}}$; g是量子比特与NAMR的耦合强度; 量子比特的外加驱动场振幅是$\varOmega $, 频率是${\omega _{\rm{q}}}$; NAMR的外加驱动场振幅是$\varepsilon $, 频率是${\omega _{\rm{a}}}$; $\phi $是两个外加驱动场的相位差.
假设两个外加驱动场的频率相同且量子比特与NAMR的频率相同, 即${\omega _{\rm{q}}} = {\omega _{\rm{a}}} = {\omega _{\rm{d}}}$和${\omega _0} = {\omega _{\rm{m}}}$, 在以入射光频率旋转的框架中, 得出哈密顿量:
$\begin{split}{H_2} =\; & \varDelta {\sigma _ + }{\sigma _ - } + \varDelta {a^ + }a + g\left( {{\sigma _ + }a + {a^ + }{\sigma _ - }} \right) \\ &+ \left( {\varOmega {{\rm{e}}^{ - {\rm{i}}\phi }}{\sigma _ + } + \varepsilon {a^ + } + {\rm{H}}{\rm{.c}}{\rm{.}}} \right), \end{split}$
这里失谐量$\varDelta \equiv {\omega _0} - {\omega _{\rm d}} = {\omega _{\rm m}} - {\omega _{\rm d}}$.
3.非厄米哈密顿量作用下的声子阻塞效应的零本征值
2
3.1.声子阻塞的条件分析
-->

3.1.声子阻塞的条件分析

考虑机械模衰减率$\kappa $和量子比特自发辐射速率$\gamma $. 所以, 整个系统可以用非厄米哈密顿量描述为
${H_3} = {H_2} - {\rm{i}}\frac{\kappa }{2}{a^ + }a - {\rm{i}}\frac{\gamma }{2}{\sigma _ + }{\sigma _ - }.$
考虑系统总激发数为$n + 1$, 在基矢$\left| {n, e} \right\rangle $和$\left| {n + 1, g} \right\rangle $下, 可得到非厄米哈密顿量的矩阵形式:
${\widetilde H^{\left( {n + 1} \right)}} = \left[\!\!\!{\begin{array}{*{20}{c}} {n\!\left( {\varDelta - {\rm{i}}\dfrac{\kappa }{2}} \right) \!+\! \varDelta \!-\! {\rm{i}}\dfrac{\gamma }{2}}&{g\sqrt {n + 1} } \\ {g\sqrt {n + 1} }&{\left( {n + 1} \right)\!\left( {\varDelta - {\rm i}\dfrac{\kappa }{2}} \right)} \end{array}}\!\!\!\!\right].$
考虑弱驱动的条件, 即系统参数$\varOmega , \varepsilon \ll \left| \kappa \right|$.
分析(4)式中的矩阵, 可得知其两个特征值:
$\begin{split}\omega _ \pm ^{\left( {n + 1} \right)} = & \frac{1}{4}\bigg[- {\rm{i}}\left( {\kappa + 2n\kappa + \gamma } \right) + 4\varDelta + 4n\varDelta \bigg. \\ &\pm\left.\sqrt {16{g^2}\left( {1 + n} \right) - {{\left( {\kappa - \gamma } \right)}^2}}\right].\end{split}$
在弱驱动条件下, 考虑系统激发数为1, 即$n = 0$. 而且要使得(5)式特征值取到实数, 则满足条件
$\begin{split}& \kappa = - \gamma, ~~ \left| g \right| \geqslant \frac{\gamma }{2}.\end{split} $
此时, 两个特征值可以化简为
$\omega _ \pm = \varDelta \pm \frac{1}{4}\sqrt {16{g^2} - 4{\gamma ^2}} .$
当$\left| g \right| = \dfrac{\gamma }{2}$时, 两个特征值相等, 即${\omega _ \pm } = \varDelta $. 所有特征值取0时需满足
$\begin{split}& \varDelta = 0, \\& \kappa = - \gamma, \\& \left| g \right| = \frac{\gamma }{2}.\end{split} $
当$g \gg \dfrac{\gamma }{2}$时, 有$\sqrt {16{g^2} - 4{\gamma ^2}} \simeq 4g$, 此时特征值$\omega _ \pm \simeq \varDelta \pm g$.
当$\varDelta = \mp g$时, 这一对特征值相等并且取得最小值0. 此时可以知道特征值为零的条件为:
$\begin{split}& \varDelta = \left| g \right|, \\& \kappa = - \gamma, \\& \left| g \right| \gg \frac{\gamma }{2}.\end{split} $
由图2中的系统能级可知, 系统参数条件满足(8)式或(9)式时, 可出现声子阻塞现象.
图 2 外加驱动场共振的情况下声子阻塞现象的能级图 (a)满足(8)式时的能级图, 其中$\omega _ \pm ^{\left( 1 \right)} = {\omega _{\rm{m}}}$和$\omega _ \pm ^{\left( 2 \right)} = 2{\omega _{\rm{m}}} \pm $ $2g + {\rm{i}}\dfrac{\gamma }{2}$; (b)满足(9)式时的能级图, 其中$\omega _ \pm ^{\left( 1 \right)} = {\omega _{\rm{m}}} \pm g$和$\omega _ \pm ^{\left( 2 \right)} = 2{\omega _{\rm{m}}} \pm \sqrt 2 g + {\rm{i}}\dfrac{\gamma }{2}$
Figure2. Schematic energy-level diagram explaining the occurrence of the phonon blockade by the driving field satisfying the resonance condition: (a) Energy-level diagram when Eqs. (8) are satisfied, where $\omega _ \pm ^{\left( 1 \right)} = {\omega _{\rm{m}}}$ and $\omega _ \pm ^{\left( 2 \right)} = 2{\omega _{\rm{m}}} \pm$ $ 2g + {\rm{i}}\dfrac{\gamma }{2}$; (b) energy-level diagram when Eqs. (9) are satisfied, where $\omega _ \pm ^{\left( 1 \right)} = {\omega _{\rm{m}}} \pm g$ and $\omega _ \pm ^{\left( 2 \right)} = 2{\omega _{\rm{m}}} \pm \sqrt 2 g + {\rm{i}}\dfrac{\gamma }{2}$.

在我们的方案中, $\gamma $设为正数, $\kappa $设为负数. 事实上, 在光学腔中, 腔模衰减率${\kappa _{\rm{C}}} < 0$已被广泛地研究[45,46], 而且已经有实验实现了这个参数[47-51]. 本文基于此, 做出了相应的假设, 即声子场衰减系数$\kappa < 0$. 实际上已有实验通过相干驱动声子技术[5,52-54]、引入声子激光[55-57]和直接驱动机械振子[58,59]等方式来提供机械增益. 这里简要讨论直接驱动纳米机械谐振器来产生力学增益的可能性[60]. 外场驱动在可分辨边带条件下, 对应于纳米机械谐振器中声子冷却的斯托克斯散射过程被极大抑制, 给纳米机械谐振器提供能量的反斯托克斯过程得到了增强并占主要优势.
2
3.2.数值模拟方法
-->

3.2.数值模拟方法

通过计算系统的主方程来研究二阶关联函数. 系统的密度矩阵为:
$\frac{{\partial \rho }}{{\partial t}} = - {\rm{i}}\left[ {H, \rho } \right] + \kappa \ell \left( a \right)\rho + \gamma \ell \left( \sigma \right)\rho , $
这里的超算符$\ell \left( {\hat o} \right)\hat \rho = \hat o\hat \rho {\hat o^ + } - \dfrac{1}{2}{\hat o^ + }\hat o\hat \rho - \dfrac{1}{2}\hat \rho {\hat o^ + }\hat o$.
只需关注稳态时的情况, 即$\partial \rho /\partial t = 0$. 通过二阶关联函数${g^{\left( 2 \right)}}\left( 0 \right)$来研究声子的统计特性.
${g^{\left( 2 \right)}}\left( 0 \right) = \frac{{\left\langle {{a^ + }{a^ + }aa} \right\rangle }}{{{{\left\langle {{a^ + }a} \right\rangle }^2}}}.$
通过计算主方程的稳态解$\left( {\dfrac{{\partial \rho }}{{\partial t}} = 0} \right)$得出${g^{\left( 2 \right)}}\left( 0 \right)$. 二阶关联函数${g^{\left( 2 \right)}}\left( 0 \right) < 1$对应于亚泊松分布的统计特性. 通过计算关联函数, 可以知道声子阻塞效应是否会发生.
2
3.3.结果分析
-->

3.3.结果分析

下面通过计算机数值模拟分析研究声子阻塞效应, 比较和分析(8)式和(9)式. 系统中腔模的希尔伯特空间被截取到5维, 以耗散率$\gamma $作为基本单位来设定其他参数. 图3给出了不同条件下($g/\gamma = 0.5$, 2.0, 4.0) ${g^{\left( 2 \right)}}\left( 0 \right)$与$\varDelta /\gamma $的关系图, 可以得出出现声子阻塞的特点.
图 3 在不同的$g/\gamma $关系下, 二阶关联函数${g^{\left( 2 \right)}}\left( 0 \right)$与$\varDelta /\gamma $的关系图　(a) $\kappa /\gamma = - 1$, $\varepsilon /\gamma = 0.1$, $\varOmega = 0$, $\varepsilon \ne 0$; (b) $\kappa /\gamma = - 1$, $\varOmega /\gamma = 0.1$, $\varepsilon = 0$, $\phi = 0$
Figure3. Logarithmic plot (of base e) of the zero-delay-time second-order correlation functions ${g^{\left( 2 \right)}}\left( 0 \right)$ as a function of $\varDelta /\gamma $ for different $g/\gamma $ with $\kappa /\gamma = - 1$: (a) $\varepsilon /\gamma = 0.1$, $\varOmega = 0$, $\varepsilon \ne 0$; (b) $\varOmega /\gamma = 0.1$, $\varepsilon = 0$, $\phi = 0$.

1) $g/\gamma = 0.5$, 当出现${g^{\left( 2 \right)}}\left( 0 \right)$的波谷, 即${g^{\left( 2 \right)}}\left( 0 \right) \ll 1$时, 对应于强声子反聚束现象. 此时$\varDelta /\gamma = 0$, 与(8)式相对应.
2) $g/\gamma = 2$, 当出现${g^{\left( 2 \right)}}\left( 0 \right) \ll 1$时, 此时$\varDelta /\gamma = $$\pm 2$, 对应于(9)式. 也就是说, 在条件$g/\gamma = 0.5$下比在条件$g/\gamma = 2$下出现的反聚束现象更强. 这是因为这里只取了$\left| g \right| \gg \dfrac{\gamma }{2}$的一个近似值. 只要随着g的增大, 近似效果越明显, 反聚束效应就会越强, 例如图3中的$g/\gamma = 4$.
2
3.4.与CPNB和UPNB比较
-->

3.4.与CPNB和UPNB比较

现在讨论在非厄米特哈密顿量作用下产生声子阻塞效应在零本征值下的物理机理.
图4描述了基于不同路径干涉的情况下的声子阻塞现象的能级情况. 绿色箭头描述的是第一条路径$\left| {1g} \right\rangle \to \left| {2g} \right\rangle $, 振幅为$\sqrt 2 \varepsilon $; 蓝色箭头是第二条路径$\left| {1g} \right\rangle \to \left| {0e} \right\rangle \to \left| {1e} \right\rangle \to \left| {2g} \right\rangle $, 振幅$\dfrac{{g \times \varepsilon \times \sqrt 2 g}}{{\dfrac{\kappa }{2} \times \left( {\dfrac{\kappa }{2} + \dfrac{\gamma }{2}} \right)}}$, 其中$\dfrac{\kappa }{2}$和$\left( {\dfrac{\kappa }{2} + \dfrac{\gamma }{2}} \right)$分别是$\left| {0e} \right\rangle $衰减到$\left| {0g} \right\rangle $和$\left| {1e} \right\rangle $衰减到$\left| {1g} \right\rangle $的衰减率.
图 4 UPNB能级图
Figure4. Schematic energy-level diagram of unconventional phonon blockade.

物理原理上分析, 能级分裂决定了声子阻塞效应的类型是UPNB还是CPNB. 对于(8)式或者(9)式的情况, 驱动场频率${\omega _{\rm{d}}}$如果等于$\omega _ \pm ^{\left( 1 \right)}$, 那就满足了单激发数共振条件, 同时产生单声子的可能性就大大增加. 在这里, 由于失谐, 产生双声子的可能性并没有像单声子那样大.
在少声子数近似下, 系统的态矢量可写作[16]:
$\begin{split}\left| \psi \right\rangle =\;& {C_{0g}}\left| {0g} \right\rangle + {C_{0e}}\left| {0e} \right\rangle + {C_{1g}}\left| {1g} \right\rangle \\ &+ {C_{1e}}\left| {1e} \right\rangle + {C_{2g}}\left| {2g} \right\rangle + {C_{2e}}\left| {2e} \right\rangle .\end{split}$
在弱驱动的条件下, 二阶关联函数
${g^{\left( 2 \right)}}\left( 0 \right) = \frac{{2{{\left| {{C_{2g}}} \right|}^2} + 2{{\left| {{C_{2e}}} \right|}^2}}}{{{{\left( {{{\left| {{C_{1e}}} \right|}^2} + {{\left| {{C_{1g}}} \right|}^2} + 2{{\left| {{C_{2g}}} \right|}^2} + 2{{\left| {{C_{2e}}} \right|}^2}} \right)}^2}}}$
可以写成${g^{\left( 2 \right)}}\left( 0 \right) \approx 2{p_2}/p_1^2$, 这里${p_2} = {\left| {{C_{2g}}} \right|^2} + {\left| {{C_{2e}}} \right|^2}$是产生双声子的几率, ${p_1} = {\left| {{C_{1e}}} \right|^2} + {\left| {{C_{1g}}} \right|^2}$是产生单声子的几率. 由于是在弱驱动的条件下, 所以${p_1} \gg {p_2}$. 在单激发数下的零本征值的条件, ${p_1}$可以急剧增大, 导致的结果是${g^{\left( 2 \right)}}\left( 0 \right) \approx 0$. 这样, 就可以出现声子阻塞现象.
图5比较了CPNB, UPNB和本文方案的区别. 不管是在强非线性还是弱非线性的限制下, 本文方案都可以实现声子的反聚束现象. 实现该方案的物理机理和CPNB相同, 都是由于失谐所造成的声子的反聚束. 实现CPNB的最佳条件如下:
图 5 本文方案与CPNB和UPNB的区别
Figure5. Different parameter range for the CPNB, the UPNB, and the present scheme.

$\begin{split}& \varDelta = \left| g \right|, ~~ \kappa > 0, ~~ \gamma > 0, ~~ \left| g \right| \gg \gamma .\end{split} $
这些最佳条件同时也分析了$\kappa $与$\gamma $之间的定量关系. 这是因为仅仅考虑了(2)式中的特征频率, 而没有引入$\kappa $与$\gamma $. 这就导致了在$\varDelta = \left| g \right|$的条件下, 不能得到$\left| g \right| ={\gamma }/{2}$这个条件(见图6), 而只能注意到$\varDelta = \left| g \right| = 0$. 本文方案最显著的特点就是加入了$\kappa /\gamma = - 1$这个条件. 事实上, 如果(9)式中的条件变为$\kappa > 0$和$\gamma > 0$, 那么就变成了(13)式, 也就是说变为了CPNB.
图 6 二阶关联函数${g^{\left( 2 \right)}}\left( 0 \right)$ 和平均声子数$N = \left\langle {{a^ + }a} \right\rangle = \left[ {{\rm{Tr}}\left( {{\rho _s}{a^ + }a} \right)} \right]$的对数图像　(a), (b) $\varepsilon /\gamma = 0.1$, $\varOmega = 0$; (c), (d) $\varOmega /\gamma = 0.1$, $\varepsilon = 0$, $\phi = 0$. 其他参数为$\kappa /\gamma = - 1$.
Figure6. The logarithmic plot (of base e) of the zero-delay-time second-order correlation functions ${g^{\left( 2 \right)}}\left( 0 \right)$ and average photon numberN: (a) and (b) $\varepsilon /\gamma = 0.1$, $\varOmega = 0$; (c) and (d) $\varOmega /\gamma = 0.1$, $\varepsilon = 0$, $\phi = 0$. The shared parameters: $\kappa /\gamma = - 1$.

图7和图8比较了本文方案和CPNB, UPNB的不同. 结果表明, 在设定相应的参数下, 本文方案得出的声子阻塞效果比其他两种方案要好. 这是因为非厄米哈密顿量的作用方式更贴合主方程. 这就导致了在其他条件相同时, 本文方案会比直接使用(2)式的更精确. CPNB的触发是由于声子场频率与驱动场频率不同步(失谐)而导致的, 而UPNB的触发是通过两个不同路径的量子干涉产生的. 因此, 当作用方式是量子干涉时, 声子的反聚束会对干涉中的相位差极其敏感, 所以干涉路径的改变会使反聚束的区域发生剧烈变化. 然而, 在本文方案中, 驱动场的改变导致量子干涉路径的改变并不会对反聚束的物理环境产生影响. 所以本方案中的强声子反聚束与UPNB方案的物理机理不同.
图 7 二阶关联函数${g^{\left( 2 \right)}}\left( 0 \right)$(对数)与$\varDelta /\gamma $的关系图, 其中$\varepsilon /\gamma = 0.1$, $\varOmega = 0$ (a) $g/\gamma = 0.5$, $\kappa /\gamma = - 1$(本文方案), $\kappa /\gamma = 1$(UPNB); (b) $g/\gamma = 5$, $\kappa /\gamma = - 1$(本文方案), $\kappa /\gamma = 1$(CPNB)
Figure7. The logarithmic plot (of base e) of the zero-delay-time second-order correlation functions ${g^{\left( 2 \right)}}\left( 0 \right)$ as a function of $\varDelta /\gamma $. The shared parameters: $\varepsilon /\gamma = 0.1$, $\varOmega = 0$. Other parameters: (a) $g/\gamma = 0.5$, $\kappa /\gamma = - 1$(our scheme), $\kappa /\gamma = 1$ (UPNB); (b) $g/\gamma = 5$, $\kappa /\gamma = - 1$(our scheme), $\kappa /\gamma = 1$(CPNB).

图 8 二阶关联函数${g^{\left( 2 \right)}}\left( 0 \right)$(对数)与$\varDelta /\gamma $的关系图, $\varOmega \ne 0$, $\varepsilon = 0$, $\phi = 0$ (a) $\varOmega /\gamma = 0.1$, $g/\gamma = 0.5$, $\kappa /\gamma = - 1$(本文方案), $\kappa /\gamma = 1$(UPNB); (b) $\varOmega /\gamma = 2$, $g/\gamma = 10$, $\kappa /\gamma = - 1$(本文方案), $\kappa /\gamma = 1$(CPNB)
Figure8. The logarithmic plot (of base e) of the zero-delay-time second-order correlation functions ${g^{\left( 2 \right)}}\left( 0 \right)$ as a function of $\varDelta /\gamma $: $\varOmega \ne 0$, $\varepsilon = 0$, $\phi = 0$. Other parameters: (a) $\varOmega /\gamma = 0.1$, $g/\gamma = 0.5$, $\kappa /\gamma = - 1$(our scheme), $\kappa /\gamma = 1$ (UPNB); (b) $\varOmega /\gamma = 2$, $g/\gamma = 10$, $\kappa /\gamma = - 1$(our scheme), $\kappa /\gamma = 1$(CPNB).

4.声子阻塞效应的探测
通过测量光学腔场的输出光子的相关性, 可以间接观测到NAMR中的声子阻塞现象. 对于NAMR-量子比特系统, 我们参照了文献中的测量方法[16]. 在最近的实验中, 一个量子比特耦合一个NAMR以后, NAMR再耦合一个超导腔组成了一个三体系统[61]. 假设这个超导光腔通过辐射压相互作用耦合到NAMR, 就可以通过研究超导光腔中的光子统计特性, 来研究NAMR中的声子统计特性, 从而观察NAMR中的声子阻塞现象. 声子统计特性探测模型见图9.
图 9 声子统计特性探测模型
Figure9. The detecting model of the statistical properties of the phonons.

假设超导腔的驱动频率为${\omega _{\rm b}}$, 在以入射光频率${\omega _{\rm{b}}}$旋转的框架中, 整个系统的总哈密顿量为
${H_{\rm{T}}} = {H_{{\rm{mq}}}} + {H_{{\rm{om}}}}, $
这里
$\begin{split}{H_{{\rm{mq}}}} =\;& \varDelta {\sigma _+}{\sigma _-} + \varDelta {a^ + }a + g\left( {{\sigma _ + }a + {a^ + }{\sigma _ -}} \right) \\ &+ \left( {\varOmega {e^{ - {\rm{i}}\phi }}{\sigma _ + } + \varepsilon {a^ + } + {\rm{H}}{\rm{.c}}.} \right), \\H_{{\rm{om}}} =\;& {\varDelta _{\rm{b}}}{b^ + }b + J{b^ + }b\left( {a + {a^ + }} \right) + \left( {{\varOmega _{\rm{c}}}{b^ + } + {\rm{H}}{\rm{.c}}.} \right),\end{split}$
其中b和${b^ + }$表示超导腔中光子湮灭和产生算符, 超导腔的共振频率为${\omega _{\rm{c}}}$; J是超导腔与NAMR的耦合强度; ${\varOmega _{\rm{c}}}$(假定为实数)是超导腔的外加驱动场振幅; ${\varDelta _{\rm{b}}}$满足共振条件${\varDelta _{\rm{b}}} \equiv {\omega _{\rm{c}}} - {\omega _{\rm{b}}} = {\omega _{\rm{m}}}$.
超导腔的腔模算符可以写成量子波动算符和稳态平均场之和: $b \to b + \alpha $. 这里$\alpha $在频率为${\omega _{\rm{b}}}$驱动光场中的稳态平均场. 在强驱动的条件下, $\left| \alpha \right| \gg 1$. 同样地, 可以将${H_{{\rm{om}}}}$再在旋波近似下, 哈密顿量变为$H_{{\rm{om}}}^\prime = {\omega _{\rm{m}}}{b^ + }b + G{a^ + }b + {G^*}a{b^ + }, $其中光力系统有效的耦合强度$G = g\alpha $, 同样假定$\left| G \right| \ll {\omega _{\rm{m}}}$, 在相互作用哈密顿量的演化下, 声子场和腔场之间形成纠缠态, 探测腔场光子统计特性可以实现对声子阻塞或声子统计特性的探测[61].
5.结　论
研究了在非厄米哈密顿量作用下的NAMR-量子比特系统的零本征值条件下的声子统计特性. 充分利用了零本征值方法, 通过分析和计算, 得到了强声子反聚束的最佳条件, 同时引入二阶关联函数讨论了声子阻塞效应. 通过主方程的稳态求解, 分析了各种条件对声子阻塞的影响, 结果发现都能很好地匹配声子阻塞现象. 比较了本方案与CPNB和UPNB的区别, 发现本文方案在强非线性和弱非线性条件下都能实现强声子反聚束现象. 阐明了这种声子阻塞的物理机理是不同于UPNB的失谐原理. 最后, 简要分析了利用外加超导腔去间接地探测NAMR中的声子统计特性的方法. 本文的研究结果将在量子信息处理以及精密测量等研究领域具有重要的意义和应用前景.
相关话题/系统 机械 方案 物理 统计
  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效，10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材，产品类型包含电子书、题库、全套资料以及视频，无论您是考研复习、考证刷题，还是考前冲刺等，不同类型的产品可满足您学习上的不同需求。 ...
    考试优惠券 本站小编 Free壹佰分学习网 2022-09-19
  • 自旋轨道耦合系统中的整数量子霍尔效应
    摘要:研究了Rashba自旋轨道耦合作用下的二维无限长条形样品中的电子输运,计算了样品的霍尔电导和纵向电阻,得到了完整的整数量子霍尔效应.在一定强磁场范围内,由于样品两边缘的限制,能级在大波矢范围快速上升,在小波矢范围形成平坦的朗道能级.强磁场下自旋轨道耦合完全解除自旋简并.位于朗道能级上升和下降区 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于指数加权-核在线序列极限学习机的混沌系统动态重构研究
    摘要:利用指数加权在线核序列极限学习机(exponentialweightedonlinesequentialextremelearningmachinewithkernel,EW-KOSELM)辨识算法,开展了针对混沌动力学系统的动态重构研究.EW-KOSELM算法将核递归最小二乘(kernelr ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于数据库进行乏燃料鉴别的多元统计分析研究
    摘要:近年来,随着国际核军控形势的变化,包含防扩散、防核恐及核安保的多边国际军控合作越来越受到重视.核取证技术作为防扩散、防核恐及核安保的一项核心技术,在对涉核非法活动的威慑、阻止以及响应方面具有重要作用,值得深入研究.目前针对核取证技术的研究较多,主要集中于材料的表征和数据的解读.其中解读作为核取 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 曲波域统计量自适应阈值探地雷达数据去噪技术
    摘要:非线性、非平稳探地雷达数据常掺杂各种复杂噪声源,其对精确提取弱反射波信号、识别绕射波双曲线同相轴特征具有严重影响,忽略噪声影响给探地雷达探测数据全波形偏移成像及后续解译造成较大误差.采用传统阈值函数的曲波变换去噪需要根据数据噪声水平人为确定合理阈值控制系数.对此,本文开展自适应阈值函数的曲波变 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 一种K分布强湍流下的测量设备无关量子密钥分发方案
    摘要:研究了K分布强湍流下自由空间测量设备无关量子密钥分发协议模型,采用阈值后选择方法来减少大气湍流对密钥生成率的影响,对比分析了使用阈值后选择方法前后协议的密钥率和湍流强度之间的关系.仿真结果表明,使用阈值后选择方法可以有效地提高协议的密钥生成率,尤其是在高损耗和强湍流区域,而且其最佳阈值与湍流强 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 线极化Bell-Bloom测磁系统中抽运光对磁场灵敏度的影响
    摘要:利用适用于线极化Bell-Bloom测磁系统的布洛赫方程和含有自旋弛豫的速率方程,以铯原子为研究对象,分析了抽运光对磁场灵敏度的影响,并在实验上分别采用与铯原子D1线和D2线共振的线偏光作为抽运光和探测光,用充有缓冲气体的气室进行了实验.实验结果与理论分析一致,均表明只有在一定的光强范围内,增 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 级联四波混频系统中纠缠增强的量子操控
    摘要:各类系统中的纠缠操控是量子信息科学的重要问题之一.本文研究了热原子蒸气的级联四波混频过程中产生的纠缠增强及纠缠增强的相位敏感特性.研究表明,该级联四波混频过程第二级输出的探针光和共轭光的量子纠缠较第一级明显增强,最大可达5dB以上,且随着强度因子的增大可实现完美纠缠.文中还详细讨论了量子关联类 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于共心球透镜的多尺度广域高分辨率计算成像系统设计
    摘要:针对实时广域高分辨率成像需求,充分利用具有对称结构的多层共心球透镜视场大且各轴外视场成像效果一致性好的特点,设计基于共心球透镜的多尺度广域高分辨率计算成像系统.该系统基于计算成像原理,通过构建像差优化函数获得光学系统设计参数,结合球形分布的次级相机阵列进行全局性优化,提高系统性能的同时有效简化 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 惯性约束聚变装置中靶面光场特性的统计表征方法
    摘要:在激光驱动的惯性约束聚变装置中,常采用多种束匀滑手段对焦斑的时空特性进行调控.光传输链路中涉及的光学元件众多、传输变换复杂,往往导致光传输模型复杂,且在运用衍射光学方法分析焦斑形态和特征时面临大量的数据处理和计算,致使出现计算量大、计算效率低等问题,亟需寻求快速而简便的新方法来描述焦斑的统计特 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 具有多物理特性的X射线脉冲星导航地面验证系统
    摘要:导航地面验证是X射线脉冲星导航研究必不可少的环节.针对导航算法验证需要真实连续的脉冲星信号的需求,同时避免X射线调制及探测难度大、成本高的问题,提出了一种基于可见光源的X射线脉冲星导航地面验证系统.该系统利用太阳系质心处脉冲星信号模型和航天器轨道信息,建立航天器处实时光子到达速率函数,再通过硬 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29

Free考研考试FreeKaoYan.Com
欢迎来到Free考研考试，"为实现人生的Free而奋斗"

© 2020 FreeKaoYan! . All rights reserved.