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

光学腔中一维玻色-哈伯德模型的奇异超固相

本站小编 Free考研考试/2021-12-29
摘要:利用密度矩阵重整化群计算了光学腔中一维无自旋玻色-哈伯德模型的基态. 通过研究超流序、局域密度分布、二阶和三阶关联函数, 发现该系统出现了超越平均场理论的两个奇异超固相. 这两个超固相同时具备对角和非对角长程序, 其中一个展现出包络形式的密度调制振荡, 另一个展现出均匀的密度分布. 另外, 结合光场的超辐射序参量和腔内的平均光子数, 发现奇异超固相与腔光场的涨落存在密切关系. 该工作给出了光学腔内玻色哈伯德模型的超越平均场理论的新物理, 并提供了探索光学腔内光与物质集体物态的完整计算方法.
关键词: 量子相变/
密度矩阵重整化群/
光学腔

English Abstract


全文HTML

--> --> -->
1.引　言
超冷量子气体为多体相互作用系统的量子模拟提供了独特的环境[1]. 将超冷量子气体俘获在光学晶格中可以模拟强关联晶格模型, 其中最典型的模型是玻色-哈伯德模型. 这一模型中存在超流相和莫特绝缘体之间的量子相变[2-4]. 在超冷原子实验中, 原子间具有天然的有限长程相互作用, 它可以通过Feshbach共振技术进行调节[3,4]. 而凝聚态系统中还存在大量由无限长程相互作用驱动的多体现象, 这为基于超冷原子的量子模拟提出了新的挑战. 幸运的是, 光学腔中的光子可以充当原子间沟通的中继站, 进而诱导出有效的无限长程相互作用[5-8]. 在晶格系统中, 该无限长程相互作用会与原子动能和在位相互作用能竞争, 进而衍生出周期是晶格周期两倍的密度关联, 并驱动奇偶格点分布差异的原子密度波[9]. 以上特征使得光学腔内的玻色-哈伯德模型具备更加丰富的基态属性[10-20].
超固相是一种反直觉的物质状态, 它同时具有超流体的无摩擦流与固体的周期性密度结晶[21]. 它也是一个对角序和非对角序(相位相干和空间序)共存的热力学稳态[22-27]. 在平均场理论下, 光学腔中玻色原子的超固相通常出现在光场的超辐射区间, 且晶格的密度分布呈现奇偶振荡[12-18]. 然而, 由于该理论在原理上忽略了物理量的高阶涨落, 故不足以完整描述量子场的密度关联. 特别是当腔光场未发生超辐射时, 腔场涨落可能剧烈增强, 此时平均场理论变得更加不可靠[19,20]. 因此, 为了更加精确地理解光学腔中玻色原子的超固相, 有必要超越传统的平均场理论, 系统计算原子的局域密度分布、二阶和三阶关联函数.
本文考虑装载在光学腔内的一维超冷玻色原子, 通过对腔光场做绝热近似, 得到一种拓展的玻色-哈伯德模型. 这一模型的典型特征是其中包含在位相互作用和无限长程相互作用. 利用密度矩阵重整化群方法[28,29](DMRG)数值计算系统基态的单体超流序、局域密度分布、二阶和三阶关联函数, 探究系统的非对角和对角序. 发现了超越平均场理论的两个奇异超固相, 其中一个展现出包络形式的密度调制振荡, 另一个展现出均匀的密度分布. 另外, 结合光场的超辐射序参量和腔内的平均光子数, 发现奇异超固相与腔光场的涨落存在密切关系.
2.模型和序参量
基于最近实现的光学腔实验平台[10], 本文考虑的系统是装载在准一维(1D)背景光晶格中的超冷玻色气与光学腔的耦合, 如图1所示. 其中光晶格的晶格深度为${V_0}$, 腔的频率为${\omega _{\rm{c}}}$. 所有共振频率为ω的二能级原子不仅与$\hat x$方向的腔场耦合, 还被$\hat z$方向的泵浦光所驱动. 泵浦光驱动频率和振幅分别为${\omega _{\rm{p}}}$和Ω. 为了有效地降低原子发热, 这里只考虑红失谐腔, 即${\varDelta _{\rm{c}}} = {\omega _{\rm{c}}} - {\omega _{\rm{p}}} < 0$. 在大失谐极限时, 激发态可以被绝热消除. 考虑紧束缚近似后的有效哈密顿量为[30]
图 1 左图: 玻色原子沿着腔轴方向$\hat x$被俘获在准一维(1D)背景光学晶格中, 费米气被两束圆偏振的横向(沿着$\hat z$方向) 抽运激光驱动. 右图: 隧穿系数t, 接触型相互作用${U_{\rm{s}}}$和无限长程相互作用${U_{\rm{l}}}$
Figure1. Left: Proposed experimental setup that the bosonic atoms trapped in a quasi-1D optical lattice interact with an optical cavity. Right: Illustration of the competing terms among the hopping t, the contact interaction ${U_{\rm{s}}}$ and the global-range interaction ${U_{\rm{l}}}$.

$\begin{split}H = \;&- t\sum\limits_i {\left( {\hat b_i^{\dagger} {{\hat b}_{i + 1}} + {\rm{H}}.{\rm{c}}{{.}}} \right)} + \frac{{{U_{\rm{s}}}}}{2}\sum\limits_i {{{\hat n}_i}\left( {{{\hat n}_i} - 1} \right)} \\ & - \hbar \left( {{\varDelta _{\rm{c}}} - \delta } \right){{\hat a}^{\dagger} }\hat a \\ &+ \hbar \eta {M_0}\left( {{{\hat a}^{\dagger} } + \hat a} \right)\left( {\sum\limits_{\rm{e}} {{{\hat n}_{\rm{e}}} - \sum\limits_{\rm{o}} {{{\hat n}_{\rm{o}}}} } } \right), \\[-17pt]\end{split} $
式中, $\hat b_i^{\dagger} $和${\hat b_i}$为基态能级的玻色子产生和湮灭算符, 粒子数算符${\hat n_i} = \hat b_i^{\dagger} \hat b_i^{}$. ${\hat a^{\dagger} }$和$\hat a$是量子化腔场的产生和湮灭算符. o和e分别是奇偶格点指标. t为跃迁系数, ${U_{\rm{s}}}$是玻色子排斥相互作用强度. $\delta = {M_1}N$是腔的色散偏移, 其中${M_0} = \displaystyle \int {\rm{d}}xW_i^*\left( x \right) \cos \left( {{k_0}x} \right) $$ {W_i}\left( x \right)$, ${M_1} = \left( {{g^2}/\varDelta } \right)\displaystyle \int {{\rm{d}}xW_i^*\left( x \right){{\cos }^2}\left( {{k_0}x} \right){W_i}\left( x \right)}$, N是原子总数, ${W_i}\left( x \right)$是局域在格点的瓦尼儿函数, ${k_0}$为腔场波矢, g为腔内单光子耦合常数, $\eta = $$ g\varOmega /\varDelta$是有效的光与原子耦合强度.
根据方程(1)可以写出光场消灭算符所满足的量子郎之万方程:
$\begin{split}{\partial _t}\hat a(t) =\;& {\rm{i}}({\varDelta _{\rm{c}}} - \delta )\hat a(t) - {\rm{i}}\eta {M_0}\left( {\sum\limits_{\rm{e}} {{{\hat n}_{\rm{e}}}} - \sum\limits_{\rm{o}} {{{\hat n}_{\rm{o}}}} } \right) \\&- \kappa \hat a(t) - \sqrt {2\kappa } \hat \varGamma (t),\\[-10pt]\end{split}$
方程最后一项为库环境引起的量子噪声. 考虑腔场的耗散κ远大于原子反冲能, 此时可作绝热近似, 即忽略噪声项, 认为光场始终驰豫在稳态[30]. 从方程(2)出发得到光场的稳态形式为[10]
$\hat a = \frac{{\eta {M_0}}}{{{\varDelta _{\rm{c}}} - \delta + {\rm{i}}\kappa }}\left( {\sum\limits_{\rm{e}} {{{\hat n}_{\rm{e}}} - \sum\limits_{\rm{o}} {{{\hat n}_{\rm{o}}}} } } \right).$
将方程(3)引入(1)式的哈密顿量, 可以得到有效哈密顿量为
$\begin{split} H = \;&- t\sum\limits_i {\left( {\hat b_i^{\dagger} {{\hat b}_{i + 1}} + {\rm{H}}.{\rm{c}}{{.}}} \right)} + \frac{{{U_{\rm{s}}}}}{2}\sum\limits_i {{{\hat n}_i}\left( {{{\hat n}_i} - 1} \right)} \\& - \frac{{{U_{\rm{l}}}}}{L}{\left( {\sum\limits_{\rm{e}} {{{\hat n}_{\rm{e}}} - \sum\limits_{\rm{o}} {{{\hat n}_{\rm{o}}}} } } \right)^2},\\[-18pt] \end{split} $
其中, ${U_{\rm{l}}} = - L\hbar {\left| {\eta {M_0}} \right|^2}( {{\varDelta _{\rm{c}}} - \delta } )/\left[ {{{\left( {{\varDelta _{\rm{c}}} - \delta } \right)}^2} + {\kappa ^2}} \right]$[10], L是晶格的尺寸. 哈密顿量(4)中的${U_{\rm{l}}}$项清楚地展示了光学腔中以光子为媒介诱导的原子原子之间无限长程相互作用, 它可以通过改变泵浦激光的大小来独立调整. 至此就得到了一个隧穿、在位相互作用和无限长程相互作用共存且竞争的有趣模型, 如图1所示.
为了探究超固相的属性, 应该考虑非对角序和对角序的实验可观测序参量: 超流序(非对角序)
$G(r) = \frac{1}{L}\sum\limits_i {\big\langle {\hat b_i^{\dagger} {{\hat b}_{i + r}}} \big\rangle } .$
它描述的是相位相干. 傅里叶形式为$n(k) = \dfrac1L\displaystyle\sum\nolimits_r $$ {{{\rm{e}}^{{\rm{i}}kr}}G(r)}$. $n(k = 0)$是凝聚系数, 表示占据超流态的玻色子的密度. 如果$G(r)$具有长程序, 则$n( k = $$ 0 ) > 0$, 此时系统表现出相位相干. 局域密度分布为
${D^1}\left( j \right) = \left\langle {{{\hat n}_j}} \right\rangle .$
它的傅里叶变换为$\nu \left( k \right) = \dfrac1L\displaystyle\sum\nolimits_j {{{\rm{e}}^{{\rm{i}}kj}}\left\langle {{{\hat n}_j}} \right\rangle }$. $\nu \left( k \right)$的尖峰位置所对应的动量值${\nu _{{\rm{peak}}}}$描述的是密度分布的振荡周期. 如果密度分布$\left\langle {{{\hat n}_j}} \right\rangle $呈现的是奇偶晶格振荡, 则${\nu _{{\rm{peak}}}} = {\text{π}}$; 如果密度分布$\left\langle {{{\hat n}_j}} \right\rangle $呈现的是均匀分布, 则${\nu _{{\rm{peak}}}} = 0$; 如果密度分布$\left\langle {{{\hat n}_j}} \right\rangle $呈现的是包络型的非公度振荡分布, 则$\nu \left( k \right)$会出现两个尖峰, 且每个尖峰都满足$0 < {\nu _{{\rm{peak}}}} < {\text{π}}$. 二阶关联函数(对角序)必须去除掉局域密度分布$\left\langle {{{\hat n}_j}} \right\rangle $的贡献, 即
${D^2}\left( {l,j} \right) = \left\langle {{{\hat n}_l}{{\hat n}_j}} \right\rangle - \left\langle {{{\hat n}_l}} \right\rangle \left\langle {{{\hat n}_j}} \right\rangle .$
它的傅里叶变换为$S(k) = \dfrac1L\displaystyle\sum\nolimits_{l, j} {{{\rm{e}}^{{\rm{i}}k\left( {l - j} \right)}}\left\langle {{{\hat n}_l}{{\hat n}_j}} \right\rangle }$. 当$\left\langle {{{\hat n}_l}{{\hat n}_j}} \right\rangle $表现为奇偶振荡, 则$S(k = {\text{π}}) > 0$, 此时系统存在对角序. 三阶关联函数必须减去来自二阶关联${D^2}\left( {l, j} \right)$和局域密度分布$\left\langle {{{\hat n}_j}} \right\rangle $的贡献, 它的定义为[31-35]
$\begin{split} {D^3}\left( {i,l,j} \right) = \;&\left\langle {{{\hat n}_i}{{\hat n}_l}{{\hat n}_j}} \right\rangle - {D^2}\left( {l,j} \right)\left\langle {{{\hat n}_i}} \right\rangle - {D^2}\left( {i,j} \right)\left\langle {{{\hat n}_l}} \right\rangle \\ & - {D^2}\left( {i,l} \right)\left\langle {{{\hat n}_j}} \right\rangle - \left\langle {{{\hat n}_i}} \right\rangle \left\langle {{{\hat n}_l}} \right\rangle \left\langle {{{\hat n}_j}} \right\rangle . \\[-10pt]\end{split} $
目前关于光学腔中玻色子的研究大多只关注原子的量子效应, 而较少探究光子的性质[19,20]. 事实上, 光学腔中的光子是原子之间的耦合媒介, 它在系统的集体动力学中往往发挥着重要作用. 这里重点考察两个与光子相关的物理量, 第一个是基于平均场理论的超辐射序参量:
${\left| {\left\langle {\hat a} \right\rangle } \right|^2} = \frac{{{U_{\rm{l}}}}}{L}{\left| {\left\langle {\sum\limits_{\rm{e}} {{{\hat n}_{\rm{e}}} - \sum\limits_{\rm{o}} {{{\hat n}_{\rm{o}}}} } } \right\rangle } \right|^2}.$
另一个是腔中光子数:
$\left\langle {{{\hat a}^{\dagger} }\hat a} \right\rangle = {U_{\rm{l}}}LS\left( {k = {\rm{\pi }}} \right).$
在平均场理论中, 超辐射序参量和腔中的光子数相等, 即${\left| {\left\langle {\hat a} \right\rangle } \right|^2} = \left\langle {{{\hat a}^{\dagger} }\hat a} \right\rangle $[36]. 在超辐射区间, ${\left| {\left\langle {\hat a} \right\rangle } \right|^2} > 0$, 腔场涨落较小, 平均场理论一般能较好描述系统的动力学性质. 而对于非超辐射区间, 超辐射序参量为零, 即${\left| {\left\langle {\hat a} \right\rangle } \right|^2} = 0$, 此时腔场涨落可能会对系统状态有较大影响[36].
本文用DMRG有限尺寸算法计算系统开边界情况下的多体基态. 选取的格点长度L最大为96, DMRG最大态数为200, 扫描次数为30, 此时的截断误差最大为${10^{ - 8}}$. 令$t = 1$作为能量单位.
3.量子多体相
首先, 将一个较大的无限长程相互作用(例如${U_{\rm{l}}} = 6$)引入到${U_{\rm{s}}} = 10$的玻色-哈伯德模型中. DMRG数值结果显示, 超流序$G(r)$随着格点间隔r的增大呈现幂次衰减(幂率值为0.7)且奇偶振荡, 此时$n(k = 0) > 0$, 如图2(a1)所示. 这说明此时体系存在相位相干, 即存在非对角长程序. 密度分布$\left\langle {{{\hat n}_i}} \right\rangle $呈现奇偶振荡, 此时${\nu _{{\rm{peak}}}} = {\text{π}}$, 如图2(a2)所示. 消除了密度分布$\left\langle {{{\hat n}_i}} \right\rangle $的贡献, 二阶关联${D^2}(l, j)$接近0, 如图2(a3)所示. 二阶关联的贡献全部来自奇偶振荡的密度分布$\left\langle {{{\hat n}_i}} \right\rangle $. 这一参数下系统对应的基态与平均场近似理论的结果一致. 我们又计算了这一参数下的三阶关联函数${D^3}(i, l, j)$. 消除了密度分布$\left\langle {{{\hat n}_i}} \right\rangle $和二阶关联${D^2}(l, j)$的贡献, 三阶关联${D^3}(i, l, j)$接近0, 如图2(a4)所示. 本文将这一传统的超固相称为超固0相.
图 2 (a1), (b1), (c1)超流序$G(r)$和插图$n(k)$; (a2), (b2), (c2)密度分布$\left\langle {{{\hat n}_j}} \right\rangle $; (a3), (b3), (c3)二阶关联${D^2}(l, j)$; (a4), (b4), (c4)三阶关联${D^3}\left( {i, l, j} \right)$. (a1)—(a3) 超固0相, 相互作用参数${U_{\rm{s}}} = 10$和${U_{\rm{l}}} = 6$; (b1)—(b3) 超固1相, 相互作用参数${U_{\rm{s}}} = 10$和${U_{\rm{l}}} = 4.4$; (c1)—(c3) 超固2相, 相互作用参数${U_{\rm{s}}} = 2$和${U_{\rm{l}}} = 4.4$. 所有图中, 共有参数$L = 96$, $\rho = 0.4375$
Figure2. (a1), (b1), (c1) Superfluid order $G(r)$ and $n(k)$; (a2), (b2), (c2) density profile $\left\langle {{{\hat n}_j}} \right\rangle $; (a3), (b3), (c3) two order correlation ${D^2}\left( {l, j} \right)$; (a4), (b4), (c4) three order correlation ${D^3}\left( {i, l, j} \right)$ for $i = j$ of (a1)–(a3) supersolid 0 phase with ${U_{\rm{s}}} = 10$ and ${U_{\rm{l}}} = 6$, (b1)–(b3) supersolid 1 phase with ${U_{\rm{s}}} = 10$ and ${U_{\rm{l}}} = 4.4$, and (c1)–(c3) supersolid 2 phase with ${U_{\rm{s}}} = 2$ and ${U_{\rm{l}}} = 4.4$. In all subfigure, we have $L = 96$ and $\rho = 0.4375$.

当${U_{\rm{s}}} = 10$的玻色-哈伯德模型中引入一个大小适中的无限长程相互作用(以${U_{\rm{l}}} = 4.4$为例). DMRG数值结果显示, 超流序$G(r)$随着格点间隔r的增大幂率衰减(幂率值为0.56), 且此时$n(k = 0) > 0$, 如图2(b1)所示. 密度分布$\left\langle {{{\hat n}_i}} \right\rangle $呈现包络型振荡, 此时$0{{ < }}{\nu _{{\rm{peak}}}} < {\text{π}}$, 如图2(b2)所示. 二阶关联${D^2}(l, j)$呈现奇偶振荡且不随格点间距增大而衰减, 如图2(b3)所示. 二阶关联的贡献全部来自关联效应, 这一数值结果完全超越了平均场理论. 三阶关联函数${D^3}(i, l, j)$呈现奇偶振荡且不随格点间距增大而衰减, 如图2(b4)所示. 这一包络型分布的奇偶格点原子总数相等, 即$\left\langle { \displaystyle\sum\nolimits_{\rm{e}} {{{\hat n}_{\rm{e}}}} - \sum\nolimits_{\rm{o}} {{{\hat n}_{\rm{o}}}} } \right\rangle = 0$. 然而, 无限长程相互作用能$\left\langle {{{\Big(\displaystyle\sum\nolimits_{\rm{e}} {{{\hat n}_{\rm{e}}}} - \sum\nolimits_{\rm{o}} {{{\hat n}_{\rm{o}}}} \Big)}^2}} \right\rangle > 0$. 这两个值之间的差异主要源自于量子涨落. 本文中, 具有$0{{ < }}{\nu _{{\rm{peak}}}} < {\text{π}}$特征的超固相被称为超固1相.
降低在位相互作用大小(以${U_{\rm{s}}} = 2$为例), 但无限长程相互作用不变, 仍为${U_{\rm{l}}}{{ = }}4.4$. $G(r)$随着格点间隔r的增大幂率衰减(幂率值为0.24), 且此时$n(k = 0) > 0$, 如图2(c1)所示. 此时体系存在相位相干, 即存在非对角长程序. 这一参数对应的基态同样为超固相, 但是这一超固相不同于超固0相和超固1相. 密度分布$\left\langle {{{\hat n}_i}} \right\rangle $是均匀的, 此时${\nu _{{\rm{peak}}}} = 0$, 如图2(c2)所示. 二阶关联${D^2}(l, j)$和三阶关联${D^3}(i, l, j)$呈现奇偶振荡且不随格点间距增大而衰减, 如图2(b3)和图2(b4)所示. 本文中, 特征为${\nu _{{\rm{peak}}}} = 0$的超固相被称超固2相.
为了更好地展示超流序$G(r)$的存在, 计算了多个尺寸的数据, 并做了有限尺寸分析, 如图3所示. 在有限尺寸下, 超固1相和超固2相的超流序$G(r)$随着距离呈幂率衰减, 将幂率做有限尺寸分析发现, 在热力学极限下($L \to \infty $)幂率在$\left( {0, 1} \right)$区间, 如图3(a1)和图3(b1)所示. 对$n(k = 0){{/}}L$做有限尺寸分析可以发现, 在热力学极限下($L \to \infty $)超固1相和超固2相的$n(k = 0){{/}}L > 0$, 如图3(a2)和图3(b2)所示. 综上可得, 超固1相和超固2相的超流序是固然存在的, 不因尺寸变化而消失.
图 3 (a1), (b1)不同尺寸下的超流序$G(r)$, 插图为冥率的有限尺寸分析; (a2), (b2) $n(k = 0){{/}}L$ 的有限尺寸分析. (a1), (a2)超固1相, 参数${U_{\rm{s}}} = 10$和${U_{\rm{l}}} = 4.4$; (b1), (b2) 超固2相, 参数${U_{\rm{s}}} = 2$和${U_{\rm{l}}} = 4.4$. 所有图中, 共有参数$\rho = 0.4375$
Figure3. (a1), (b1) Superfluid order $G(r)$ of several lattice length. The inset is the finite-size scaling of the decay rate. (a2), (b2) Finite-size scaling of $n(k = 0){{/}}L$. (a1), (a2) Supersolid 1 phase with ${U_{\rm{s}}} = 10$ and ${U_{\rm{l}}} = 4.4$. (b1), (b2) Supersolid 2 phase with ${U_{\rm{s}}} = 2$ and ${U_{\rm{l}}} = 4.4$. In all subfigure, we have $\rho = 0.4375$.

超固1相和超固2相源自于在位相互作用和无限长程相互作用的共存和竞争. 超固1相在不同尺寸情况下都会出现相同周期的包络型的密度分布, 如图4(a1)—(a3)所示. 且不同尺寸的${\nu _{{\rm{peak}}}}$相同, 通过有限尺寸分析可以发现, 热力学极限下($L \to $$ \infty$)的${\nu _{{\rm{peak}}}}$与有限尺寸的${\nu _{{\rm{peak}}}}$相等, 如图4(a4)所示. 超固2相在不同尺寸时的密度分布$\left\langle {{{\hat n}_j}} \right\rangle $都是均匀的, 如图4(b1)—(b3)所示. 不同尺寸包括热力学极限下($L \to \infty $)的${\nu _{{\rm{peak}}}} \equiv 0$, 如图4(b4)所示.
图 4 (a1)—(a3)和(b1)—(b3)密度分布$\left\langle {{{\hat n}_j}} \right\rangle $; (a4), (b4) ${\nu _{{\rm{peak}}}}$的有限尺寸分析. (a1)—(a4) 超固1相, 参数${U_{\rm{s}}} = 10$和${U_{\rm{l}}} = 4.4$; (b) 超固2相, 参数${U_{\rm{s}}} = 2$和${U_{\rm{l}}} = 4.4$. (a1), (b1) $L = 80$; (a2), (b2) $L = 64$; (a3), (b3) $L = 48$. 所有图中, 共有参数$\rho = 0.4375$
Figure4. (a1) –(a3) and (b1) –(b3) the density profile $\left\langle {{{\hat n}_j}} \right\rangle $; (a4), (b4) the finite-size scaling of the ${\nu _{{\rm{peak}}}}$. (a1)–(a4) Supersolid 1 phase with ${U_{\rm{s}}} = 10$ and ${U_{\rm{l}}} = 4.4$; (b1)–(b4) supersolid 2 phase with ${U_{\rm{s}}} = 2$ and ${U_{\rm{l}}} = 4.4$. (a1), (b1) $L = 80$; (a2), (b2) $L = 64$; (a3), (b3) $L = 48$. In all subfigure, we have $\rho = 0.4375$.

本文利用$n(k = 0)$和$S(k = {\text{π}})$的值, 即可辨别超流、莫特绝缘体、电荷密度波和超固相. 但是, 需要使用${\nu _{{\rm{peak}}}}$来区分超固0相、超固1相和超固2相, 如表1所列.
序/相简写超流SF莫特绝缘体MI电荷密度波CDW超固0相SS0超固1相SS1超固2相SS2
$n~(k = 0)$> 000> 0> 0> 0
$S~(k = {\text{π} })$00> 0> 0> 0> 0
${\nu _{{\rm{peak}}}}$00ππ$(0, {\text{π}})$0


表1量子相对应的序参量
Table1.Corresponding orders of the phases.

为了清楚地探究上述超固相之间的相变, 固定在位相互作用为${U_{\rm{s}}} = 10$和填充因子$\rho = 0.4375$, 将无限长程相互作用强度${U_{\rm{l}}}$从0不断增大. 这一过程中$n(k = 0) > 0$. 当$0 < {U_{\rm{l}}} < 3.3$时, $S(k = {\text{π}}) = 0$和${\nu _{{\rm{peak}}}} = 0$, 对应的量子相为超流; 当$3.3 < {U_{\rm{l}}} < $$ 4.6$时, $S(k = {\text{π}}) > 0$和${\nu _{{\rm{peak}}}} \in (0, {\text{π}})$, 对应的量子相为超固1相; 当${U_{\rm{l}}} > 4.6$时, $S(k = {\text{π}}) > 0$和${\nu _{{\rm{peak}}}} = {\text{π}}$对应的量子相为超固0相, 如图5(a)所示. 基态从超流相变为超固1相, 最终又相变为超固0相, 临界点为$U_{\rm{l}}^{\rm{c}} = 3.3$和4.6. 超辐射相变同样发生在无限长程相互作用${U_{\rm{l}}}$较大时, 如图5(b)所示. 超固0相是超辐射驱动的, 在此区域中DMRG的结果显示${\left| {\left\langle {\hat a} \right\rangle } \right|^2} = \left\langle {{{\hat a}^{\dagger} }\hat a} \right\rangle > 0$. 然而, 超固1相不是超辐射驱动, 是腔光场涨落驱动的. 此时${\left| {\left\langle {\hat a} \right\rangle } \right|^2} = 0$且${\left| {\left\langle {\hat a} \right\rangle } \right|^2} \ne \left\langle {{{\hat a}^{\dagger} }\hat a} \right\rangle $, 如图5(b)所示.
图 5 $S(k = {\text{π}})$(黑色实线)和 ${\nu _{{\rm{peak}}}}$(蓝色方块)　(a) 关于${U_{\rm{l}}}$的变化, 固定${U_{\rm{s}}} = 10$; (c) 关于${U_{\rm{s}}}$ 的变化, 固定${U_{\rm{l}}} = 4.4$. 平均光子数${\left| {\left\langle {\hat a} \right\rangle } \right|^2}$和腔中的光子数$\left\langle {{{\hat a}^{\dagger} }\hat a} \right\rangle $关于 (b) 关于${U_{\rm{l}}}$的变化, 固定${U_{\rm{s}}} = 10$; (d) 关于${U_{\rm{s}}}$的变化, 固定${U_{\rm{l}}} = 4.4$. 所有图中, 共有参数$L = 96$和$\rho = 0.4375$. SF表示超流, MI表示莫特绝缘体, CDW表示电荷密度波, SS0表示超固0相, SS1表示超固1相, SS2表示超固2相
Figure5. The $S(k = {\text{π}})$(black solid line) and ${\nu _{{\rm{peak}}}}$(blue square) as function of (a) ${U_{\rm{l}}}$ with ${U_{\rm{s}}} = 10$, and (c) ${U_{\rm{s}}}$ with ${U_{\rm{l}}} = 4.4$. The number of cavity photons $\left\langle {{{\hat a}^{\dagger} }\hat a} \right\rangle $ and mean cavity-field ${\left| {\left\langle {\hat a} \right\rangle } \right|^2}$ as a function of (b) ${U_{\rm{l}}}$ with ${U_{\rm{s}}} = 10$, (d) ${U_{\rm{s}}}$ with ${U_{\rm{l}}} = 4.4$. In all subfigure, we have $L = 96$ and $\rho = 0.4375$. SF denotes superfluid, MI denotes Mott insulator, CDW denotes charge density wave, SS0 denotes supersolid 0, SS1 denotes supersolid 1, SS2 denotes supersolid 2.

当固定无限长程相互作用强度${U_{\rm{l}}} = 4.4$和填充因子$\rho = 0.4375$, 将在位相互作用强度${U_{\rm{s}}}$从0不断增大. 这一过程中$n(k = 0) > 0$和$S(k = {\text{π}}) > 0$. 当$0 < {U_{\rm{s}}} < 5.5$时${\nu _{{\rm{peak}}}} = 0$, 对应的量子相为超固2相; 当$5.5 < {U_{\rm{s}}} < 11.5$时${\nu _{{\rm{peak}}}} \in (0, {\text{π}})$, 对应的量子相为超固1相; 当${U_{\rm{s}}} > 11.5$时${\nu _{{\rm{peak}}}} = {\text{π}}$, 对应的量子相为超固0相, 如图5(c)所示. 我们同样计算了这一过程的超辐射序参量和光子数, 发现超固1相和超固2相这两个相的${\left| {\left\langle {\hat a} \right\rangle } \right|^2} = 0$且${\left| {\left\langle {\hat a} \right\rangle } \right|^2} \ne \left\langle {{{\hat a}^{\dagger} }\hat a} \right\rangle $, 如图5(d)所示. 综上, 可以得出结论, 超固1相和超固2相是腔光场涨落驱动的.
基于上述描述的相变, 勾勒了哈密顿量(3)的相图, 如图6所示. 图6(a)—(c)给出了${U_{\rm{l}}}$-ρ平面的相图, 其中(a) ${U_{\rm{s}}} = 0$, (b) ${U_{\rm{s}}} = 5$和(c) ${U_{\rm{s}}} = 10$. 当无在位相互作用和无限长程相互作用时, 体系基态在任何填充率ρ下都是超流. 当引入无限长程相互作用${U_{\rm{l}}}$后体系变为超固2相, 继续增强${U_{\rm{l}}}$体系变为超固0相, 如图6(a)所示. 当${U_{\rm{s}}} = 5$时, 无限长程相互作用${U_{\rm{l}}}$驱动体系产生超固0相、超固1相和超固2相. 只有$\rho = 0.5$和1时, 体系会出现莫特绝缘体和电荷密度波, 如图6(b)所示. 当${U_{\rm{s}}} = 10$时, 超固2相被压缩而消失, 相图中仍然有超流、莫特绝缘体、电荷密度波、超固0相和超固1相, 如图6(c)所示. 为了更直观地展示在位相互作用和无限长程相互作用的竞争, 给出了${U_{\rm{l}}}$-${U_{\rm{s}}}$平面的相图, 如图6(d)所示. 无限长程相互作用${U_{\rm{l}}}$驱动产生超固相, ${U_{\rm{l}}}$很大时体系为超固0相. 当无限长程相互作用${U_{\rm{l}}} = 4$左右时, 在位相互作用${U_{\rm{s}}}$驱动超固2相变为超固1相.
图 6 (a)—(c) ${U_{\rm{l}}}$-ρ平面的相图　(a) ${U_{\rm{s}}} = 0$; (b) ${U_{\rm{s}}} = 5$; (c) ${U_{\rm{s}}} = 10$. (d) ${U_{\rm{l}}}$-${U_{\rm{s}}}$平面的相图, 参数$\rho = 0.4375$. 所有图中, 共有参数$L = 96$. SF表示超流, MI表示莫特绝缘体, CDW表示电荷密度波, SS0表示超固0相, SS1表示超固1相, SS2表示超固2相
Figure6. (a)–(c) Phase-diagram between ${U_{\rm{l}}}$ and ρ with (a) ${U_{\rm{s}}} = 0$, (b) ${U_{\rm{s}}} = 5$ and (c) ${U_{\rm{s}}} = 10$. (d) Phase-diagram between ${U_{\rm{l}}}$ and ${U_{\rm{s}}}$, with filling $\rho = 0.4375$. In all subfigure, we have $L = 96$. SF denotes superfluid, MI denotes Mott insulator, CDW denotes charge density wave, SS0 denotes supersolid 0, SS1 denotes supersolid 1, SS2 denotes supersolid 2.

4.总　结
上述量子相在目前的实验平台上都是可以观测的, 动量分布可以通过时间飞行探测, 结构因子可以通过布拉格散射探测[37,38], 密度分布可以通过单格点分辨进行探测[39].
本文研究了光学腔中一维无自旋玻色-哈伯德模型. 利用密度矩阵重整化群方法, 计算了基态的单体超流序、局域密度分布、二阶和三阶关联函数. 发现系统出现了存在对角和非对角的长程序并超越平均场理论的两个奇异超固相. 这两个超固相分别展现出包络型的密度调制振荡和均匀的密度分布. 通过计算光子数和腔场平均, 发现奇异超固相主要由腔场涨落诱导. 相互作用能够驱动奇异超固相之间的相变. 本文的工作给出了光学腔内玻色-哈伯德模型的超越平均场理论的新物理, 并提供了充分探索光学腔的完整计算方法.
相关话题/光学 系统 光子 辐射 计算
  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效，10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材，产品类型包含电子书、题库、全套资料以及视频，无论您是考研复习、考证刷题，还是考前冲刺等，不同类型的产品可满足您学习上的不同需求。 ...
    考试优惠券 本站小编 Free壹佰分学习网 2022-09-19
  • X射线荧光CT成像中荧光产额、退激时间、散射、偏振等关键物理问题计算与分析
    摘要:X射线荧光CT(X-rayfluorescencecomputedtomography,XFCT)是一种使用X射线荧光(X-rayfluorescence,XRF)实现功能性成像的新技术,在生物医学成像中表现出较大潜力.但是,X射线穿过生物体的同时还会产生大量康普顿散射光子,对XRF信号的采集 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • Si<sub><i>n</i></sub>团簇/石墨烯(<i>n</i> ≤ 6)结构稳定性和储锂性能的第一性原理计算
    摘要:目前,硅/碳复合材料是锂离子电池最有潜在应用前景的高容量负极材料之一,硅与碳材料的界面状态是影响其电化学性能的重要因素.本文在作为碳材料结构单元的石墨烯表面构建了Sin(n≤6)团簇,采用基于密度泛函理论(DFT)的第一性原理方法研究了Sin团簇/石墨烯(Sin/Gr)的几何构型、结构稳定性和 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 量子密钥分发系统中抗扰动偏振编码模式的实验研究
    摘要:基于移动设备的自由空间量子密钥分发系统能够为实时全覆盖多节点网络提供有效的解决方案,然而该系统目前依旧存在着稳定性不足的问题,偏振编码器的抗扰动性在移动设备应用中非常重要.在扰动情况下保偏光纤的偏振保持特性将受到很大程度的影响,进而使得基于该特性的编码器的稳定性也随之受到影响.为解决扰动情况下 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 双光梳非线性异步光学采样测距中关键参数的数值分析
    摘要:双光梳异步光学采样的绝对测距方法具有量程大、测速快和精度高等特点,在几何量精密测量领域具有广泛的应用前景.特别地,结合异步光学采样和非线性强度互相关的倍频信号时域探测方法,可以有效避免测量过程中载波包络偏移频率对测距精度的影响.本文针对双光梳非线性异步光学采样绝对测距系统,对影响其测距精度的理 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • “强光一号”加速器短<i>γ</i>二极管径向箍缩率计算方法
    摘要:电子束流箍缩是强流电子束二极管工作过程中广泛存在的物理现象.束流径向箍缩率定义为靶面形成的束斑(环)面积随时间的变化率,是判断二极管的束流箍缩情况和工作特性的重要指标,目前对其的研究方法以光学诊断和针对特定二极管的理论估算为主.在现有研究的基础上,以“强光一号”强箍缩短γ二极管为对象,给出了适 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于介电常数近零模式与间隙表面等离激元强耦合的增强非线性光学效应
    摘要:非线性光学效应在光通信、光探测、量子信息等领域发挥着举足轻重的作用,然而天然材料的光学非线性响应通常很弱.本文利用氧化铟锡(ITO)薄膜激发的介电常数近零模式,与金属-介质-金属结构激发的间隙表面等离激元发生强耦合,在近红外波段实现宽带(约1000nm)的增强非线性光学效应,其非线性折射率n2 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 高效、偏振不敏感超导纳米线单光子探测器
    摘要:超导纳米线单光子探测器(SNSPD)因其优异的综合性能被广泛应用于量子通信等众多领域,然而其独特的线性结构会导致SNSPD的探测效率对入射光的偏振态具有依赖性,从而限制了SNSPD在非常规光纤链路或其他非相干光探测环境中的应用.本文基于传统的回形纳米线结构设计制备了一种新型偏振不敏感SNSPD ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 百微米精度的单光子测距
    摘要:本文发展了一种基于高精度单光子探测器的激光测距方法,实现了百微米量级精度的非合作目标激光测距.单光子测距系统引入参考位置,有效地抑制了系统延时漂移,使光子飞行时间测量精度达到0.5ps,在2m测距距离处,单光子测距系统的测距精度达到65μm@RMS.这项工作达到了当前脉冲飞行时间测距最高精度水 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于弹性变分模态提取的时间相关单光子计数信号去噪
    摘要:单光子激光雷达的回波信号具有极低的信噪比,有效地消除噪声和提取出回波信号特征是提升单光子激光雷达测距精度的关键,变分模态分解算法需要使用者依据经验确定分解本征模态函数数量,不具有适用性和通用性.为此,本文基于时间相关单光子计数信号特点,提出了在变分模态分解中让信号按照指定频率进行聚类分解的变分 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 一种基于束缚态的可调等离子体光子晶体窄带滤波器
    摘要:利用一维光子晶体和二维等离子体光子晶体构建了一种基于束缚态的可调窄带滤波器,滤波器的工作频率位于两个光子晶体的共同禁带内.使用COMSOLMultiphysics有限元仿真软件研究了一维光子晶体的几何参数和等离子体参数对滤波器性能的影响.研究发现两个禁带的中心频率和深度越接近,则滤波器的峰值透 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29

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

© 2020 FreeKaoYan! . All rights reserved.