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

两体相互作用费米系统在自旋轨道耦合和塞曼场中的基态转变

本站小编 Free考研考试/2021-12-29

摘要:在超冷费米系统中实现人造规范势的突破, 吸引了许多新问题的研究, 展现了许多新奇的物理现象. 本文研究了在环阱中, 具有自旋轨道耦合和塞曼作用的两体相互作用费米模型. 通过平面波展开的方法, 解析求解了两体费米系统的本征能态. 系统的总动量为守恒量, 可以在不同总动量空间中研究能谱. 研究发现: 随着塞曼相互作用增大, 在不同总动量空间, 两体费米系统的本征能量均逐渐降低, 系统基态从总动量为零空间转变到有限值空间. 从吸引到排斥相互作用, 无塞曼相互作用时, 基态总动量始终为零, 有塞曼相互作用时, 基态总动量从零转变为有限值. 通过单粒子和基态动量分布研究, 本文直观地揭示了由塞曼能级劈裂引起的基态转变.
关键词: 两体费米系统/
自旋轨道耦合/
塞曼相互作用

English Abstract


--> --> -->
自旋轨道耦合是粒子自旋内禀自由度和轨道自由度之间的相互作用, 在原子分子物理、量子光学和凝聚态物理中诱导产生了许多重要的物理现象[1-3]. 冷原子系统可以实现原子及其相互作用的精确操控, 是研究凝聚态系统及复杂物理系统的理想平台. 虽然原子是中性的且没有电极和磁极, 但仍可利用激光与原子的相互作用, 在超冷原子中实现人造自旋轨道耦合, 最早美国国家标准与技术研究院Speilman研究组[4]在超冷玻色气体中实现了一维自旋轨道耦合, 山西大学Zhang研究组[5]和美国麻省理工学院Zwierlein研究组[6]均实现了一维自旋轨道耦合费米气体, 山西大学Zhang研究组[7]进一步实现了二维自旋轨道耦合费米气体, 中国科学技术大学的研究组[8]在超冷玻色气体中实现了二维自旋轨道耦合, 这些突破为研究自旋轨道耦合系统提供了新的思路.
自旋轨道耦合系统的研究展现了许多新的物理现象[9-33]. 研究发现, 自旋轨道耦合费米系统有很多新奇的量子态, 例如拓扑超流态、Fulde-Ferell(FF)态等[9-16]. 考虑塞曼场时, 自旋轨道耦合和塞曼场共同作用诱导产生了非平庸的拓扑相[17-22]. 少体系统的研究为多体系统探索提供了直观的物理图像[23-33]. 二维自旋轨道耦合两体系统研究预示了多体系统可能出现FF配对态[23]. Rashba自旋轨道耦合两体精确解揭示了一种新的束缚态Rashbon, 在弱相互作用下仍能诱导粒子从弱配对态到玻色爱因斯坦凝聚态转变[24,25]. 在强相互作用下, 谐振子阱中自旋轨道耦合两体精确解发现反宇称态是玻色系统基态的重要组成部分[26]等. 该类研究为进一步探索多体系统物理现象奠定了基础.
本文研究了在一维环阱中, 具有自旋轨道耦合和塞曼场的两体相互作用费米模型. 在冷原子系统中, 通过激光耦合173Yb原子的两个超精细基态, 实现了有效的一维自旋轨道耦合和塞曼相互作用的两分量费米系统[34]. 同时, 通过在另外两个方向上施加强势阱囚禁可以得到有效的一维原子相互作用[35]. 对于两体问题, 通常可以将两原子坐标转换为质心和相对坐标. 质心运动不依赖于粒子间相互作用, 质心动量为守恒量, 可以在不同质心动量空间分析能谱. 通过将质心和相对运动分开, 许多两体系统能够被约化为单体问题[36].
本文研究的两体费米系统, 由于相对运动与质心运动耦合使得精确求解更为困难. 本文利用平面波展开法, 获得了解析的本征能量和本征态, 实现了两体费米系统随着自旋轨道耦合、塞曼相互作用和接触相互作用系数变化时物理性质的精确研究.
在环阱中, 考虑一个具有自旋轨道耦合、塞曼相互作用和粒子间接触相互作用的费米系统, 哈密顿量表示为
$ \begin{split}\hat{{H}} =\;&\!\! \int {{\rm{d}}x {\hat{{a}}} _x^\dagger }\! \left[ { - \frac{{{\hbar ^2}}}{{2m}}\frac{{{\partial ^2}}}{{\partial {x^2}}} + \alpha {{{\hat {{\sigma}}} }^x}\left( { - {\rm{i}}\frac{\partial }{{\partial x}}} \right) + \beta {{{\hat{{\sigma}}} }^z}} \right]{{\hat {{a}}} _x}\\&+ g\int {{\hat{{a}}} _{x \uparrow }^\dagger } {\hat {{a}}} _{x \downarrow }^\dagger {{\hat {{a}}} _{x \downarrow }}{{\hat {{a}}} _{x \uparrow }}{\rm{d}}x, \\[-15pt] \end{split}$
其中, ${{\hat{{a}}} _x} = {\big( {{{{\hat{a}} }_{x \uparrow }}, {{{\hat{a}} }_{x \downarrow }}} \big)^{\rm{T}}}$是两分量的场算符, 在环阱中, 算符需满足${{\hat{{a}}} _x} = {{\hat{{a}}} _{x + L}}$, 其中L是环阱周长; 第一项表示单粒子哈密顿量, 其中m为粒子质量, ${{\hat{{\sigma}}} ^j}$为泡利矩阵算符, $\alpha $为自旋轨道耦合系数, $\beta $为塞曼相互作用系数; 第二项表示粒子间的接触相互作用, g为相互作用系数.
两体费米系统的波函数表示为
$\begin{split} & {{\varPsi}} \left( {{x_1},{x_2}} \right) =\\& \int \big[ {\psi _{ \uparrow \uparrow }} (x_1,x_2) a_{{x_1} \uparrow }^\dagger a_{{x_2} \uparrow }^\dagger + \psi _{ \uparrow \downarrow } (x_1,x_2) a_{{x_1}\uparrow }^\dagger a_{{x_2} \downarrow }^\dagger \\ & + \psi_{\downarrow \uparrow} (x_1,x_2)a_{{x_1} \downarrow }^\dagger a_{{x_2} \uparrow }^\dagger + \psi _{ \downarrow \downarrow} ( x_1,x_2) a_{{x_1} \downarrow }^\dagger a_{{x_2} \downarrow}^\dagger \big]\\ & \times {\rm{d}}{x_1}{\rm{d}}{x_2} |0\rangle. \end{split} $
将哈密顿量(1)和波函数(2)代入薛定谔方程${\hat{{H}}} {{\varPsi}} \left( {{x_1}, {x_2}} \right) = E{{\varPsi}} \left( {{x_1}, {x_2}} \right)$, 两体费米系统的方程约化为$ H{{\psi}} \left( {{x_1}, {x_2}} \right) = E{{\psi}} \left( {{x_1}, {x_2}} \right) $, 其中,
$\begin{split}{{H}} = \;&\left( {\sum\limits_{i = 1,2} { - \frac{{{\hbar ^2}}}{{2{m_i}}}\frac{{{\partial ^2}}}{{\partial x_i^2}} - {\rm{i}}\alpha \sigma _i^x\frac{\partial }{{\partial {x_i}}} + \beta \sigma _i^z} } \right) \\&+ g\delta \left( {{x_1} - {x_2}} \right)\delta_{\sigma\bar\sigma},\end{split}$
其中, $ \sigma $$ \bar\sigma $表示自旋$ \uparrow $$ \downarrow $$ \bar\sigma\ne \sigma $; ${{\psi}} ( {{x_1}, {x_2}} ) = $$ ( {\psi _{ \uparrow \uparrow }}( {{x_1}, {x_2}} ), {\psi _{ \uparrow \downarrow }}( {{x_1}, {x_2}} ), {\psi _{ \downarrow \uparrow }}({x_1}, {x_2}), {\psi _{\downarrow \downarrow}}({x_1}, {x_2}))^{\rm{T}}$.
在环阱中, 两体费米系统的波函数需满足周期性边界条件, 表示为
$ {{\psi}} \left( {{x_1},{x_2}} \right) = {{\psi}} \left( {{x_1} + L,{x_2}} \right) = {{\psi}} \left( {{x_1},{x_2} + L} \right). $
由于自旋轨道耦合和塞曼相互作用, 自旋分量不再是好量子数, 不能在某一个自旋分量空间仅求解波函数的坐标部分, 而需要同时考虑自旋分量与坐标之间的耦合, 考虑所有自旋态时的本征波函数.
引入质心坐标$X = \left( {{x_1} + {x_2}} \right)/2$和相对坐标$x = {x_1} - {x_2}$, 哈密顿量(3)重新表示为
$\begin{split}{{H}} =\;& - \frac{{{\hbar ^2}}}{{2m}}\left( {\frac{1}{2}\frac{{{\partial ^2}}}{{\partial {X^2}}} - 2\frac{{{\partial ^2}}}{{\partial {x^2}}}} \right) - {\rm{i}}\alpha {\sigma ^X}\frac{\partial }{{\partial X}} \\&-{\rm{i}}\alpha \frac{{\sigma _1^x - \sigma _2^x}}{2}\frac{\partial }{{\partial x}} + \beta {\sigma ^Z} + g\delta \left( x \right)\delta_{\sigma\bar\sigma},\end{split}$
其中${\sigma ^X} = \sigma _1^x + \sigma _2^x$${\sigma ^Z} = \sigma _1^z + \sigma _2^z$. 在质心及相对坐标空间, 两体费米模型仍满足周期性边界条件, 该条件限定了波函数中动量的取值范围. 由(5)式的哈密顿量表达式可以看出, 自旋轨道耦合项与质心坐标和相对坐标都有关系, 而粒子间接触相互作用仅与相对运动有关, 与质心运动无关. 质心动量${\hat{{K}}} = - {\rm{i}}\dfrac{\partial }{{\partial X}}$为两费米子的总动量, 与哈密顿量对易, 是好量子数, 两费米子系统的本征能谱能够在不同总动量空间中求解.
将本征波函数用平面波展开, 表示为
${{\psi}} \left( {X,x} \right) = A\sum\limits_k {{{\rm{e}}^{{\rm{i}}KX + {\rm{i}}kx}}} \left( {\begin{array}{*{20}{c}} {{\psi _{ \uparrow \uparrow }}\left( {K,k} \right)} \\ {{\psi _{ \uparrow \downarrow }}\left( {K,k} \right)} \\ {{\psi _{ \downarrow \uparrow }}\left( {K,k} \right)} \\ {{\psi _{ \downarrow \downarrow }}\left( {K,k} \right)} \end{array}} \right),$
其中A是归一化系数, K为总动量, k是相对动量的平面波分量. 取$\sigma $$\sigma '$表示自旋$ \uparrow $$ \downarrow $, 波函数每一个分量${\psi _{\sigma \sigma '}}\left( {K, k} \right)$表示不同自旋态时平面波的展开系数. 通过求解薛定谔方程, 获得${\psi _{\sigma \sigma '}}\left( {K, k} \right)$的解析表达式, 以下将${\psi _{\sigma \sigma '}}\left( {K, k} \right)$简写为${\psi _{\sigma \sigma '}}$.
求解波函数之前, 首先通过周期性边界条件获得波函数中Kk的范围. 将${{\psi}} \left( {X, x} \right)$中的Xx变换到${x_1}, {x_2}$, 获得波函数$\psi \left( {{x_1}, {x_2}} \right)$表达式, 依据周期性边界条件(4)式, 总动量K和相对动量分量k需满足
${K_n} = \frac{{{\rm{2}}{\rm{\pi }}n}}{L},\tag{7a}$
$k_n^{n'} = \left\{ {\begin{aligned}&{\frac{2n'{\rm{\pi }}}{L}},&{n{\text{为偶数}},}\\&{\frac{(2n'{+1}){\rm{\pi }}}{L}},&{n{\text{为奇数}},}\end{aligned}} \right.\tag{7b}$
(6)式中的Kk分别用${K_n}$$k_n^{n'}$表示, 其中n$n'$是整数, n表示了总动量的量子数, 而$n'$描述了相对动量分量的量子数. 从(7b)式得出, $k_n^{n'}$与总动量量子数n的奇偶有关, 是一系列离散的数值. 对于给定的n, 即给定总动量${K_n}$, 相应的本征能量E表示为${E_n}$.
为了解析求解系统的本征波函数, 引入接触势$\delta \left( x \right) = \dfrac{1}{L}\sum\nolimits_{n'} {{{\rm{e}}^{{\rm{i}}k_n^{n'}x}}}$, 将波函数(6)式代入薛定谔方程${{H}}{{\psi}} \left( {X, x} \right) = {E_n}{{\psi}} \left( {X, x} \right)$, 展开为
$\begin{split}& \left( {\varepsilon _{n,n'}^0 + 2\beta } \right){\psi _{ \uparrow \uparrow }} \\=\;& \alpha \left( {\frac{{{K_n}}}{2} - k_n^{n'}} \right){\psi _{ \uparrow \downarrow }} + \alpha \left( {\frac{{{K_n}}}{2} + k_n^{n'}} \right){\psi _{ \downarrow \uparrow }},\end{split}\tag{8a}$
$\begin{split}\;&\left( {\varepsilon _{n,n'}^0 - 2\beta } \right){\psi _{ \downarrow \downarrow }}\\ =\;& \alpha \left( {\frac{{{K_n}}}{2} + k_n^{n'}} \right){\psi _{ \uparrow \downarrow }} + \alpha \left( {\frac{{{K_n}}}{2} - k_n^{n'}} \right){\psi _{ \downarrow \uparrow }},\end{split}\tag{8b}$
$\begin{split}\varepsilon _{n,n'}^0{\psi _{ \uparrow \downarrow }} =\;& \alpha \left( {\frac{{{K_n}}}{2} + k_n^{n'}} \right){\psi _{ \uparrow \uparrow }}+ \frac{g}{{2L}}\sum\limits_{n'} {\left( {{\psi _{ \uparrow \downarrow }} - {\psi _{ \downarrow \uparrow }}} \right)} \\&+ \alpha \left( {\frac{{{K_n}}}{2} - k_n^{n'}} \right){\psi _{ \downarrow \downarrow }},\\[-20pt]\end{split}\tag{8c}$
$\begin{split}\varepsilon _{n,n'}^0{\psi _{ \downarrow \uparrow }} =\;& \alpha \left( {\frac{{{K_n}}}{2} - k_n^{n'}} \right){\psi _{ \uparrow \uparrow }} - \frac{g}{{2L}}\sum\limits_{n'} {\left( {{\psi _{ \uparrow \downarrow }} - {\psi _{ \downarrow \uparrow }}} \right)} \\&+ \alpha \left( {\frac{{{K_n}}}{2} + k_n^{n'}} \right){\psi _{ \downarrow \downarrow }},\\[-20pt]\end{split}\tag{8d}$
其中$\varepsilon _{n, n'}^0 = {E_n} - \dfrac{{{\hbar ^2}}}{{2 m}}\left[ {\dfrac{{K_n^2}}{2} + 2{{\left( {k_n^{n'}} \right)}^2}} \right]$.
引入${\psi _{\rm{s}}} = \dfrac{1}{{\sqrt 2 }}\left( {{\psi _{ \uparrow \downarrow }} - {\psi _{ \downarrow \uparrow }}} \right)$${\psi _{\rm{a}}} = \dfrac{1}{{\sqrt 2 }}( {\psi _{ \uparrow \downarrow }} + $$ {\psi _{ \downarrow \uparrow }})$, 综合方程(8a)—(8d), 可以得到${\psi _{\rm{s}}}$${\psi _{\rm{a}}}$的表达式,
$\begin{split}&{\psi _{\rm{s}}} =\\&\frac{{g\displaystyle\sum\nolimits_{n'} {{\psi _{\rm{s}}}} }}{{L\left( {\varepsilon _{n,n'}^0 - 2\varepsilon _{n,n'}^ + {{\left( {\alpha k_n^{n'}} \right)}^2} - \dfrac{{{{\left( {{\alpha ^2}{K_n}k_n^{n'}\varepsilon _{n,n'}^ - } \right)}^2}}}{{\varepsilon _{n,n'}^0 - \dfrac{1}{2}\varepsilon _{n,n'}^ + {{\left( {\alpha {K_n}} \right)}^2}}}} \right)}},\end{split}\tag{9a}$
${\psi _{\rm{a}}} = - \frac{{{\alpha ^2}{K_n}k_n^{n'}\varepsilon _{n,n'}^ - {\psi _{\rm{s}}}}}{{\varepsilon _{n,n'}^0 - \dfrac{1}{2}{{\left( {\alpha {K_n}} \right)}^2}\varepsilon _{n,n'}^ + }},\tag{9b}$
其中
$\varepsilon _{n,n'}^ \pm = \frac{1}{{\varepsilon _{n,n'}^0 - 2\beta }} \pm \frac{1}{{\varepsilon _{n,n'}^0 + 2\beta }}.$
${\psi _{\rm{s}}}$的方程式(9a)推得本征能量的自洽方程, 表示为
$\frac{L}{g} = \sum\limits_{n'} {\frac{1}{{\varepsilon _{n,n'}^0 \!-\! 2\varepsilon _{n,n'}^ + {{\left( {\alpha k_n^{n'}} \right)}^2} \!-\! \dfrac{{{{\left( {{\alpha ^2}{K_n}k_n^{n'}\varepsilon _{n,n'}^ - } \right)}^2}}}{{\varepsilon _{n,n'}^0 - \dfrac{1}{2}\varepsilon _{n,n'}^ + {{\left( {\alpha {K_n}} \right)}^2}}}}}} .$
通过数值求解(11)式的自洽方程, 获得了两体费米模型精确的本征能量${E_n}$. 由方程(11)可得出, 本征能量关于$ \pm {K_n}$是简并的, 以下分析以${K_n} > 0$为代表. 将解出的本征能量代入方程(9a)和方程(9b)中, 获得精确的${\psi _{\rm{s}}}$${\psi _{\rm{a}}}$. 同时, ${\psi _{ \uparrow \downarrow }}$${\psi _{ \downarrow \uparrow }}$通过${\psi _{\rm{s}}}$${\psi _{\rm{a}}}$推出, 另外两个分量${\psi _{ \uparrow \uparrow }}$${\psi _{ \downarrow \downarrow }}$通过${\psi _{ \uparrow \downarrow }}$${\psi _{ \downarrow \uparrow }}$获得, 归一化后得到精确的波函数. 通过精确的本征能量和波函数, 实现了两体费米系统物理特性的精确研究.
针对不同的自旋轨道耦合、塞曼相互作用和接触相互作用系数, 本文分析了两体费米系统的基态能谱和动量分布. 取$\hbar = 1$, $m = 1$$L = 1$. 总动量${K_n}$是好量子数, 可以用n表征. 当$n = 0$$n = {\rm{1}}$, 总动量分别为${K_{\rm{0}}}{{ = 0}}$${K_{\rm{1}}} = 2{\text{π}}$.
首先, 考虑仅有接触相互作用的费米系统, 即当$\alpha {{ = 0}}$$\beta {{ = 0}}$时, 本征能量的自洽方程(11)式约化为
$ \frac{1}{g} = \sum\limits_{n'} {\frac{1}{{{E_n} - \dfrac{1}{2}\left[ {\dfrac{{K_n^2}}{2} + 2{{\left( {k_n^{n'}} \right)}^2}} \right]}}}. $
通过简化得到
$ \frac{1}{g} = \left\{ {\begin{aligned} &{{\rm{cot}}\left(\sqrt{{E}_{n}-\frac{{K}_{n}^{2}}{4}}\frac{1}{2}\right)\Bigg/\left(2\sqrt{{E}_{n}-\frac{{K}_{n}^{2}}{4}}\right), }\\ &\qquad\qquad\qquad\qquad\qquad\qquad\qquad{n{\text{为偶数}},}\\ &{-{\rm{tan}}\left(\sqrt{{E}_{n}-\frac{{K}_{n}^{2}}{4}}\frac{1}{2}\right)\Bigg/\left(2\sqrt{{E}_{n}-\frac{{K}_{n}^{2}}{4}}\right),}\\ &\qquad\qquad\qquad\qquad\qquad\qquad\qquad{n{\text{为奇数}},} \end{aligned}} \right. $
根据该式可得
$ \begin{split} &{k'=(2I{{ - }}1){\rm{\pi}} - 2{\rm{arctan}} \left( {{{k'} / {c'}}} \right),}\quad{n{\text{为偶数}},}\\ &{k'=2I{\rm{\pi }} - 2{\rm{arctan}} \left( {{{k'} / {c'}}} \right),}\;\;\qquad\quad{n{\text{为奇数}},} \end{split} $
其中$k' = \sqrt {{E_n} - {{K_n^2}}/{4}}$, $c' = g/2$. 当I为整数时, (13)式与两体杂核系统精确解[33]中的(12)式结果一致.
考虑仅具有自旋轨道耦合的费米相互作用系统. 即当$\beta = 0$, 本征能量的自洽方程(11)式约化为
$\begin{split}& {1}/{g} =\\ & \sum\limits_{n'} {\frac{1}{{{E_n} \!-\! \dfrac{1}{2}\Big[ \dfrac{{K_n^2}}{2} \!+\! 2{{( {k_n^{n'}})}^2} \Big] \!-\! \dfrac{{4{\alpha ^2}{{( {k_n^{n'}} )}^2}}}{{{E_n} \!-\! \dfrac{1}{2}\Big[{\dfrac{{K_n^2}}{2} + 2{{( {k_n^{n'}} )}^2}} \Big]}}}}} . \end{split}$
由(14)式可得, 当$g \to {\rm{0}}$, 本征能量${E_n} \to {{K_n^2}}/{4} + $$ {\left( {k_n^{n'}} \right)^2} \pm 2\alpha k_n^{n'}$, 约化为单粒子的本征能量. 给定$g = 1$, 图1给出了当$n = 0$时, 在总动量${K_0} = 0$空间, 最低能量${E_{\rm{0}}}$随着自旋轨道耦合系数$\alpha $变化的图像. 随着$\alpha $的增大, 本征能量逐渐变小. 图1中的插图给出了能量${E_{\rm{0}}} + {\alpha ^2}$随着$\alpha $增加而周期性的变化. 在无塞曼相互作用情况下, 研究发现自旋轨道耦合项通过幺正变换, 能够变换到周期性边界条件中, 展现出周期性地依赖于自旋轨道耦合系数[37].
图 1 最低能量$ {E_{\rm{0}}}$随着自旋轨道耦合参量$ \alpha$的变化. 插图为$ {E_0} + {\alpha ^2}$$ \alpha$的变化, $ n = 0$, $ g = 1$, $ \beta = 0$
Figure1. The lowest energy ${E_{\rm{0}}}$ versus spin-orbit coupling parameter $\alpha $. Insert is ${E_0} + {\alpha ^2}$ versus $\alpha $. $n = 0$, $g = 1$, $\beta = 0$.

考虑具有自旋轨道耦合和塞曼相互作用的费米系统, 对于不同的${K_n}$, 本征能量${E_n}$通过数值求解自洽方程(11)式获得. 特别地, 对于${K_{\rm{0}}} = 0$, 本征能量的自洽方程(11)式可以约化为
$ \begin{split} \frac{1}{g} =\;& \sum\limits_{n'} {\frac{1}{{2\left( {{\beta ^2} + {\alpha ^2}k_n^{n'2}} \right)}}} \left[ \frac{{2{\beta ^2}}}{{{E_n} - k_n^{n'2}}}\right. \\ &+ \frac{{{\alpha ^2}k_n^{n'2}}}{{{E_n} - k_n^{n'2} - 2\sqrt {{\beta ^2} + {\alpha ^2}k_n^{n'2}} }} \\ &\left.+ \frac{{{\alpha ^2}k_n^{n'2}}}{{{E_n} - k_n^{n'2} + 2\sqrt {{\beta ^2} + {\alpha ^2}k_n^{n'2}} }} \right],\end{split} $
从该式可得, 当相互作用$g \to 0$时, 本征能量趋于$k_n^{n'2}$$k_n^{n'2} \pm 2\sqrt {{\beta ^2} + {\alpha ^2}k_n^{n'2}} $.
图2给出了在较弱相互作用$g = 1$$\alpha = 1.2{\text{π}}$条件下, 随着塞曼相互作用系数$\beta $变化, 在不同总动量${K_n}$空间中最低能量${E_n}$变化的图像. 如图2所示, 随着塞曼相互作用$\beta $增大, 能量${E_n}$逐渐减小, 在两个总动量空间的最低能级间发生了能级交叉. 当$\beta $较小时, ${E_{\rm{0}}}< {E_{\rm{1}}}$, 如取$\beta = {\text{π}}$时, 两体费米系统基态能量为${E_{\rm{0}}}$, 基态的总动量为${K_{\rm{0}}}$. 当$\beta $足够大时, 相反地, ${E_{\rm{1}}} < {E_{\rm{0}}}$, 如取$\beta = 2{\text{π}}$时, ${E_{\rm{1}}} < {E_{\rm{0}}}$, 两体费米系统基态能量为${E_{\rm{1}}}$, 基态的总动量为${K_{\rm{1}}}$. 综上所述, 随着$\beta $的增加, 两体费米系统的基态从${K_{\rm{0}}}$转移至${K_{\rm{1}}}$空间. 通过增加塞曼相互作用强度, 系统基态从质心静止转移到了质心运动的状态.
图 2 在弱相互作用条件下最低能量$ {E_{\rm{0}}}$$ {E_1}$随塞曼系数$ \beta $的变化, $ g = 1$, $ \alpha = 1.2{\text{π}}$
Figure2. The lowest energies ${E_{\rm{0}}}$ and ${E_1}$ versus Zeeman interaction parameter $\beta $ in the weak interaction condition. $g = 1$, $\alpha = 1.2{\text{π}}$.

图3给出了在任意相互作用强度g下本征能量${E_n}$的变化. 横坐标表示费米子间的相互作用系数, 从负到正表示从吸引相互作用到排斥相互作用. 给定自旋轨道耦合系数$\alpha = 1.2{\text{π}}$, 图3(a)图3(c)对比了不同塞曼相互作用强度$\beta $时, 能量${E_n}$随粒子间相互作用系数g的变化. 图3(a)呈现了在无塞曼相互作用($\beta = 0$)情况下, 能量${E_{\rm{0}}}$总是小于${E_{\rm{1}}}$, 即基态的总动量始终为零. 考虑塞曼相互作用时, 图3(b)图3(c)均呈现了能级交叉现象. 在图3(b)中, 当塞曼相互作用系数$\beta = {\text{π}}$时, 在强吸引相互作用下基态能量为${E_{\rm{0}}}$, 强排斥下为${E_{\rm{1}}}$, 而在图2中, 相同$\beta $时在弱排斥作用下基态能量为${E_{\rm{0}}}$, 与强排斥下不同. 类似地, 图3(c)呈现了$\beta = 2{\text{π}}$时本征能量的图像, 强吸引下基态能量为${E_{\rm{0}}}$, 强排斥下为${E_{\rm{1}}}$, 对比图2中在弱排斥作用下基态能量仍为${E_{\rm{1}}}$. 综上所述, 具有塞曼相互作用时, 两体费米系统能够在强相互作用下呈现和弱相互作用下不同的基态特性. 随着相互作用从吸引到排斥转变, 具有塞曼相互作用的费米系统基态能量从${E_{\rm{0}}}$转变到${E_{\rm{1}}}$, 基态总动量从零变到有限值.
图 3 能量${E_0}$${E_1}$随着接触相互作用g的变化, $\alpha = 1.2{\text{π}}$ (a)$\beta = 0$; (b)$\beta = {\text{π}}$; (c) $\beta = 2{\text{π}}$
Figure3. The energies ${E_{\rm{0}}}$ and ${E_1}$ versus contact interaction parameter g, $\alpha = 1.2{\text{π}}$: (a)$\beta = 0$; (b)$\beta = {\text{π}}$; (c)$\beta = 2{\text{π}}$.

为了直观理解基态的性质, 我们讨论了无接触相互作用($g = 0$)时, 单费米子的本征能量和本征态. 单粒子哈密顿量为${H_0} = - \dfrac{{\rm{1}}}{{\rm{2}}}\dfrac{{{\partial ^2}}}{{\partial {x^2}}} + \alpha {\sigma ^x}\times $$ \left( { - {\rm{i}}\dfrac{\partial }{{\partial x}}} \right) + \beta {\sigma ^z}$. 不考虑自旋轨道耦合和塞曼相互作用($\alpha = 0$, $\beta = 0$)时, 单费米子能谱不依赖于自旋, 基态动量${p_{\rm{0}}} = 0$. 两个费米子分别占据不同自旋态的最低能级, 此时总动量为零.
仅有自旋轨道耦合($\beta = 0$)系统, 单粒子本征能量为$e_{{p_{\rm{0}}}}^ \rightleftarrows = {{p_{\rm{0}}^{\rm{2}}}}/{2} \pm \alpha {p_{\rm{0}}}$, 自旋轨道耦合项是动量${p_{\rm{0}}}$的线性项, $ \leftarrow $$ \to $代表${\sigma _x}$的自旋态. 在自由空间时, 单粒子基态从${p_{\rm{0}}} = 0$转变到有限动量值$ {p}_{{\rm{0}}, \pm }=\pm \alpha $, 符号取决于自旋分量. 同时改变动量和自旋 (${p_{\rm{0}}} \to - {p_{\rm{0}}}$${\sigma _x} \to - {\sigma _x}$), 哈密顿量保持不变, 系统具有Z2对称性, 此时系统基态为两个简并态, 动量值相等但符号相反. 因此, 当自旋轨道耦合占主导时, 两费米子占据两个简并基态, 总动量为零. 在周期性边界条件下, 动量为一系列离散值${p_{\rm{0}}} = 2{\text{π}}{n_0}$, ${n_{\rm{0}}}$为整数, 此时基态动量取与$ \pm \alpha $最相近的值, 即$\left| {{p_{\rm{0}}} \pm \alpha } \right|$最小. 当$- {\text{π}} + 2{\text{π}}{n_{\rm{0}}} < \alpha < {\text{π}} $$ + 2{\text{π}} {n_{\rm{0}}}$, 基态相应的动量为$ {p}_{0, \pm }=\pm 2{\text{π}}{n}_{\rm{0}}$. 当$\alpha = 0.2{\text{π}}$, 则有${n_0} = {\rm{0}}$, 基态为${p_{\rm{0}}} = {\rm{0}}$. 当$\alpha = 1.2{\text{π}}$, 则有${n_0} = 1$, 基态为$ {p}_{0, \pm }=\pm 2{\text{π}}$.
最后, 考虑具有塞曼相互作用和自旋轨道耦合的费米子系统. 仅考虑塞曼相互作用($\alpha = 0$)时, 自旋向上和自旋向下的单粒子能谱发生劈裂, 产生能隙, 且能隙的大小与塞曼相互作用的系数$\beta $有关. 当$\beta $为正时, 自旋向下的费米子能级降低, 且最低能级仍处于${p_{\rm{0}}} = {\rm{0}}$, 而自旋向上的能级升高, 多体费米系统将先填充自旋向下的费米子能级, 再填充自旋向上的能级, 因此多体费米系统基态的总动量由所填充的能级决定. 考虑两费米子基态, 一个费米子将填充最低能级${p_{\rm{0}}} = {\rm{0}}$, 另一个填充${p_{\rm{0}}} \ne {\rm{0}}$的态, 此时总动量为有限值. 当塞曼效应占主导时, 两费米子的总动量为有限值. 考虑塞曼相互作用和自旋轨道耦合作用, 单粒子本征能量为$e_{{p_0}}^ + = \dfrac{{p_{\rm{0}}^{\rm{2}}}}{2} + \sqrt {{\alpha ^{\rm{2}}}p_{\rm{0}}^{\rm{2}} + {\beta ^{\rm{2}}}}$, $e_{{p_{\rm{0}}}}^ - = \dfrac{{p_{\rm{0}}^{\rm{2}}}}{2} - \sqrt {{\alpha ^{\rm{2}}}p_{\rm{0}}^{\rm{2}} + {\beta ^{\rm{2}}}}$. 在自由空间时, 当${\alpha ^2} > \beta $时, 自旋轨道耦合效应占主导作用时, 单粒子基态是简并的, 相应的动量为${p}_{0, \pm }\!=\! \pm \alpha \sqrt{1-\left( {{\beta }^{2}}/{{\alpha }^{4}}\right)}$, 能量为$e_{{p_{0, \pm }}}^ - \!=\! - (1/2)[ {\alpha ^{\rm{2}}} + $$ {{\left( { {\beta }/{\alpha }} \right)}^{\rm{2}}} ]$. 随着$\beta $的增大, 两个极小值${p_{{\rm{0}}, \pm }}$逐渐靠近, 当$\beta = {\alpha ^2}$时, 这两个极小值融合为单个最小值${p_{\rm{0}}} = 0$. 当$\beta > {\alpha ^2}$时, 塞曼相互作用占主导作用, 基态${p_{\rm{0}}} = 0$, 相应的基态能量为$ e_{{0}}^ - = - \beta $. 在周期性边界条件下, 类似地, 动量仅能取离散值${p_{\rm{0}}} = $$ 2{\text{π}}{n_0}$. 当自旋轨道耦合占主导作用时, 需满足$e_{{n_0} \ne {\rm{0}}}^ - < e_{\rm{0}}^ - $条件, 两体费米子能够处于简并基态, 此时基态对应的动量为$ {p}_{0, \pm }= $$ \pm 2{\text{π}}{n}_{\rm{0}}$. 当塞曼相互作用占主导时, 则有$e_{\rm{0}}^ - < e_{{n_0} \ne {\rm{0}}}^ - $, 基态动量为${p_0} = 0$, 次低能级的动量${p_0}' \ne 0$. 例如, 当${n_{\rm{0}}} = 1$时, 依据$e_{\rm{0}}^ - < e_{2{\text{π}}}^ - $推出$\beta > {\alpha ^2} - {{\text{π}}^2}$, 令${\beta _{\rm{c}}} = {\alpha ^2} - {{\text{π}}^2}$. 当$\beta < {\beta _{\rm{c}}}$, 自旋轨道耦合为主导作用, 两体费米子占据两个简并基态, 总动量为零. 当$\beta > {\beta _{\rm{c}}}$, 塞曼相互作用为主导作用, 两体费米子依次填充最低能级、零动量和有限动量的能级, 总动量为有限值.
图4给出了不同自旋轨道耦合和塞曼相互作用参量时的单粒子能级结构. 图4(a)是自旋轨道耦合和塞曼相互作用参量较小($\alpha = 0.2{\text{π}}$$\beta = $$ 0.1{\text{π}}$)时的能级结构图. 可以看出, 能谱上下分支$e_{{p_0}}^ \pm $最低能级均为${p_{\rm{0}}} = 0$, 两体无相互作用费米子占据这两个态. 考虑$\alpha = 1.2{\text{π}}$, 则有${n_{\rm{0}}} = 1$, 此时${\beta _{\rm{c}}} = 0.{\rm{44}}{{\text{π}}^{\rm{2}}} \approx $$ {\rm{4}}.{\rm{34}}$. 图4(b)是自旋轨道耦合和塞曼相互作用参量较大时($\alpha = {\rm{1}}.2{\text{π}}$$\beta = {\text{π}} < {\beta _{\rm{c}}}$)的能级图, 两体无相互作用费米子占据$e_{{p_0}}^ - $的两个简并基态, 两个动量取${p_{{\rm{0}} \pm }} = $$ \pm 2{\text{π}}$, 对应的总动量仍为零. 图4(c)是塞曼相互作用更大时($\alpha = {\rm{1}}.2{\text{π}}$$\beta = 2{\text{π}} > {\beta _{\rm{c}}}$)的能级图, 两体无相互作用费米子占据$e_{{p_0}}^ - $的两个最低能级, 基态动量为${p_{\rm{0}}} = 0$, 次之为${p_{\rm{0}}}' \ne {\rm{0}}$, 对应的总动量为有限值. 当塞曼相互作用增大时, 总动量从零转变到有限值, 单粒子分析很好地解释了图2中弱排斥相互作用时基态转变特性.
图 4 单粒子能级e和e+ (a) $ \alpha = 0.2{\text{π}}$, $ \beta = 0.1{\text{π}}$; (b) $ \alpha = 1.2{\text{π}}$, $ \beta = {\text{π}}$; (c) $ \alpha = 1.2{\text{π}}$, $ \beta = 2{\text{π}}$
Figure4. The single particle eigenenergies with two branches e and e+: (a) $\alpha = 0.2{\text{π}}$, $\beta = 0.1{\text{π}}$; (b) $\alpha = 1.2{\text{π}}$, $\beta = {\text{π}}$; (c) $\alpha = 1.2{\text{π}}$, $\beta = 2{\text{π}}$.

两体费米系统的动量分布可以表示为
$n\left( {{p_1}} \right) = \int {{\psi ^\dagger }\left( {{p_1},{p_2}} \right)\psi \left( {{p_1},{p_2}} \right){\rm{d}}{p_2}} ,\tag{15a}$
$n\left( {{p_2}} \right) = \int {{\psi ^\dagger }\left( {{p_1},{p_2}} \right)\psi \left( {{p_1},{p_2}} \right){\rm{d}}{p_1}} ,\tag{15b}$
其中, $\psi \left( {{p_1}, {p_2}} \right)$$\psi \left( {{x_1}, {x_2}} \right)$的傅里叶变换而来. 图5给出了在较弱相互作用时对于不同自旋轨道耦合$\alpha $和塞曼相互作用$\beta $下基态的动量分布, 且自旋轨道耦合与塞曼相互作用系数取值与图4中一致. 因为$n\left( {{p_1}} \right) = n\left( {{p_2}} \right)$, 这里以$n\left( {{p_1}} \right)$作为代表. 当$\alpha = 0.2{\text{π}}$$\beta = 0.1{\text{π}}$时, 粒子动量分布集中于$p = 0$. 当$\alpha = $$ 1.2{\text{π}}$$\beta = {\text{π}}$时, 动量分布出现两个峰, 为$p = \pm 2{\text{π}}$. 当$\alpha = 1.2{\text{π}}$$\beta = {\rm{2\pi }}$时, 动量分布出现两个峰, 为$p = {\rm{0}}$$p = 2{\text{π}}$. 图5是在弱排斥作用情况下的基态动量分布, 与图4中无相互作用情况下, 两费米子动量分布的分析一致. 通过单粒子能谱分析和动量分布研究, 直观地理解了随着塞曼相互作用变化, 费米子系统基态的转变过程, 为理解多体系统的物理现象提供了基础.
图 5 对于不同的$ \alpha $$ \beta $, 基态的动量分布, $ g = 1$
Figure5. The momentum distributions of ground states for different $\alpha $ and $\beta $, $g = 1$.

本文解析求解了环阱中具有自旋轨道耦合、塞曼作用和接触相互作用的两体费米系统. 研究发现: 首先, 仅考虑自旋轨道耦合作用时, 随着自旋轨道耦合系数增加, 本征能量将降低, 主要依赖于与自旋轨道耦合相关的常数项. 同时, 不考虑上述常数项时, 本征能量周期性地依赖于自旋轨道耦合系数. 进一步地, 考虑自旋轨道耦合和塞曼场作用, 随着塞曼相互作用增大, 不同总动量空间的最低能量均有降低, 并且在这些最低能级间发生了能级交叉, 基态从零转变到有限值动量空间. 最后, 随着两体接触相互作用从吸引到排斥的变化, 无塞曼作用时基态总动量始终为零, 考虑塞曼作用时, 基态总动量从零转变为有限值. 通过单粒子能谱分析, 发现在塞曼相互作用为主导时, 塞曼能级劈裂引起了基态的这种转变. 两体费米系统的解析解展现了丰富的物理现象, 是研究多体系统的基础. 在本研究基础上, 未来可以进一步研究多体相互作用费米系统、杂核系统、光晶格势阱系统、多旋量费米系统的物理图像.
相关话题/系统 空间 物理 运动 动量

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于铌酸锂双折射晶体的皮秒拍瓦激光系统光谱整形
    摘要:为补偿皮秒拍瓦激光系统中钕玻璃宽带放大引起的增益窄化,提出了一种基于铌酸锂双折射晶体的高能量光谱整形方法.在相同强度调制下,对比了BBO、铌酸锂和石英3种晶体,针对1053nm激光,选用了高双折射率、大口径且不易潮解的铌酸锂作为整形晶体.理论分析了晶体厚度、倾斜角、面内旋转角对强度调制的影响, ...
    本站小编 Free考研考试 2021-12-29
  • 自旋轨道耦合Su-Schrieffer-Heeger原子链系统的电子输运特性
    摘要:在Su-Schrieffer-Heeger(SSH)原子链中,电子在胞内和胞间的跳跃依赖于其自旋时,即SSH原子链存在自旋轨道耦合作用时,存在不同缠绕数的非平庸拓扑边缘态.如何探测自旋轨道耦合SSH原子链不同缠绕数的边缘态是一个重要问题.本文在紧束缚近似下研究了自旋轨道耦合SSH原子链的非平庸 ...
    本站小编 Free考研考试 2021-12-29
  • 基于对偶单元法的三维集成微系统电热耦合分析
    摘要:随着三维集成微系统集成度和功率密度的提高,同时考察电设计与热管理的多场耦合分析势在必行.本文面向三维集成微处理器系统,通过改进的对偶单元法(dualcellmethod,DCM)实现了系统的快速电热分析.该方法通过引入泄漏功率、材料系数随温度的耦合,相比于传统有限元法在更新以及组装本构矩阵上有 ...
    本站小编 Free考研考试 2021-12-29
  • 单壁碳纳米管受限空间内水的分布
    摘要:碳纳米管管腔作为分子物质的纳米通道,其储存或输送水的能力具有重要研究价值.为了研究碳纳米管管腔受限空间对水分子团簇结构和分布的影响,本文采用分子动力学方法探究了管径、手性和温度对单壁碳纳米管管腔内水的结构和分布的影响.结果表明:在常温下,管径尺寸范围为1.018—1.253nm的单壁碳纳米管管 ...
    本站小编 Free考研考试 2021-12-29
  • 铝基薄膜忆阻器作为感觉神经系统的习惯化特性
    摘要:感觉神经系统可在外界刺激与生物体反应之间建立联系.感觉神经系统中的最小单位神经元可直接将外界刺激传递至中枢神经,再由中枢神经通过控制和调节生物体对外界刺激作出反应.神经突触连接了相邻神经元进行脉冲信息传递功能.习惯化是神经突触在信息传递中过滤外界无关信息时的一个基本特性,可以让感觉神经系统更快 ...
    本站小编 Free考研考试 2021-12-29
  • 考虑空间电荷层效应的氧离子导体电解质内载流子传输特性
    摘要:晶界或异质界面诱发的空间电荷层(spacechargelayer,SCL)效应,被认为是氧离子导体电解质内界面附近区域载流子传输特性显著区别于体相区域的关键原因之一.现有研究多采用Poisson-Boltzmann(PB)方程预测SCL效应的影响规律,但其基于载流子电化学平衡假设,无法用于载流 ...
    本站小编 Free考研考试 2021-12-29
  • Duffing系统的主-超谐联合共振
    摘要:非线性动力系统极易发生共振,在多频激励下可能发生联合共振或组合共振,目前关于非线性系统的主-超谐联合共振的研究少见报道.本文以Duffing系统为对象,研究系统在主-超谐联合共振时的周期运动和通往混沌的道路.应用多尺度法得到系统的近似解析解,并利用数值方法对解析解进行验证,结果吻合良好.基于L ...
    本站小编 Free考研考试 2021-12-29
  • 四稳系统的双重随机共振特性
    摘要:提出了一类8次势函数并讨论了其分岔特性,得到由左、右2个小尺度双稳势和中间势垒构成的对称四稳系统.建立了在周期力和随机力共同作用下四稳系统输出响应的近似解析表达式,并从能量角度引入功这一过程量来刻画大、小不同尺度双稳势之间的作功能力,发现四稳势中存在着双重随机共振现象.理论分析与数值仿真结果表 ...
    本站小编 Free考研考试 2021-12-29
  • 一种用于Z箍缩实验的软X射线成像系统
    摘要:基于塑料闪烁体转换和光学条纹相机的方法建立了一套用于Z箍缩实验中的软X射线条纹图像诊断系统,解决了以往实验中使用的X射线条纹相机易被电磁环境干扰以及相机电极部件易被实验产生的高速粒子损伤的问题.诊断系统的光谱响应范围主要集中在0.2—10keV,系统的空间分辨率经过理论评估小于120μm,通过 ...
    本站小编 Free考研考试 2021-12-29
  • 降雪对地表附近自由空间量子信道的影响及参数仿真
    摘要:量子通信具有覆盖面广、安全保密的优势,是当前通信领域国内外的研究热点.在自由空间量子通信过程中,光量子信号需要在地表上空一定高度进行传输,因此各种环境因素,例如降雪、沙尘暴、降雨、雾霾、浮尘等,不可避免地会影响量子通信性能.然而,迄今为止,降雪对地表附近自由空间量子信道影响的研究尚未展开.为此 ...
    本站小编 Free考研考试 2021-12-29