搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

碟形玻色-爱因斯坦凝聚体中(0, 0, 2)剪刀模的朗道阻尼和频移

赵军亚 李晨旭 马晓栋

引用本文:
Citation:

碟形玻色-爱因斯坦凝聚体中(0, 0, 2)剪刀模的朗道阻尼和频移

赵军亚, 李晨旭, 马晓栋

Landau damping and frequency-shift of (0, 0, 2) scissors mode in a disc-shaped Bose-Einstein condensate

Zhao Jun-Ya, Li Chen-Xu, Ma Xiao-Dong
PDF
HTML
导出引用
计量
  • 文章访问数:  4332
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-05-03
  • 修回日期:  2019-09-21
  • 上网日期:  2019-11-26
  • 刊出日期:  2019-12-05

碟形玻色-爱因斯坦凝聚体中(0, 0, 2)剪刀模的朗道阻尼和频移

  • 1. 新疆师范大学物理与电子工程学院,乌鲁木齐 830054
  • 2. 新疆医科大学医学工程技术学院,乌鲁木齐 830011
  • 通信作者: 马晓栋, xdma07@aliyun.com
    基金项目: 国家自然科学基金(批准号: 100864006, 11264039)、新疆维吾尔自治区高校科研计划重点项目(批准号: XJEDU20141029)、新疆师范大学“物理学”特色专业和新疆师范大学“物理学”重点学科资助的课题

摘要: 应用哈特里-福克-博戈留波夫平均场理论近似和基于托马斯-费米近似的解析方法, 研究碟形玻色-爱因斯坦凝聚体中(0, 0, 2)剪刀模的朗道阻尼和频移, 计算阻尼系数和频移大小以及它们的温度依赖. 计算中, 在集体激发本征频移微扰关系中考虑元激发弛豫及其弛豫之间的正交关系以获得阻尼和频移的计算公式, 把凝聚体基态波函数取为高斯分布函数的一级近似以消除托马斯-费米近似中三模耦合矩阵元的发散. 采用与相关实验研究相同的粒子数、囚禁频率和各向异性参量, 理论计算结果与相关实验测量结果相符合. 由于理论的复杂性和计算的困难性, 在大多数基于平均场理论的单分量和两分量玻色-爱因斯坦凝聚集体激发阻尼和频移的研究中采用半经典近似, 把准粒子激发能谱看成是连续的来积分计算各个准粒子跃迁对阻尼和频移的贡献, 而本文和本文前期工作按分立的准粒子激发频谱计算阻尼或频移, 并在研究过程中提出了考虑元激发弛豫及弛豫之间正交关系的改进方法, 希望这种方法对今后的工作有一定参考价值.

English Abstract

    • 玻色-爱因斯坦凝聚(Bose-Einstein-conden-sation, BEC)的集体激发作为一个基本问题对于多体问题的研究显然是至关重要的. 元激发是统计与凝聚态物理学的基本研究内容之一, 集体激发的研究是捕陷玻色凝聚气体的实验和理论研究活动的主要兴趣领域之一[1-23]. 研究集体激发对于了解凝聚体的基态、热力学性质及其超流特性都有十分重要的意义, 阻尼和频移是BEC中集体激发的主要特征, 其中有四个实验[24-27]对BEC系统中激发模阻尼和频移有完整的温度依赖测量. 粒子间相互作用导致集体激发振幅的衰减(阻尼)和频率的改变(频移). 阻尼和频移的机制有朗道(Landau)和巴利耶夫(Beliaev)两种, 朗道机制是一个准粒子激发吸收一个集体激发变为另一个准粒子激发, 巴利耶夫机制是一个集体激发变为两个准粒子激发. 其中巴利耶夫机制在能级分立的系统中对低能量集体激发阻尼和频移的贡献很小.

      对于BEC动力学的研究方法主要有高温高密度下的二流体理论方法[28-31]和低温低密度下的平均场理论方法[32-48]. 其中大多数实验研究的是低温低密度系统.

      平均场理论框架下理论工作具体的探究方法也不尽相同, 有文献[3235]的微扰理论、文献[36,37]的格林函数方法、文献[38]的二阶量子场理论和文献[3944]的哈特里-福克-博戈留波夫(Hartree-Fock-Bogoliubov, HFB)理论. 文献[3440]都聚焦于对实验[27]中集体模阻尼和频移温度依赖的理论解释. 这些工作都有一个共同点, 在求和计算各个准粒子跃迁对集体激发阻尼和频移的贡献时, 都采用了半经典近似, 即把准粒子激发能谱看成是连续的而进行积分运算. 与文献[3440]不同, 文献[32,33,4144]按分立的准粒子激发频谱来求和计算朗道阻尼或频移. 其中文献[4144]是本文的前期工作, 文献[41]应用的方法与文献[32]相同, 而文献[4244]则在文献[32,33]方法的基础上进行了发展和改进.

      文献[32]在文献[34]的频率微扰理论公式中考虑元激发的弛豫而引入洛伦兹宽度$ \varDelta $, 得到计算集体激发朗道阻尼的公式

      $ {\gamma} = \sum\nolimits_{ij}{\gamma}_{ij}\dfrac{1}{\text{π}} \frac{{\varDelta}/2}{\left({\omega}_{0}+{\omega}_i-{\omega}_j\right)^{2} +\left({\varDelta}/2\right)^{2}}. $

      文献[32]计算阻尼的公式中含有洛伦兹宽度, 应用公式计算时, 不考虑阻尼强度大的准粒子共振跃迁、只考虑一部分阻尼强度小的准粒子背景跃迁对阻尼的贡献, 在阻尼$ {\gamma} $随洛伦兹宽度$ {\varDelta} $的缓变范围内取阻尼的计算值. 由于在某些系统中不出现阻尼随洛伦兹宽度缓变的情况, 文献[33]把文献[32]阻尼公式中的洛伦兹半宽认定为集体激发的阻尼

      $ {\varDelta}/2 = \displaystyle\sum\nolimits_{ij}{\gamma}_{ij} \dfrac{1}{\text{π}}\dfrac{{\varDelta}/2}{\left({\omega}_{0} + {\omega}_i -{\omega}_j\right)^{2}+\left({\varDelta}/2\right)^{2}} $

      用迭代的方法计算朗道阻尼. 文献[32,33]有一个共同点, 它们在考虑元激发的实际弛豫时, 把三个元激发(朗道机制中相互耦合的一个集体激发和两个准粒子激发)的弛豫简单地加在一起, 把$ {\varDelta}/2 $看作是加在一起的结果. 另外, 阻尼和频移是相伴的物理现象, 但文献[32]只计算了阻尼, 而文献[33]虽然阻尼和频移都计算了, 但在其计算中对阻尼和频移的物理考虑有比较大的差别.

      与文献[33]相同, 文献[4244]也把洛伦兹半宽认定为集体激发的阻尼, 也用迭代的方法计算阻尼. 但与文献[32,33]不同, 文献[4244]在考虑三个元激发的弛豫时, 还考虑了三个元激发弛豫的正交关系, 理由是三个元激发的本征函数是正交的. 应用文献[4244]提出的考虑元激发弛豫及其各弛豫间正交关系的方法, 不仅给出了阻尼的计算公式, 而且同时给出了频移的计算公式, 先用迭代的方法计算出阻尼, 再根据阻尼的计算结果计算频移, 分别得到了与实验相符或与相近实验对比合理的理论计算结果. 这些问题将在2.3节详细说明.

      总体上看, BEC中集体激发阻尼和频移的理论研究开展不足. 如前所述文献[3440]采用处理连续本征值元激发谱的半经典近似方法, 与实验[24-27]中本征值分立的轴对称系统有很大的差别. 迄今为止, 只有几个关于两分量玻色-爱因斯坦凝聚(two-component Bose-Einstein condensations, 2BECs)集体激发阻尼和频移的理论研究[45-48], 它们都是采用HFB平均场理论, 也都采用了半经典近似. 由于问题的复杂性, 计算球对称和轴对称系统中集体激发的阻尼和频移, 首先要计算大量分立的元激发本征函数和本征值, 数值模拟工作相当困难, 目前只有文献[32]这一个数值模拟工作, 而且研究的是较为简单的球对称系统. 这是理论工作开展不足的一个原因, 但这不是主要原因. 其主要原因是文献[3241,4548]在考虑元激发的弛豫时, 把三个元激发的弛豫简单地相加在一起, 这是值得商榷的. 因为如前所述三个元激发的本征函数是正交的, 例如在HFB平均场理论中, 元激发本征函数和本征频率是通过对角化系统巨正则哈密顿得到的.

      需要说明一点, 文献[4144]是解析研究工作, 在阻尼或频移的计算中, 分别应用了文献[22,23]的球对称和轴对称BEC元激发分立的本征函数和本征函数解析解, 而这种解析方法是在HFB平均场理论[39,40]框架下, 忽略正常和反常准粒子平衡密度, 并通过改进托马斯-费米近似(Thomas-Fermi approximation, TFA)得到的. 这个问题将在2.4节详细说明.

      本文研究蝶形BEC中剪刀模的朗道阻尼和频移, 计算中采用文献[27]实验的系统参量, 理论计算结果与实验的测量数据相对比, 进一步检验文献[4244]提出的考虑元激发弛豫及其正交关系的方法, 并进行更详细的说明, 希望这种方法能够得到推广和发展, 对研究BEC和2BECs中集体激发阻尼和频移的工作开展有一定的参考价值.

    • 采用基于HFB平均场理论的博戈留波夫-德热纳(Bogoliubov-de Gennes, BdG)方程组本征函数集[22,23]和集体激发本征频率微扰关系[39,40], 分别在2.1节和2.2节简要介绍.

      采用文献[4244]在集体激发本征频率微扰关系中考虑元激发弛豫及弛豫间正交关系而获得的阻尼和频移计算公式, 在2.3节详细讨论.

      采用文献[4244]在凝聚体基态波函数引入参量$ q = 1 $以消除TFA三模耦合矩阵元发散的方法, 在2.4节详细讨论.

    • 稀薄囚禁BEC原子气体在轴对称外势阱$ V_{\rm ext}\left({ r}\right) = m\omega_{\rm ho}^{2}r^{2}/2 $ ($ r^{2} = s^{2}+{\lambda}^{2}z^{2} $, $ s^{2} = x^{2}+y^{2}, $ m是原子质量, $ \omega_{\rm ho} $是横向谐振频率, $ \lambda $是各向异性参量, 即纵向横向谐振频率之比$ \lambda = \omega_{z} / \omega_{\rm ho} $)中的巨正则哈密顿量为

      $\begin{split} K =\, & H-\mu N = \int {\rm d} { r}\psi^{†}\left({ r},t\right)H_{0}\psi\left({ r},t\right)\\ &+\frac{g} {2}\int {\rm d}{ r}\psi^{†}\left({ r},t\right)\psi^{†}\left({ r},t\right)\psi\left({ r}, t\right)\psi\left({ r},t\right),\end{split} $

      其中$ H_{0} = -\hbar^{2}\nabla^{2}/\left(2 m\right)+V_{\rm ext}\left({ r}\right)-\mu $, $ \mu $是系统化学势、耦合常数$ g = 4{\text{π}}\hbar^{2}a_{\rm sc}/m $,$ a_{\rm sc} $是s波散射长度. 玻色场算符$ \psi({ r}, t) $的运动方程为

      $ {\rm i}\hbar\frac{\partial\psi\left({ r},t\right)}{\partial t} \!=\! H_{0}\psi\left({ r},t\right) \!+\!{g}\psi^{†}\left({ r},t\right)\psi\left({ r},t\right)\psi\left({ r},t\right). $

      玻色场算符$ \psi $可分解为凝聚部分$ \phi = \langle\psi\rangle $和非凝聚部分$ \widetilde{\psi} $, 即$ \psi({ r}, t) = \phi ({ r}, t)+\widetilde{\psi}\left({ r}, t\right) $, 凝聚部分又分解为静态部分$ \phi_{0} = \left\langle\psi\right\rangle_{0} $和激发态部分$ {\text{δ}} \phi $, 即$ \phi\left({ r}, t\right) = \phi_0\left({ r}\right) + {\text{δ}}\phi\left({ r}, t\right) $, 非凝聚部分的正常、反常热粒子密度$ \widetilde{n} = \big\langle\widetilde{\psi}^{†}\widetilde{\psi}\big\rangle $,$ \widetilde{m} = \big\langle\widetilde{\psi}\widetilde{\psi}\big\rangle $也分别分为静态部分$ \tilde{n}^{0} = \big\langle\tilde{\psi}^{†}\tilde{\psi} \big\rangle_{0} $,$ \tilde{m}^{0} = \big\langle\tilde{\psi}\tilde{\psi}\big \rangle_{0} $和激发部分$ {\text{δ}}\tilde{n} $,$ {\text{δ}}\tilde{m} $, 即$ \widetilde{n}({ r}, t) = \widetilde{n}^{0}({ r})+{\text{δ}}\widetilde{n}({ r}, t) $,$ \widetilde{m}({ r}, t) = \widetilde{m}^{0}({ r})+ $${\text{δ}}\widetilde{m}({ r}, t) $, 其中$ \langle\cdots\rangle $$ \langle\cdots\rangle_0 $分别代表非平衡平均和平衡平均.

      采用近似: 1) $ \langle\tilde{\psi}\rangle = 0 $; 2) $ \tilde{n}^{0} = 0 $, $ \tilde{m}^{0} = 0$; 3) $ \langle\tilde{\psi}\tilde{\psi}\tilde{\psi}\rangle = \langle\tilde{\psi}^{†}\tilde{\psi} \tilde{\psi}\rangle = 0 $; 4) $ \tilde{\psi}^{†}\tilde{\psi}^{†}\tilde{\psi}\tilde{\psi} = 4\tilde{n}\tilde{\psi}^{†}\tilde{\psi}+$ $ \tilde{m}\tilde{\psi}^{†}\tilde{\psi}^{†}+ \tilde{m}^{*}\tilde{\psi}\tilde{\psi} $.

      利用上述分解和近似, 可以得到凝聚部分波函数满足的方程

      $ \begin{split} & {\rm i}\hbar\frac{\partial\phi\left({ r},t\right)}{\partial t} \\=\, & H_{0}\phi\left({ r},t\right)+g|\phi\left({ r},t\right)|^{2}\phi\left({ r},t\right)\\ & +2g\phi\left({ r},t\right){{\text{δ}}}\widetilde{n} \left({ r},t\right)+g\phi^{*}\left({ r},t\right){\text{δ}} \widetilde{m}\left({ r},t\right).\end{split} $

      忽略(3)式中的${\text{δ}}\tilde{n} $${\text{δ}}\tilde{m} $可以得到凝聚静态部分基态波函数满足的格罗斯-皮塔耶夫斯基(Gross-Pitaevskii, GP)方程

      $ \left[H_{0}+gn_{0}\left({ r}\right)\right]\phi_{0}\left({ r}\right) = 0, $

      其中$ n_{0} = \phi_{0}^{2} $为凝聚体密度.

      采用博戈留波夫变换

      $ \widetilde{\psi}({ r}, t) = \sum\limits_{k} \left[u_{k}\left({ r} \right) \alpha_{k}^{†}(t) + v^{*}_{k}\left({ r}\right)\alpha_{k}(t)\right], $

      ($\displaystyle \int {\rm d} { r}\left[u_{k}^{*}\left({ r}\right)u_{k'}\left({ r}\right)-\right.$$\left. v_{k}\left({ r}\right)v^{*}_{k'}\left({ r}\right)\right] = \delta_{kk'} $ 为其归一化条件, 而玻色算符$ \alpha_{k} $满足$ \left[\alpha^{†}_{k}\right., $$\Big.\alpha_{k'}\Big] = \delta_{kk'} $), 并对角化哈密顿量(1), 得到非凝聚激发部分准粒子本征函数集满足的BdG方程组

      $ \begin{split} & \mathcal{L}u_{j}\left({ r}\right)+gn_{0}v_{j}\left({ r}\right) = \epsilon_{j}u_{j}\left({ r}\right), \\ & \mathcal{L}v_{j}\left({ r}\right)+gn_{0}u_{j}\left({ r}\right) = -\epsilon_{j}v_{j}\left({ r}\right), \end{split} $

      其中算子$ \mathcal{L} = H_{0}+2 gn_{0} $.

      由(3)式还可得到描述凝聚激发部分的运动方程

      $ \begin{split} & {\rm i}\hbar \frac{\partial{\text{δ}}\phi\left({ r},t\right)}{\partial t} \\= \,&\left[H_{0} +2gn_{0}\right]{\text{δ}}\phi\left({ r},t\right)+gn_{0}{\text{δ}}\phi^{*}\left({ r},t\right)\\ & +2g\phi_{0}{\text{δ}}\widetilde{n}\left({ r},t\right)+g\phi_{0}^{*}{\text{δ}}\widetilde{m}\left({ r},t\right), \end{split}$

      略去(6)式中$ {\text{δ}}\widetilde{n} $$ {\text{δ}}\widetilde{m} $, 容易看出凝聚体中的集体激发$ {\text{δ}}\phi\left({ r}, t\right) = u_{0}\left({ r}\right)\exp \left(-{\rm i}\epsilon_{0}t/\hbar\right) $$ {\text{δ}}\phi^{*}\left({ r}, t\right) = $$v_{0}\left({ r}\right)\exp\left(-{\rm i}\epsilon_{0}t/\hbar\right) $也决定于BdG方程, 即$ u_{0}\left({ r}\right) $,$ v_{0}\left({ r}\right) $$ \epsilon_{0} $满足于(5)式.

      通过无量纲化$ \overline{r} = r/R_{0} $,$ \overline{\nabla} = R_{0} \nabla $ ($R_{0} = $$ \sqrt{2\mu/m\omega_{\rm ho}^{2}} $是凝聚体特征半径)并引入$ \zeta = \hbar\omega_{h0}/(2\mu) $, (4)式变为无量纲形式

      $ \zeta^{2}\sigma \left(\overline{{ r}}\right) + \overline{r}^{2}-1+|\phi_{0}\left(\overline{{ r}}\right)/\phi_{0}\left(0\right)|^{2} = 0, $

      其中$ \sigma\left(\overline{{ r}}\right) = -\big[\overline{\nabla}^{2} \phi_{0} \left(\overline{{ r}}\right)\big]/\phi_{0}(0) $是一个正比于凝聚静态部分动能的量. 通过无量纲化元激发频率$ \overline{\omega} = \omega/\omega_{\rm ho} $ ($ \omega = \epsilon/\hbar $)并定义$ \phi^{\pm}_j = u_{j}\pm v_{j} $, (5)式变为无量纲形式

      $ \begin{split} & -\bar{\nabla}^{2} \left(1-\bar{r}^{2}\right)\phi^{+}_j -\left(1-\bar{r}^{2}\right)\sigma\phi^{+}_j \\ & +\frac{\zeta^{2}}{2} \left[\bar{\nabla}^{4}+3\bar{\nabla}^{2}\sigma+\sigma\bar{\nabla}^{2} +3\sigma^{2}\right] \phi^{+}_j = 2\bar{\omega}_{j}^{2}\phi^{+}_j,\\ & -\left(1-\bar{r}^{2}\right)\bar{\nabla}^{2}\phi^{-}_j-\left(1-\bar{r}^{2}\right)\sigma\phi^{-}_j\\ &+ \frac{\zeta^{2}}{2}\left[\bar{\nabla}^{4} + \bar{\nabla}^{2}\sigma+3\sigma\bar{\nabla}^{2} +3\sigma^{2}\right]\phi^{-}_j = 2\bar{\omega}_{j}^{2}\phi^{-}_j. \end{split} $

      (7)式改进的托马斯-费米近似(beyond Thomas-Fermi-approximation, bTFA)解为

      $ \phi_{0}\left(\bar{{ r}}\right) = \phi_{0} \left({0}\right) \left(1-\bar{s}^{2} - \lambda^{2}\bar{z}^{2}\right)^{(q+1)/2}, $

      其中$ \bar{s} = s/R_0 $, $ \bar{z} = z/R_0 $, $ \phi_{0}\left({0}\right) = \lambda\left[N_{0}/\left(2{\text{π}} R_0^{3}\times\right.\right. $$B\left.\left.\left(3/2, 2+q\right)\right)\right]^{1/2} $是归一化常数, 而$ B\left(3/2, 2+q\right) $是贝塔函数,$ N_0 = R_0^3 \displaystyle \int {\rm d} {\bar{ r}} |\varPhi_0 \left({\bar{ r}}\right)|^2 $是凝聚体粒子数, 化学势$ \mu \!=\! \hbar\omega_{\rm ho}\left[4\lambda P/B\left(3/2, 2+q\right)\right]^{2/5}\!/2, $$ P= $$N_{0}a_{\rm sc}/a_{\rm ho} $是原子相互作用强度且$ a_{\rm ho} = \sqrt{\hbar/m\omega_{\rm ho}} $是线性谐振子的特征半径.

      (8)式元激发本征函数的bTFA解为

      $\begin{split} \phi^{\pm}_{n_zn_{\rm s}m}\left(\bar{{ r}}\right) ={}& \frac{\left(\zeta\bar{\omega}_{n_zn_{\rm s}m}\right)^{\pm 1/2}}{\sqrt{2{\text{π}} R_0^{3}I_{n_zn_{\rm s}m}}}\left(1-\bar{s}^{2}-{\lambda}^{2} \bar{z}^{2}\right)^{(q\mp1)/2} \\ & \times \bar{s}^{m}P_{n_{\rm p}}^{(2n_{\rm s})} \left(\bar{s},\bar{z}\right) {\rm e}^{{\rm i} m\varphi},\\[-16pt]\end{split}$

      其中$ \bar{\omega}_{n_zn_{\rm s}m} $是本征值, $ I_{n_zn_{\rm s}m} $是归一化系数, $n_{\rm p} = $$ 0, 1, 2, \cdots $是主量子数, $ n_{\rm s} = 0, 1, 2, \cdots $, int[$ n_{\rm p}/2] $是径向量子数, $ n_{z} = n_{\rm p} - 2 n_{\rm s} $是轴向量子数, m = 0, $ \pm1 $, $ \pm2, \cdots $是方位角量子数. 径向和轴向耦合函数$ P_{n_{\rm p}}^{(2 n_{\rm s})} $满足一个二维微分方程, 解的形式为$ P_{n_{\rm p}}^{(2 n_{\rm s})} (\bar{s}, \bar{z}) = \displaystyle\sum\limits_{k = 0}^{n_{\rm p}} \sum\limits_{n_{\rm s} = 0}^{{\rm int} [{k}/2]}b_{k, n_{\rm s}}\bar{z}^{{k}-2 n_{\rm s}}\bar{s}^{2 n_{\rm s}} $. 元激发的本征值$ \bar{\omega}_{n_zn_{\rm s}m} $由如下封闭方程决定:

      $ \begin{split} -1 = \, &\dfrac{g_{0}}{1+\dfrac{g_{1}}{1+ \dfrac{g_{2}}{1+\cdots}}} \\ &\quad\quad\ddots\\ & 1+\dfrac{g_{Z-2}}{1+0}, \end{split} $

      其中

      $ \begin{split} g_{n_{\rm s}} =\, & \frac{-4\lambda^{2}\left(n_{\rm s}+1\right)\left(n_{\rm s}+|m|+1\right)\left(n_{\rm p}-2n_{\rm s}-1\right)\left(n_{\rm p}-2n_{\rm s}\right)}{2\gamma-2|m|\left(1+q\right)-4n_{\rm s}\left(n_{\rm s} +|m|+1+q\right) -\lambda^{2}\left(n_{\rm p}-2n_{\rm s}\right)\left(n_{\rm p}-2n_{\rm s}+1 +2q\right)}\\ & \times\big[2\gamma-2|m|\left(1+q\right)-4\left(n_{\rm s}+1\right)\left(n_{\rm s}+|m|+2+q\right) \\ &-\lambda^{2}\left(n_{\rm p}-2n_{\rm s}-2\right)\left(n_{\rm p}-2n_{\rm s}-1+2q\right)\big]^{-1}, \end{split} $

      $ Z = 1+ {\rm int} [n_{\rm p}/2] $$ \gamma = \bar{\omega}_{n_zn_{\rm s}m}^{2} $. 系数$ b_{k, n_{\rm s}} $由如下递推关系决定

      $ \begin{split} & 4 (n_{\rm s}+1) (n_{\rm s}+|m|+1) b_{k+2,n_{\rm s}+1} + (k-2n+2) (k-2n_{\rm s}+1) b_{k+2,n_{\rm s}}\\ =\, & 4\lambda^{2}(n_{\rm s}+1) (n_{\rm s}+|m|+1 )b_{k,n_{\rm s}+1}-t[2\gamma-2|m| (1+q)-4n_{\rm s} (n_{\rm s}+|m|+1+q)\\ & -\lambda^{2} (k-2n_{\rm s}) (k-2n_{\rm s}+1+2q)]b_{k,n_{\rm s}}+ (k-2n_{\rm s}+2)(k-2n_{\rm s}+1)b_{k,n_{\rm s}-1}. \end{split} $

      以上各式中的q是为了消除TFA三模耦合矩阵元的发散而引入的参量, 将在2.4节说明.

    • 通过博戈留波夫变换还可以得到反映${\text{δ}}\tilde{n} $$ {\text{δ}}\tilde{m} $随时间的演化方程

      $ \begin{split} & {\rm i}\hbar\frac{\partial f_{ij}\left(t\right)}{\partial t} = \big\langle\big[\alpha_{i}^{†}\alpha_{j}, K\big]\big\rangle,\\ &{\rm i}\hbar\frac{\partial g_{ij}\left(t\right)}{\partial t} = \left\langle\left[\alpha_{i}\alpha_{j},K\right]\right\rangle, \end{split} $

      其中$ f_{ij} = \big\langle\alpha_{i}^{†}\alpha_{j}\big\rangle-f_{i}^{0} \delta_{ij}$$ g_{ij} = \left\langle\alpha_{i}\alpha_{j} \right\rangle $分别是正常和反常准粒子分布函数而$ f_{j}^{0} = \big\langle\alpha^{†}_{j} \alpha_{j}\big\rangle_{0} = $$1/\left[\exp \left(2\zeta\overline\omega_{j}/\overline{T}\right)-1\right] $是准粒子平衡占有数, $ \overline{T} = $$k_{\rm B}T/\mu $是无量纲温度, $ k_{\rm B} $是玻尔兹曼常数.

      凝聚体中的集体模, 设其频率为$ \overline{\omega}_{0} $. 对(6)式和(14)式进行傅里叶变换可以得到朗道机制的集体模本征频率微扰关系

      $ \overline{\omega} = \overline{\omega}_{0}+\sum\limits_{ij} \overline{\gamma}_{ij}\frac{1}{\text{π}} \frac{1}{\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j}, $

      其中

      $ \bar{\gamma}_{ij} = 4{\text{π}}\left(4{\text{π}} a_{\rm sc}/a_{\rm ho}\right)^{2}|\overline{A}_{ij}|^2\left(f_{i}^{0}-f_{j}^{0}\right), $

      是以$ \omega_{\rm ho}^{2} $为标度的无量纲阻尼强度, 而

      $ \begin{split} \bar{A}_{ij}\, & =a^{3}_{\rm ho}\int {\rm d}{ r} \phi_{0} \big[u_{0}\left(u_{i}u^{*}_{j}+v_{i}v^{*}_{j}+v_{i}u^{*}_{j}\right) \\ & \quad +v_{0}\left(u_{i}u^{*}_{j}+v_{i}v^{*}_{j}+u_{i}v^{*}_{j}\right)\big]\\ & = \frac{\left[4P/B(3/2, 2+q)\right]^{1/10}}{8{\text{π}}\lambda^{9/10}\left[I_{0}I_{i}I_{j}\bar{\omega}_{0}\bar{\omega}_{i} \bar{\omega}_{j}\right]^{1/2}}\left(\frac{N_{0}}{P}\right)^{1/2}\\ & \quad \times\int_{0}^{1}{\rm d}x\int_{0}^{1}{\rm d}y W_{0}W_{i}W_{j}^{*}F\left(x,y\right), \end{split} $

      是以$ 1/a^{3}_{\rm ho} $为标度的无量纲三模耦合矩阵元. 在此定义

      $ x = \bar{s}^{2} ,\; y = \lambda\bar{z}/\sqrt{1-x}, \;W_j = \bar{s}^{m}P_{n_{\rm p}}^{\left(2n_{\rm s}\right)}\left(\bar{s},\bar{z}\right), $

      $ \begin{split} F\left(x,y\right) =\, & 3\xi^{2}\bar{\omega}_{0}\bar{\omega}_{i} \bar{\omega}_{j}\left(1-x\right)^{2q-1/2}\left(1-y^{2}\right)^{2q-1}\\ & + (\bar{\omega}_{0} \!+\! \bar{\omega}_{i} \!-\! \bar{\omega}_{j} ) (1-x)^{2q+3/2} (1 \!-\! y^{2})^{2q+1}. \end{split}$

    • 本节介绍本文前期工作[42-44]提出的考虑元激发弛豫及其正交关系的方法, 并进行更为详细的说明.

      在文献[32-48]中, 都有集体激发频率扰动公式

      $ \overline{\omega} = \overline{\omega}_{0}+\Delta\overline{\omega}- {\rm i}\overline{\gamma}_{0}, $

      其中$ \overline{\omega}_{0} $$ \overline{\omega} $分别是与外界驱动共振的集体激发频率和外驱撤销后扰动的集体激发频率, 而$ \Delta\overline{\omega} $$ \overline{\gamma}_{0} $分别是集体激发的频移和阻尼系数. 这个公式是唯象的, 当不考虑阻尼时, $ \overline{\omega}-\overline{\omega}_{0} $正好等于频移$ \Delta\overline{\omega} $, 而当不考虑频移时, $ \overline{\omega} = \overline{\omega}_{0}- {\rm i} \overline{\gamma}_{0} $, 这样集体激发随时间变化的部分$ {{\rm e}^{-{\rm i}{\overline{\omega} t}} = {\rm e}^{-{\overline{\gamma}_{0} t}}{\rm e}^{-{\rm i}\overline{\omega}_{0}t}} $, 从中可以看出$ \overline{\gamma}_{0} $是阻尼系数, 也就是说虽然与阻尼和频移有关的两项都一起写在(18)式中, 但是在实际理解和应用阻尼(频移)时, 需要令频移(阻尼)为零.

      以下依次考虑三个元激发的弛豫. 在计算集体激发的阻尼和频移时, 朗道机制的正向和反向跃迁都需要考虑, 所谓元激发的弛豫就是指元激发的衰减或增益, 在正(反)向跃迁中, 集体激发$ \overline{\omega}_{0} $和准粒子激发$ \overline{\omega}_{i} $是衰减(增益)的, 准粒子$ \overline{\omega}_{j} $是增益(衰减)的, 而增益和衰减分别对应于元激发振幅随时间变化因子$ {\rm e}^{\pm{\overline{\gamma}t}} $中的正负号.

      第一步根据上述唯象公式(18)式考虑集体激发$ \overline{\omega}_{0} $的弛豫, 考虑弛豫时不考虑频移, 把(18)式改写为$ \overline{\omega} = {\overline{\omega}_{0}}\pm{{\rm i}\overline{\gamma}_{0}} $替换(15)式右边的$ \overline{\omega}_{0} $, (15)式变为

      $ \begin{split} \overline{\omega} =\; & \overline{\omega}_{0}+\sum\limits_{ij} \overline{\gamma}_{ij}\frac{1}{\text{π}} \frac{\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j} {\left(\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j\right)^{2} +\overline{\gamma}_{0}^{2}}\\ &\pm {\rm i} {\left[\overline{\gamma}_{0}-\sum\limits_{ij} \overline{{\gamma}}_{ij}\frac{1}{\text{π}}\frac{\overline{\gamma}_{0}} {\left(\overline{\omega}_{0}+\overline{\omega}_i- \overline{\omega}_j\right)^{2}+\overline{\gamma}_{0}^{2}}\right]}. \end{split}$

      此时如果考虑频移, 则不考虑弛豫, 所以可以令(19)式的虚部为零得到

      $ \overline{\gamma}_{0} = \sum\limits_{ij} \overline{{\gamma}}_{ij}\frac{1}{\text{π}}\frac{{\overline{\gamma}_{0}}} {\left(\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j\right)^{2}+\overline{\gamma}_{0}^{2}}, $

      同时得到

      $ \Delta\omega = \overline{\omega}-\overline{\omega}_{0} = \displaystyle\sum\limits_{ij} \overline{{\gamma}}_{ij}\dfrac{1}{\text{π}}\dfrac{\overline{\omega}_{0} +\overline{\omega}_i-\overline{\omega}_j} {\left(\overline{\omega}_{0}+\overline{\omega}_i- \overline{\omega}_j\right)^{2}+\overline{\gamma}_{0}^{2}}.$

      第二步考虑准粒子$ \overline{\omega}_{i} $的弛豫, 将(20)和(21)式中的$ \overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j $替换为$ \left[\left(\overline{\omega}_{0}+\overline{\omega}_i -\overline{\omega}_j\right)^{2}+\right. $$\left.\overline{\gamma}_{i}^{2}\right]/ \overline{\gamma}_i $, 其中$ \bar{\gamma}_{i} $是准粒子$ \overline{\omega}_{i} $的朗道阻尼系数, 分别得到

      $ \overline{\gamma}_{0} =\displaystyle\sum\limits_{ij} \overline{{\gamma}}_{ij}\dfrac{1}{\text{π}}\dfrac{\overline{\gamma}_{0}}{ \dfrac{\left[\left(\overline{\omega}_{0}+\overline{\omega}_i- \overline{\omega}_j\right)^{2} +\overline{\gamma}_i^{2}\right]^{2}}{\overline{\gamma}_{i}^{2}} +\overline{\gamma}_{0}^{2}}$

      $\begin{split} & \overline{\omega}-\overline{\omega}_{0} =\\ & \displaystyle\sum\limits_{ij} \overline{{\gamma}}_{ij}\dfrac{1}{\text{π}} \frac{\dfrac{\left(\overline{\omega}_{0}+\overline{\omega}_i -\overline{\omega}_j\right)^{2} +\overline{\gamma}_{i}^{2}}{\overline{\gamma}_{i}} \dfrac{\overline{\omega}_{0}+\overline{\omega}_i- \overline{\omega}_j}{|\overline{\omega}_{0}+\overline{\omega}_i -\overline{\omega}_j|}} {\dfrac{\left[\left(\overline{\omega}_{0}+\overline{\omega}_i- \overline{\omega}_j\right)^{2} +\overline{\gamma}_{i}^{2}\right]^{2}}{\overline{\gamma}_{i}^{2}} +\overline{\gamma}_{0}^{2}}. \end{split} $

      第三步考虑准粒子$ \overline{\omega}_{j} $的弛豫, 将(22)式和(23)式中的$ \overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j $替换为$ \big[\left(\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j\right)^{2} $$ +\overline{\gamma}_{j}^{2}\big]/\overline{\gamma}_j $, 其中$ \bar{\gamma_{j}} $是准粒子$ \overline{\omega}_{j} $的朗道阻尼系数, 最终得到计算集体激发朗道阻尼系数的公式

      $ \begin{split} & \overline{\gamma}_{0} =\displaystyle\sum\limits_{ij} \overline{{\gamma}}_{ij}\frac{1}{\text{π}} \\ & \times\frac{\overline{\gamma}_{0}}{ \dfrac{\left\{\left[\left(\overline{\omega}_{0}+\overline{\omega}_i- \overline{\omega}_j\right)^{2} +\overline{\gamma}_{j}^{2}\right]^{2}+\overline{\gamma}_j^{2} \overline{\gamma}_{i}^{2}\right\}^{2}}{\overline{\gamma}_j^{4} \overline{\gamma}_{i}^{2}} +\overline{\gamma}_{0}^{2}}, \end{split} $

      以及计算集体模频移的公式

      $ \overline{\omega}-\overline{\omega}_{0} = \displaystyle\sum\limits_{ij} \overline{{\gamma}}_{ij}\frac{1}{\text{π}}\frac{\dfrac{\left[\left( \overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j\right)^{2} +\overline{\gamma}_{j}^{2}\right]^{2}+\overline{\gamma}_j^{2} \overline{\gamma}_{i}^{2}}{\overline{\gamma}_j^{2}\overline{\gamma}_{i}} \dfrac{\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j} {|\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j|}} {\dfrac{\left\{\left[\left(\overline{\omega}_{0}+\overline{\omega}_i -\overline{\omega}_j\right)^{2} +\overline{\gamma}_{j}^{2}\right]^{2}+\overline{\gamma}_j^{2} \overline{\gamma}_{i}^{2}\right\}^{2}}{\overline{\gamma}_j^{4} \overline{\gamma}_{i}^{2}} +\overline{\gamma}_{0}^{2}}.$

      第二步和第三步在考虑准粒子的弛豫时, 用$ \big[\left(\overline{\omega}_{0}+\overline{\omega}_i- \overline{\omega}\right)^{2}+\overline{\gamma}^{2}\big]/ \overline{\gamma} $替换$ \overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j $, 既考虑了准粒子的弛豫, 也考虑了各元激发弛豫之间的正交关系, 而不是把三个元激发的弛豫简单地直接加在一起的, 这是因为各元激发的本征函数是正交的.

      考虑元激发的弛豫, 弛豫使元激发能级产生宽度, 从而使跃迁的失谐增大. 把$ \overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j $替换为$ \big[\left(\overline{\omega}_{0}+\overline{\omega}_i- \overline{\omega}\right)^{2}+\overline{\gamma}^{2}\big]/\overline{\gamma} $, 相当于修改失谐, 修改后的失谐比原来的失谐大, 跃迁对阻尼和频移的贡献减小.

      如果$ \overline{\omega}_{0}\!+\!\overline{\omega}_i \!-\!\overline{\omega}_j\ll \overline{\gamma} $, $ \big[\left(\overline{{\omega}}_{0}+\overline{{\omega}}_i-\overline{{\omega}}_j \right)^{2}\!+\!\overline{{\gamma}}^{2}\big]/ $$\overline{\gamma}\rightarrow \overline{\gamma} $, 尽管跃迁的失谐很小, 但由于修改, 失谐接近于弛豫; 如果$ \overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j = \overline{\gamma} $, $ \big[\left(\overline{{\omega}}_{0}+\overline{{\omega}}_i- \overline{{\omega}}_j\right)^{2}+$$\overline{{\gamma}}^{2}\big]/\overline{\gamma} = 2\overline{\gamma} $, 修改后的失谐恰好是原来失谐的2倍; 如果$ \overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j\gg\overline{\gamma} $, $ {\left[\left(\overline{{\omega}}_{0}+\overline{{\omega}}_i -\overline{{\omega}}_j\right)^{2}+\overline{{\gamma}}^{2}\right]}/ $${\overline{\gamma}\gg{\omega}_{0}+\overline{{\omega}}_i-\overline{{\omega}}_j} $, 失谐修改后极大, 跃迁对阻尼贡献极小.

      每考虑一个元激发的弛豫, 都使失谐增大、阻尼和频移减小, 考虑三个元激发的弛豫, 最终结果使失谐更大、阻尼和频移更小.

      阻尼和频移的计算公式都是在集体激发的本征频率微扰关系(15)式中考虑元激发的弛豫及弛豫之间正交关系得到的, 而(15)式实际上就是不考虑任何元激发弛豫的频移公式. 在文献[42-44]的研究中发现这样一个规律, 即如果直接用(15)式计算集体激发的频移, 理论结果比实验结果大约大三个数量级, 如果依次用考虑一个、两个、三个元激发弛豫及其正交关系的频移公式(21)式、(23)式、(25)式计算, 理论结果依次大约减小一个数量级, 最后理论结果与实验结果相符. 阻尼的计算结果与上述考虑弛豫及其弛豫之间正交关系的依赖规律相类似.

      文献[42,43]分别研究雪茄形BEC中单极子模和蝶形BEC中四极子模的朗道阻尼和频移, 分别采用文献[26]和文献[27]实验的系统参量, 理论计算结果分别与文献[26]和文献[27]实验的测量数据相符.

      文献[32]和文献[44]研究的是系统参量完全相同的球对称BEC, 都计算了单极子的朗道阻尼. 文献[44]考虑了所有贡献大的准粒子跃迁, 计算结果比轴对称系统的实验结果[24-27]小而合理, 这是因为球对称系统的能级比轴对称的稀疏, 失谐小的准粒子跃迁比较少. 而文献[32]如前所述只考虑一部分准粒子跃迁, 计算结果却比轴对称系统的实验结果大, 这是因为文献[32]在考虑元激发的弛豫时, 直接把三个元激发弛豫简单地加了起来而没有考虑弛豫之间的正交关系.

    • 囚禁势中BEC基态波函数一般应该取高斯分布函数(谐振子基态波函数)形式$ \phi_{0}\sim {\rm e}^{-\bar{r}^{2}} $, 但把它代入BdG方程(8)式得不到元激发本征函数集的解析解.

      忽略GP方程(7)式中的凝聚静态部分动能$ \sigma\left(\bar{{ r}}\right) $, 可得TFA基态波函数

      $ \phi_{0}' \left(\bar{{ r}}\right) = \phi_{0}' \left({0}\right)\left(1-{\bar{r'}}^{2}\right)^{1/2}, $

      其中 $ \bar{r}' = r/R_{0}^{\prime} $, 归一化常数$\phi_{0}^{\prime}\left({0}\right) = $$ \lambda\left[N_{0}/\left(2{\text{π}} R_0^{\prime3}B\left(3/2, 2\right)\right)\right]^{1/2} $而凝聚体特征半径$ R_{0}^{\prime} = \sqrt{2\mu^{\prime}/({m}\omega^{2}_{h0})} $, 化学式$ \mu^{\prime} = \hbar\omega_{\rm ho}\left[4\lambda P/\right. $$\left. B\left(3/2, 2\right)\right]^{2/5}/2 $. 把(26)式代入BdG方程(8)式可以得到元激发本征函数集的解析解[22], 但由这些解析解计算出的三模耦合矩阵元$ \bar{A}_{ij} $是发散的.

      在TFA基态波函数(26)式中引入参量$ q>0 $, 得到bTFA基态波函数(9)式, 把(9)式代入BdG方程(8)式也可得到本征函数集的解析解[23], 而由此计算的$ \bar{A}_{ij} $不发散.

      文献[42-44]取$ q = 1 $, 以下从两个方面说明其中的原因.

      一方面, 取$ q = 1 $, bTFA基态波函数(9)式为谐振子基态波函数(高斯分布函数)的一阶近似形式$ \phi_{0}\sim 1-\bar{r}^{2} $, 当$ \bar{r} \rightarrow 0 $时, $ \phi_{0}\rightarrow {\rm e}^{-\bar{r}^{2}} $, 即越靠近BEC中心, 基态波函数越接近高斯分布函数, 囚禁势中BEC粒子保留着一定的谐振子特征.

      另一方面, 在凝聚体边界$ 1-\bar{s}^{2}-\lambda^{2}\bar{z}^{2} = 0 $处, 凝聚体的动能密度$ \sim-\phi_{0} \overline{\nabla}^{2} \phi_{0}\sim-(1+q)(1-\bar{s}^{2}-$$\lambda^{2}\bar{z}^{2})^{-1+q} [-2-\lambda^{2}+\bar{z}^{2}(2\lambda^{2}+q\lambda^{4}) +\bar{s}^{2}(1+q+\lambda^{2})] $$ 0<q<1 $是发散的, 在$ q = 1 $时不确定, 在$ q>1 $时等于0. 因此, 取$ q = 1+\epsilon $ ($ \epsilon $是无穷小量)是BEC基态波函数满足连续性边界条件时q最小的一种情况.

      因为TFA元激发本征值的理论计算值接近于实验值, 采用TFA的元激发本征值. 因为TFA的元激发本征函数会导致三模耦合矩阵元的发散, 采用bTFA的元激发本征函数.

      如前所述, 文献[44]与文献[32]均研究系统参量完全相同的球对称BEC. 文献[44]还进行了另一种计算, 也用上述解析方法, 但在考虑元激发的弛豫上与文献[32]方法相同, 即在本征频率微扰关系(15)式中引入洛伦兹宽度, 也只考虑一部分阻尼强度小的准粒子背景跃迁对阻尼的贡献, 在阻尼随洛伦兹宽度的缓变范围内取阻尼的计算值, 阻尼的解析计算结果与文献[32]的数值模拟计算结果相符合. 通过这样一个计算对取$ q = 1 $进行了验证.

    • 考虑频率为$ \omega_{\rm ho} = 810 $ Hz, 各向异性参数$ \lambda = 2.83 $的谐振囚禁势中原子数$ N = 6000 $的碟形87Rb原子气体($ a_{\rm sc} = 5.82\times10^{-9} $ m). 研究的集体模是剪刀模$ \left(n_{z}, n_{\rm s}, m\right) = \left(0, 0, 2\right) $, 其频率$ {\omega}_{0} $ = 1.414 (以$ \omega_{\rm ho} $为单位), 其集体模的博戈留波夫振幅$ u_{0} = \left(\phi^{+}_{002}+\phi^{-}_{002}\right)/2 $$ v_{0} = \left(\phi^{+}_{002}-\phi^{-}_{002}\right)/2 $. 以上集体模和所有参量的选取均与实验[27]相同. BEC相变临界温度为$ T_{\rm c} = 4.5\left(\omega_{\rm ho}/200{\text{π}}\right) \left(\lambda N\right)^{1/3} $ nK = $ 149.2 $ nK[1].

    • 阻尼强度采用(16)式和(17)式计算. 在图1图2及其插图中, 给出了各个准粒子跃迁$ \omega_{ij} $ ($ \omega_{ij} = \omega_{j}-\omega_{i} $, 以$ \omega_{\rm h0} $为单位)的阻尼强度$ \gamma_{ij} $ (以$ \omega_{\rm h0}^{2} $为单位). 这些准粒子跃迁服从剪刀模选择条件$ m_{i}-m_{j} = 2 $ (其中$ m_{i} $$ m_{j} $分别是准粒子$ \omega_{i} $$ \omega_{j} $的方位角量子数). 图1图2及其插图中的箭头指向集体模的频率$ \omega_{\rm o} $, 竖线的位置对应允许的跃迁频率$ \omega_{ij} $, 其高度给出$ \gamma_{ij} $的计算值.

      图  1  以跃迁频率$\gamma_{ij}$为变量的剪刀模阻尼强度$\omega_{ij}$函数线状图

      Figure 1.  Histogram of damping strength $\gamma_{ij}$ as a function of the transition $\omega_{ij}$ for the scissors mode in the condensate.

      图  2  以跃迁频率$ \gamma_{ij} $为变量的剪刀模阻尼强度$ \omega_{ij} $函数线状图 ($ \omega_{ij}/\omega_{\rm ho} $取值范围为0—1.6)

      Figure 2.  Histogram of damping strength $ \gamma_{ij} $ as a function of the transition $ \omega_{ij} $ for the scissors mode in the condensate, ($ \omega_{ij}/\omega_{\rm ho} $= 1–1.6).

      图1和其插图中的跃迁频率分别满足$ -0.6\omega_{\rm ho}<\omega_{ij}<3.4\omega_{\rm ho} $$ 1.21\omega_{\rm ho} < \omega_{ij}<1.61\omega_{\rm ho} $, 图2和其插图中的跃迁频率分别满足$1.21\omega_{\rm ho}< $$ \omega_{ij}<1.61\omega_{\rm ho} $$ 1.396\omega_{\rm ho}<\omega_{ij}<1.432\omega_{\rm ho} $.

      由阻尼强度计算公式(16)式可以看出, 当跃迁频率$ \omega_{ij}<0 $ (即$ \omega_{j}<\omega_{i}) $时, 阻尼强度$ \gamma_{ij}<0 $; 当跃迁频率$ \omega_{ij} = 0 $ (即$ \omega_{j} = \omega_{i}) $时, 阻尼强度$ \gamma_{ij} = 0 $; 当跃迁频率$ \omega_{ij}>0 $ (即$ \omega_{j}>\omega_{i}) $时, 阻尼强度$ \gamma_{ij}>0 $.

      能级量子数大的准粒子本征函数振荡很快且极大值远离凝聚中心, 另外量子数大能级能量高, 其玻色布居因子$ f^{0}_{i} $小, 因此量子数大能级之间跃迁的三模耦合矩阵元小. 忽略对阻尼和频移贡献小的跃迁, 在计算中选取($ 0, 0, m $), ($ 1, 0, m $), ($ 0, 1, m $), ($ 2, 0, m $), ($ 1, 1, m $), ($ 3, 0, m $), ($ 0, 2, m $), ($ 2, 1, m $), ($ 4, 0, m $), ($ 1, 2, m $), ($ 3, 1, m $), ($ 5, 0, m $), ($ 0, 3, m $), ($ 2, 2, m $), ($ 4, 1, m $)和($ 6, 0, m $) 而m满足$ m\leqslant 40 $的能级. 上述所选能级之间的跃迁频率范围在图1跃迁频率取值范围$ -0.6\omega_{ij} < \omega_{ij} < 3.4\omega_{0} $之内.

    • 将(24)式右边的$ \bar{\gamma_{0}} $替换成$ \bar{\gamma} $并近似地认为阻尼系数正比于频率, 把(24)式改写为

      $ \begin{array}{l} \overline{\gamma}_{0} =\displaystyle\sum\limits_{ij} \overline{\gamma}_{ij} \dfrac{1}{\text{π}} \dfrac{\overline{\gamma}}{\dfrac{\left\{\left[\left(\overline{\omega}_{0} +\overline{\omega}_i-\overline{\omega}_j\right)^{2}+\left(\dfrac{\overline{\omega}_j} {\overline{\omega}_{0}}\overline{\gamma}\right)^{2}\right]^{2} +\left(\dfrac{\overline{\omega}_j}{\overline{\omega}_{0}} \overline{\gamma}\right)^{2} \left(\dfrac{\overline{\omega}_i}{\overline{\omega}_{0}} \overline{\gamma}\right)^{2}\right\}^{2}} {\left(\dfrac{\overline{\omega}_j}{\overline{\omega}_{0}} \overline{\gamma}\right)^{4}\left(\dfrac{\overline{\omega}_i} {\overline{\omega}_{0}}\overline{\gamma}\right)^{2}} +\overline{\gamma}^{2}}, \end{array} $

      用迭代的方法计算集体模的阻尼系数. 图3给出$ \gamma_{0} $$ \gamma $变化的图($ \gamma_{0} $$ \gamma $都以$ \omega_{h0} $为单位), 图中可以看出, $ \gamma_{0} $$ \gamma $的增加而减小, 由星号表示三个点处$ \gamma_{0} = \gamma $, 温度$ T = 80 $, 100, 120 nK时的集体模阻尼系数值在此三点给出, 分别为0.02263, 0.03387, 0.05546, 相当于18.33, 27.43, 44.92$ {\rm s}^{-1} $.

      图  3  朗道阻尼系数$ \gamma_{0} $$ \gamma $变化函数图

      Figure 3.  The $ \gamma_{0} $ as a function of $ \gamma $ for the Landau damping rate.

      把(25)式改写为

      $ \begin{array}{l} \overline{\omega}-\overline{\omega}_{0} =\displaystyle\sum\limits_{ij} \overline{{\gamma}}_{ij}\dfrac{1}{\text{π}}\dfrac{\dfrac{\left[\left(\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j\right)^{2} +\left(\dfrac{\overline{\omega}_j}{\overline{\omega}_{0}}\overline{\gamma}_{0}\right)^{2}\right]^{2} +\left(\dfrac{\overline{\omega}_j}{\overline{\omega}_{0}}{\overline{\gamma}_{0}}\right)^{2}\left(\dfrac{\overline{\omega}_i}{\overline{\omega}_{0}}{\overline{\gamma}_{0}}\right)^{2}} {\left(\dfrac{\overline{\omega}_j}{\overline{\omega}_{0}}{\overline{\gamma}_{0}}\right)^{2}\dfrac{\overline{\omega}_i}{\overline{\omega}_{0}}{\overline{\gamma}_{0}}} {\dfrac{\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j}{|\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j|}}} {\dfrac{\left\{\left[\left(\overline{\omega}_{0}+\overline{\omega}_i-\overline{\omega}_j\right)^{2} +\left(\dfrac{\overline{\omega}_j}{\overline{\omega}_{0}}{\overline{\gamma}_{0}}\right)^{2}\right]^{2} +\left(\dfrac{\overline{\omega}_j}{\overline{\omega}_{0}}{\overline{\gamma}_{0}}\right)^{2} \left(\dfrac{\overline{\omega}_i}{\overline{\omega}_{0}}{\overline{\gamma}_{0}}\right)^{2}\right\}^{2}} {\left(\dfrac{\overline{\omega}_j}{\overline{\omega}_{0}}{\overline{\gamma}_{0}}\right)^{4}\left(\dfrac{\overline{\omega}_i}{\overline{\omega}_{0}}{\overline{\gamma}_{0}}\right)^{2}} +{\overline{\gamma}_{0}}^{2}}. \end{array} $

      并用阻尼系数$ \overline{\gamma}_{0} $的计算结果来计算频移大小. 图4给出了阻尼系数$ \gamma_{0} $ (单位为s–1)和频率$ \omega $ (单位为$ \omega_{\rm h0} $)随温度T (单位为nK)变化的关系. 图中平滑的曲线是我们的理论计算结果, 实心圆点和误差线表示实验值. 理论计算得出的阻尼系数和频移温度依赖结果都与文献[27]的实验结果相符.

      图  4  凝聚体集体模的频移(a)和朗道阻尼系数(b)随温度T变化

      Figure 4.  The frequency-shift (a) and the Landau damping rate $ \gamma_{0} $ (b) of the collective mode in the condensate as a function of the temperature T.

    • 如前所述, 采用元激发本征频率的TFA理论计算值, 因为它们接近于实验测量值. 图4(a)中, 理论曲线和实验数据点及其误差线错开了一段距离, 这是因为剪刀模本征频率的TFA理论计算值与实验测量值之间存在大约$ 5\% $的误差.

      如前所述, TFA基态波函数(26)式是GP方程(7)式的解, 而bTFA基态波函数(9)式是在(26)式中引入q而得到的, 因此(9)式相当于(26)式有误差, 其平均误差为

      $\begin{split} & \bar{\delta} = \\ & \frac{\displaystyle\int_{0}^{1}\!\!\bigg[\dfrac{\phi_{0}}{\phi^{\prime}_{0}} \Big(\! 1 \!-\! \dfrac{R^{\prime}_{0}}{R_{0}}\bar{r}^{\prime 2}\! \Big)^{\textstyle\frac{1+q}{2}} \!-\! (1\!-\! \bar{r}^{\prime2})^{\frac12} \bigg] (1-r^{\prime 2})^{\frac12}\bar{r}^{\prime2} {\rm d} \bar{r}^{\prime}} {\displaystyle\int_{0}^{1}\left(1-\bar{r}^{\prime 2}\right)^{\frac12}\bar{r}^{\prime 2}{\rm d} \bar{r}^{\prime}}. \end{split} $

      $ q = 1 $, 将本文选取的相关系统参量代入(29)式计算得其平均误差$ \bar{\delta} = 0.77\% $, 说明bTFA基态波函数(9)式近似满足GP方程(7)式.

      如前所述, 只有失谐很小的跃迁才对阻尼和频移有显著的贡献. 本文计算频移时考虑在图1跃迁频率范围内的跃迁, 在此范围外的跃迁对频移的贡献可以忽略不计. 图2的跃迁频率范围比图1小而与图1插图相同, 本文计算阻尼时考虑这个范围内的跃迁, 在此范围外的跃迁对阻尼的贡献可以忽略不计. 图2插图的跃迁频率取值范围只是图2$ 9\% $, 但在图2插图跃迁频率取值范围内的跃迁对阻尼的贡献是总阻尼的$ 50\% $.

    • 在HFB平均场理论近似的朗道机制频率扰动关系中, 考虑了三个元激发的弛豫及三个弛豫的正交关系, 得到计算BEC中集体激发朗道阻尼和频移的公式. 把凝聚体基态波函数取为高斯分布函数的一阶近似来消除TFA近似三模耦合矩阵元的发散而实施了解析计算. 采用同样的计算方法, 本文和两个前期工作分别研究BEC中剪刀模、单极子模和四极子模, 分别采用与相关实验相同的参量, 所有集体激发阻尼和频移的理论结果分别与相关实验结果相符合. 期望本文方法能够有助于BEC和2BECs中集体激发阻尼和频移的研究工作, 有助于探究元激发阻尼和频移的物理机制.

参考文献 (48)

目录

    /

    返回文章
    返回