搜索

x

留言板

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

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

由超辐射引起的迁移率边和重返局域化

吴瑾 陆展鹏 徐志浩 郭利平

引用本文:
Citation:

由超辐射引起的迁移率边和重返局域化

吴瑾, 陆展鹏, 徐志浩, 郭利平

Mobility edges and reentrant localization induced by superradiance

Wu Jin, Lu Zhan-Peng, Xu Zhi-Hao, Guo Li-Ping
PDF
HTML
导出引用
  • 研究了束缚在梯子形光晶格中的中性原子与高精度腔耦合的系统, 发现由超辐射引起的准周期调制可以导致迁移率边和重返局域化现象的出现. 在平均场近似下, 超辐射现象可以引起两种不同频率的准周期调制, 它们可以由腔场和泵浦场有效调节. 在本文的观测范围内, 当高频调制强度小于某一临界值时, 系统随着低频调制强度的增加经历一次局域转变. 通过数值求解分形维度、密度分布、平均参与率、平均逆参与率以及标度分析, 证明了局域转变历经的临界相区存在迁移率边; 当高频调制强度大于临界值时, 随着低频调制强度的增强, 系统依次经历完全扩展相-临界相-完全局域相-临界相-完全局域相, 这是一个典型的重返局域化现象. 最后给出了局域化相图. 该研究结果为超辐射相变和重返局域化现象的研究建立了联系, 也为重返局域化的研究搭建了新的平台.
    We study a Bose-Einstein condensate trapped by a ladder lattice in a high-fitness cavity. The ladder lattice is loaded in the $x\text-y$ plane and the cavity is along the x direction. A pump laser shines on atoms from the z direction. Under the mean-field approximation, we consider the emergence of the quasi-periodic potentials induced by superradiance in the ladder lattice, which is described by $\hat{H}_{\text{MF}}=\hat{H}_{\text{Lad}}+\hat{V}_{\text{eff}}$ with the effective potential $\hat{V}_{\text{eff}}(\alpha)={\displaystyle \sum\nolimits_{i = 1}^{N}}\displaystyle \sum\nolimits_{\sigma = 1,2}\left[\lambda_{\rm{D}}\cos({2\pi\beta i})+U_{\rm{D}}\cos^{2}({2\pi\beta i})\right]\hat{c}^{†}_{i,\sigma}\hat{c}_{i,\sigma}$. We find that the quasi-periodic potential can induce the reentrant localization transition and the regime with mobility edges. In the smaller $U_{\rm{D}}$ case, the system exhibits a localization transition. The transition is associated with an intermediate regime with mobility edges. When $U_{\rm{D}}$ goes beyond a critical value $U_{\rm{D}}^{(\rm c)}$, with the increase of $\lambda_{\rm{D}}$, the system undergoes a reentrant localization transition. This indicates that after the first transition, some of the localized eigenstates change back to the extended ones for a range of $\lambda_{\rm{D}}$. For a larger $\lambda_{\rm{D}}$, the system experiences the second localization transition, then all states become localized again. Finally, the local phase diagram of the system is also discussed. This work builds a bridge between the reentrant localization and the superradiance, and it provides a new perspective for the reentrant localization.
      通信作者: 徐志浩, xuzhihao@sxu.edu.cn ; 郭利平, guolp@sxu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11604188, 12147215, 11904216)、山西省基础研究计划(批准号: 20210302123442)、北京凝聚态物理国家研究中心开放课题、山西省高等学校科技创新项目(批准号: 2019L0097) 和山西省“1331工程”重点学科建设计划资助的课题
      Corresponding author: Xu Zhi-Hao, xuzhihao@sxu.edu.cn ; Guo Li-Ping, guolp@sxu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11604188, 12147215, 11904216), the Fundamental Research Program of Shanxi Province, China (Grant No. 20210302123442), the Open Project of Beijing National Laboratory for Condensed Matter Physics, China, the Scinence and Technology Innovation Project of Higher Education Institutions in Shanxi Province, China (Grant No. 2019L0097), and the Fund for Key Subjects of Shanxi “1331 Project”, China
    [1]

    Marklund M, Shukla P K 2006 Rev. Mod. Phys. 78 591Google Scholar

    [2]

    Chin C, Grimm R, Julienne P, Tiesinga E 2010 Rev. Mod. Phys. 82 1225Google Scholar

    [3]

    Mekhov I B, Maschler C, Ritsch H 2007 Nat. Phys. 3 319Google Scholar

    [4]

    Ritsch H, Domokos P, Brennecke F, Esslinger T 2013 Rev. Mod. Phys. 85 553Google Scholar

    [5]

    Baumann K, Guerlin C, Brennecke F, Esslinger T 2010 Nature 464 1301Google Scholar

    [6]

    Maschler C, Mekhov I B, Ritsch H 2008 Eur. Phys. J. D 46 545Google Scholar

    [7]

    Klinder J, Keβler H, Bakhtiari M R, Thorwart M, Hemmerich A 2015 Phys. Rev. Lett. 115 230403Google Scholar

    [8]

    Landig R, Hruby L, Dogra N, Landini M, Mottl R, Donner T, Esslinger T 2016 Nature 532 476Google Scholar

    [9]

    Zheng W, Cooper N R 2018 Phys. Rev. A 97 021601(RGoogle Scholar

    [10]

    Zhou L, Pu H, Zhang K Y, Zhao X D, Zhang W P 2011 Phys. Rev. A 84 043606Google Scholar

    [11]

    Habibian H, Winter A, Paganelli S, Rieger H, Morigi G 2013 Phys. Rev. Lett. 110 075304Google Scholar

    [12]

    Anderson P W 1958 Phys. Rev. 109 1492Google Scholar

    [13]

    Aubry S, André G 1980 Ann. Isr.: Phys. Soc. 3 18

    [14]

    Roati G, D’Errico C, Fallani L, Fattori M, Fort C, Zaccanti M, Modugno G, Modugno M, Inguscio M 2008 Nature 453 895Google Scholar

    [15]

    Sarma S D, Kobayashi A, Prange R E 1986 Phys. Rev. Lett. 56 1280Google Scholar

    [16]

    Sarma S D, Song H, Xie X C 1988 Phys. Rev. Lett. 61 2144Google Scholar

    [17]

    Biddle J, Sarma S D 2010 Phys. Rev. Lett. 104 070601Google Scholar

    [18]

    Ganeshan S, Pixley J H, Sarma S D 2015 Phys. Rev. Lett. 114 146601Google Scholar

    [19]

    Lanini Y, Bromberg Y, Christodoulides D N, Silberberg Y 2010 Phys. Rev. Lett. 105 163905Google Scholar

    [20]

    Xu Z H, Xia X, Chen S 2021 arXiv: 2109.02072 v1[cond-mat.dis-nn]

    [21]

    Wang Y C, Xia X, Wang Y J, Zheng Z H, Liu X J 2021 Phys. Rev. B 103 174205Google Scholar

    [22]

    徐志浩, 皇甫宏丽, 张云波 2019 物理学报 68 087201Google Scholar

    Xu Z H, Huangfu H L, Zhang Y B 2019 Acta Phys. Sin. 68 087201Google Scholar

    [23]

    Xu Z H, Xia X, Chen S 2022 Sci. China-Phys. Mech. Astron. 65 227211Google Scholar

    [24]

    Roy S, Mishra T, Tanatar B, Basu S 2021 Phys. Rev. Lett 126 106803Google Scholar

    [25]

    Jiang X P, Qiao Y, Cao J P 2021 Chin. Phys. B 30 097202Google Scholar

    [26]

    Padhan A, Giri M K, Mondal S, Mishra T 2021 arXiv: 2109.09621 v1 [cond-mat.quant-gas]

    [27]

    Zhuang W F, Geng B, Lou H G, Guo G C, Gong M 2021 Phys. Rev. A 104 053308Google Scholar

    [28]

    Keeling J, Bhaseen M J, Simons B D 2014 Phys. Rev. Lett 112 143002Google Scholar

    [29]

    Piazza F, Strack P 2014 Phys. Rev. Lett. 112 143003Google Scholar

    [30]

    Chen Y, Yu Z H, Zhai H 2014 Phys. Rev. Lett. 112 143004Google Scholar

    [31]

    Chen Y, Zhai H, Yu Z H 2015 Phys. Rev. A 91 021602(RGoogle Scholar

    [32]

    Chen Y, Yu Z H, Zhai H 2016 Phys. Rev. A 93 041601(RGoogle Scholar

    [33]

    Cai M L, Liu Z D, Zhao W D, Wu Y K, Mei Q X, Jiang Y, He L, Zhang X, Zhou Z C, Duan L M 2021 Nat. Commun. 12 1126Google Scholar

    [34]

    Xu Z H, Chen S 2021 Phys. Rev. A 103 043325Google Scholar

  • 图 1  实验装置示意图. 在高精度腔中, 中性原子被x-y平面的梯子形光晶格俘获, 并沿着z方向加入驱动泵浦场. 泵浦场的频率为$\omega_{\rm{p}}$远失谐于原子能级跃迁频率$\omega_{\rm{a}}$, 但接近于腔场频率$\omega_{\rm{c}}$

    Fig. 1.  A schematic diagram of experimental setup. In a high-finesses optical cavity, spinless atoms are trapped by a ladder lattice in x-y plane. The atoms are driven by a pump laser beam along z direction. The frequency $ \omega_{\rm{p}} $ of the pump laser is far detuned from the atomic transition line $ \omega_{\rm{a}} $ but close to the cavity-mode frequency $ \omega_{\rm{c}} $

    图 2  (a) 腔场耦合强度U取不同的值时, 光场$ |\alpha| $随着耦合强度λ变化的图像; (b) 临界耦合强度$ \lambda_{\rm{c}} $随着U的变化情况. 这里, $ L = 1974 $, $ K = 0.8 $, $ \beta = 610/987 $, $ \varDelta_{{\rm{c}}} = -2 $$ \kappa = 1.2 $

    Fig. 2.  (a) The cavity field $ |\alpha| $ as a function of the pumping strength λ for different U; (b) the critical pumping strength $ \lambda_{\rm{c }}$ as a function of U. Here $ L = 1974 $, $ K = 0.8 $, $ \beta = 610/987 $, $ \varDelta_{{\rm{c}}} = -2 $ and $ \kappa = 1.2 $

    图 3  (a) 当$ U_{\rm{D}} = 1, K = 0.8, L = 3194 $时, $ \langle \text{IPR} \rangle $(黑色实线) 和$ \langle \text{NPR} \rangle $(红色虚线) 随$ \lambda_{\rm{D}} $变化的曲线, 灰色区域代表着具有迁移率边的临界区域; (b)$ \lambda_{\rm{D}} = 0.5 $, (c)$ \lambda_{\rm{D}} = 2 $, (d)$ \lambda_{\rm{D}} = 5 $时, $ \langle \text{NPR} \rangle $随着$ L^{-1} $的变化, 其余参数的取值是$ U_{\rm{D}} = 1 $, $ K = 0.8 $

    Fig. 3.  (a) $ \langle \text{IPR} \rangle $ (the black solid line) and $ \langle \text{NPR} \rangle $ (the red dashed line) as the functions of $ \lambda_{\rm{D}} $ for $U_{\rm{D}} = 1, K = 0.8, $$ L = 3194$. The grey region denotes the critical region with mobility edges; $ \langle \text{NPR} \rangle $ as a function of $ L^{-1} $ with $ U_{\rm{D}} = 1 $ and $ K = 0.8 $ for (b) $ \lambda_{\rm{D}} = 0.5 $, (c) $ \lambda_{\rm{D}} = 2 $, and (d) $ \lambda_{\rm{D }}= 5 $

    图 4  (a1) $ \lambda_{\rm{D}} = 0.5 $, (b1) $ \lambda_{\rm{D}} = 2 $和(c1) $ \lambda_{\rm{D}} = 5 $时第300个激发态的密度分布; (a2) $ \lambda_{\rm{D }}= 0.5 $, (b2) $ \lambda_{\rm{D}} = 2 $和(c2) $ \lambda_{\rm{D}} = 5 $时第 1100 个激发态的密度分布. 这里, $ K = 0.8 $, $ U_{\rm{D}} = 1 $$ L = 1220 $

    Fig. 4.  Density distributions of the 300th excited eigenstates for (a1) $ \lambda_{\rm{D}} = 0.5 $, (b1) $ \lambda_{\rm{D}} = 2 $, and (c1) $ \lambda_{\rm{D}} = 5 $; density distributions of the 1100th excited eigenstates for (a2) $ \lambda_{\rm{D}} = 0.5 $, (b2) $ \lambda_{\rm{D }}= 2 $, and (c2) $ \lambda_{\rm{D}} = 5 $. Here, $ K = 0.8 $, $ U_{\rm{D}} = 1 $, and $ L = 1220 $

    图 5  $ U_{\rm{D}} = 1 $, $ K = 0.8 $$ L = 3194 $时, 分形维度γ随着调制强度$ \lambda_{\rm{D }}$和能量本征值$ E_n $的变化. 其中, 黑色实线和黑色虚线分别对应了两个局域化转变点$ \lambda_{\rm{D}}^{({\rm{c}}1)} $$ \lambda_{\rm{D}}^{({\rm{c}}2)} $. 图中颜色代表着γ的大小

    Fig. 5.  γ of the eigenstates as a function of the energy spectrum and $ \lambda_{\rm{D}} $ for $ U_{\rm{D}} = 1 $, $ K = 0.8 $, and $ L = 3194 $. Here the black solid line and black dashed line denote two transport points, $ \lambda_{\rm{D}}^{({\rm{c}}1)} $ and $ \lambda_{\rm{D}}^{({\rm{c}}2)} $. The color code represents the values of γ

    图 6  $ U_{\rm{D}} = 2, K = 0.8, L = 3194 $时 (a)$ \langle \text{IPR} \rangle $(黑色实线)和$ \langle \text{NPR} \rangle $ (红色虚线)随着$ \lambda_{\rm{D}} $变化的曲线, 灰色区域表示具有迁移率边的临界相区; (b)分形维度γ随着能量本征值$ E_n $和调制强度$ \lambda_{\rm{D}} $的变化. 图中颜色代表γ的大小

    Fig. 6.  (a) $ \langle \text{IPR} \rangle $ (the black solid line) and $ \langle \text{NPR} \rangle $ (the red dashed line) as the functions of $ \lambda_{\rm{D}} $ for $U_{\rm{D}} = 2,\; K = 0.8, $$ L = 3194$. The grey regions denote the intermediate regimes with mobility edges; (b) fractal dimension γ of all the eigenstates as a function of energies and $ \lambda_{\rm{D}} $ for $U_{\rm{D}} = 2,$ $ K = 0.8 $, and $ L = 3194 $. Here the color code represents the values of γ

    图 7  $ U_{\rm{D}} = 2 $, $ K = 0.8 $时平均参与率$ \langle \text{NPR} \rangle $$ 1/L $变化的曲线  (a) 临界相区; (b) 局域相区

    Fig. 7.  $ \langle \text{NPR} \rangle $ as a function of $ 1/L $ for $ U_{\rm{D}} = 2 $, $K = 0.8 $: (a) critical phase region; (b) localized phase region

    图 8  (a) $\lambda_{\rm{D}}\text{-}U_{\rm{D}}$参数平面内, 以序参量η的大小为填充颜色的相图, 其中白色区域表示完全扩展或局域相, 红色区域表示具有迁移率边的临界相. 其中, 绿色方块对应$ \lambda = 3.183 $, $ \varDelta_{{\rm{c}}} = -0.2 $, $ |\alpha| = 0.604 $, $ \lambda_{{\rm{D}}}\approx2.86 $; 黑色圆圈对应的是$ \lambda = 4.113 $, $ \varDelta_{{\rm{c}}} = -0.6 $, $ |\alpha| = 0.607 $, $ \lambda_{{\rm{D}}}\approx4.86 $; 蓝色叉号对应于$ \lambda = 5.069 $, $ \varDelta_{{\rm{c}}} = -2 $, $ |\alpha| = 0.608 $, $ \lambda_{{\rm{D}}}\approx6.08 $. 这里, $ \kappa = 1.2 $, $ U = 6 $, $ U_{\rm{D}}\approx2.21 $. 相图中的(b)绿色方块, (c) 黑色圆圈和(d) 蓝色叉号对应的参数取值下, 所有态的逆参与率随本征能变化的情况. 这里, $ K = 0.8 $$ L = 1974 $

    Fig. 8.  (a) Phase diagram in the $\lambda_{\rm{D}}\text{-} U_{\rm{D}}$ plane. The color code represents the values of η, where the white regions denote the full extended or localized phase and the red region represents the critical phase. Here, the green square corresponds to $ \lambda = 3.183 $, $\varDelta_{{\rm{c}}} = -0.2$, $ |\alpha| = 0.604 $, $ \lambda_{{\rm{D}}}\approx2.86 $, the black circle corresponds to $ \lambda = 4.113 $, $ \varDelta_{{\rm{c}}} = -0.6 $, $ |\alpha| = 0.607 $, $ \lambda_{{\rm{D}}}\approx4.86 $, and the blue cross corresponds to $ \lambda = 5.069 $, $ \varDelta_{{\rm{c}}} = -2 $, $ |\alpha| = 0.608 $, $ \lambda_{{\rm{D}}}\approx6.08 $ for $ \kappa = 1.2 $, $ U = 6 $, $ U_{\rm{D}}\approx2.21 $. The IPR of different eigenstates as a function of energies for (b) the green square, (c) the black circle, and (d) the blue cross. Here, $ K = 0.8 $, and $ L = 1974 $

    图 9  (a1)$ K = 0.5 $, (b1)$ K = 1.5 $时, 腔场耦合强度U取不同的值时, 光场$ |\alpha| $随着耦合强度λ变化的图像. 这里, $ L = 1974 $, $ \beta = 610/987 $, $\varDelta_{{\rm{c}}} = -2$$ \kappa = 1.2 $. (a2)$ K = 0.5 $, (b2)$ K = 1.5 $时, $\lambda_{\rm{D}}-U_{\rm{D}}$参数平面内, 以序参量η的大小为填充颜色的相图, 其中白色区域表示完全扩展或局域相, 红色区域表示具有迁移率边的临界相. 这里, $ L = 1974 $

    Fig. 9.  The cavity field $ |\alpha| $ as a function of the pumping strength λ for different U for (a1)$ K = 0.5 $, (b1)$ K = 1.5 $. Here $ L = 1974 $, $ \beta = 610/987 $, $\varDelta_{{\rm{c}}} = -2$ and $ \kappa = 1.2 $. Phase diagram in the $\lambda_{\rm{D}}-U_{\rm{D}}$ plane with $ L = 1974 $ for (a2)$ K = 0.5 $, (b2)$ K = 1.5 $. The color code represents the values of η, where the white regions denote the full extended or localized phase and the red region represents the critical phase

  • [1]

    Marklund M, Shukla P K 2006 Rev. Mod. Phys. 78 591Google Scholar

    [2]

    Chin C, Grimm R, Julienne P, Tiesinga E 2010 Rev. Mod. Phys. 82 1225Google Scholar

    [3]

    Mekhov I B, Maschler C, Ritsch H 2007 Nat. Phys. 3 319Google Scholar

    [4]

    Ritsch H, Domokos P, Brennecke F, Esslinger T 2013 Rev. Mod. Phys. 85 553Google Scholar

    [5]

    Baumann K, Guerlin C, Brennecke F, Esslinger T 2010 Nature 464 1301Google Scholar

    [6]

    Maschler C, Mekhov I B, Ritsch H 2008 Eur. Phys. J. D 46 545Google Scholar

    [7]

    Klinder J, Keβler H, Bakhtiari M R, Thorwart M, Hemmerich A 2015 Phys. Rev. Lett. 115 230403Google Scholar

    [8]

    Landig R, Hruby L, Dogra N, Landini M, Mottl R, Donner T, Esslinger T 2016 Nature 532 476Google Scholar

    [9]

    Zheng W, Cooper N R 2018 Phys. Rev. A 97 021601(RGoogle Scholar

    [10]

    Zhou L, Pu H, Zhang K Y, Zhao X D, Zhang W P 2011 Phys. Rev. A 84 043606Google Scholar

    [11]

    Habibian H, Winter A, Paganelli S, Rieger H, Morigi G 2013 Phys. Rev. Lett. 110 075304Google Scholar

    [12]

    Anderson P W 1958 Phys. Rev. 109 1492Google Scholar

    [13]

    Aubry S, André G 1980 Ann. Isr.: Phys. Soc. 3 18

    [14]

    Roati G, D’Errico C, Fallani L, Fattori M, Fort C, Zaccanti M, Modugno G, Modugno M, Inguscio M 2008 Nature 453 895Google Scholar

    [15]

    Sarma S D, Kobayashi A, Prange R E 1986 Phys. Rev. Lett. 56 1280Google Scholar

    [16]

    Sarma S D, Song H, Xie X C 1988 Phys. Rev. Lett. 61 2144Google Scholar

    [17]

    Biddle J, Sarma S D 2010 Phys. Rev. Lett. 104 070601Google Scholar

    [18]

    Ganeshan S, Pixley J H, Sarma S D 2015 Phys. Rev. Lett. 114 146601Google Scholar

    [19]

    Lanini Y, Bromberg Y, Christodoulides D N, Silberberg Y 2010 Phys. Rev. Lett. 105 163905Google Scholar

    [20]

    Xu Z H, Xia X, Chen S 2021 arXiv: 2109.02072 v1[cond-mat.dis-nn]

    [21]

    Wang Y C, Xia X, Wang Y J, Zheng Z H, Liu X J 2021 Phys. Rev. B 103 174205Google Scholar

    [22]

    徐志浩, 皇甫宏丽, 张云波 2019 物理学报 68 087201Google Scholar

    Xu Z H, Huangfu H L, Zhang Y B 2019 Acta Phys. Sin. 68 087201Google Scholar

    [23]

    Xu Z H, Xia X, Chen S 2022 Sci. China-Phys. Mech. Astron. 65 227211Google Scholar

    [24]

    Roy S, Mishra T, Tanatar B, Basu S 2021 Phys. Rev. Lett 126 106803Google Scholar

    [25]

    Jiang X P, Qiao Y, Cao J P 2021 Chin. Phys. B 30 097202Google Scholar

    [26]

    Padhan A, Giri M K, Mondal S, Mishra T 2021 arXiv: 2109.09621 v1 [cond-mat.quant-gas]

    [27]

    Zhuang W F, Geng B, Lou H G, Guo G C, Gong M 2021 Phys. Rev. A 104 053308Google Scholar

    [28]

    Keeling J, Bhaseen M J, Simons B D 2014 Phys. Rev. Lett 112 143002Google Scholar

    [29]

    Piazza F, Strack P 2014 Phys. Rev. Lett. 112 143003Google Scholar

    [30]

    Chen Y, Yu Z H, Zhai H 2014 Phys. Rev. Lett. 112 143004Google Scholar

    [31]

    Chen Y, Zhai H, Yu Z H 2015 Phys. Rev. A 91 021602(RGoogle Scholar

    [32]

    Chen Y, Yu Z H, Zhai H 2016 Phys. Rev. A 93 041601(RGoogle Scholar

    [33]

    Cai M L, Liu Z D, Zhao W D, Wu Y K, Mei Q X, Jiang Y, He L, Zhang X, Zhou Z C, Duan L M 2021 Nat. Commun. 12 1126Google Scholar

    [34]

    Xu Z H, Chen S 2021 Phys. Rev. A 103 043325Google Scholar

  • [1] 任国梁, 申开波, 刘永佳, 刘英光. 类石墨烯氮化碳结构(C3N)热传导机理研究. 物理学报, 2023, 72(1): 013102. doi: 10.7498/aps.72.20221441
    [2] 陈奇, 戴越, 李飞燕, 张彪, 李昊辰, 谭静柔, 汪潇涵, 何广龙, 费越, 王昊, 张蜡宝, 康琳, 陈健, 吴培亨. 5—10 µm波段超导单光子探测器设计与研制. 物理学报, 2022, 71(24): 248502. doi: 10.7498/aps.71.20221594
    [3] 刘佳琳, 庞婷方, 杨孝森, 王正岭. 无序非厄米Su-Schrieffer-Heeger中的趋肤效应. 物理学报, 2022, 71(22): 227402. doi: 10.7498/aps.71.20221151
    [4] 傅聪, 叶梦浩, 赵晖, 陈宇光, 鄢永红. 共轭聚合物链中光激发过程的无序效应. 物理学报, 2021, 70(11): 117201. doi: 10.7498/aps.70.20201801
    [5] 施婷婷, 汪六九, 王璟琨, 张威. 自旋轨道耦合量子气体中的一些新进展. 物理学报, 2020, 69(1): 016701. doi: 10.7498/aps.69.20191241
    [6] 冯彦林, 樊景涛, 陈刚, 贾锁堂. 光学微腔中一维费米气的磁性关联特性. 物理学报, 2019, 68(4): 043702. doi: 10.7498/aps.68.20181954
    [7] 徐志浩, 皇甫宏丽, 张云波. 一维准周期晶格中玻色子对的迁移率边. 物理学报, 2019, 68(8): 087201. doi: 10.7498/aps.68.20182218
    [8] 刘通, 高先龙. 具有p波超流的一维非公度晶格中迁移率边研究. 物理学报, 2016, 65(11): 117101. doi: 10.7498/aps.65.117101
    [9] 王晓, 陈立潮, 刘艳红, 石云龙, 孙勇. 纵模对光子晶体中类狄拉克点传输特性的影响. 物理学报, 2015, 64(17): 174206. doi: 10.7498/aps.64.174206
    [10] 陈再高, 王建国, 王玥, 朱湘琴, 张殿辉, 乔海亮. 相对论返波管超辐射产生与辐射的数值模拟研究. 物理学报, 2014, 63(3): 038402. doi: 10.7498/aps.63.038402
    [11] 侯碧辉, 刘凤艳, 岳明, 王克军. 纳米金属镝的传导电子定域化. 物理学报, 2011, 60(1): 017201. doi: 10.7498/aps.60.017201
    [12] 丁锐, 金亚秋, 小仓久直. 二维随机介质中柱面波传播及其局域性的随机泛函分析. 物理学报, 2010, 59(6): 3674-3685. doi: 10.7498/aps.59.3674
    [13] 李晓春, 高俊丽, 刘绍娥, 周科朝, 黄伯云. 无序对二维声子晶体平板负折射成像的影响. 物理学报, 2010, 59(1): 376-380. doi: 10.7498/aps.59.376
    [14] 曹永军, 杨旭, 姜自磊. 弹性波通过一维复合材料系统的透射性质. 物理学报, 2009, 58(11): 7735-7740. doi: 10.7498/aps.58.7735
    [15] 曹永军, 杨 旭. 广义Fibonacci准周期结构声子晶体透射性质的研究. 物理学报, 2008, 57(6): 3620-3624. doi: 10.7498/aps.57.3620
    [16] 王慧琴, 刘正东, 王 冰. 二维随机介质中的能量分布和频谱特性. 物理学报, 2008, 57(9): 5550-5557. doi: 10.7498/aps.57.5550
    [17] 王慧琴, 刘正东, 王 冰. 同材质颗粒不同填充密度的随机介质中光场的空间分布. 物理学报, 2008, 57(4): 2186-2191. doi: 10.7498/aps.57.2186
    [18] 曹永军, 董纯红, 周培勤. 一维准周期结构声子晶体透射性质的研究. 物理学报, 2006, 55(12): 6470-6475. doi: 10.7498/aps.55.6470
    [19] 许兴胜, 陈弘达, 张道中. 非晶光子晶体中的光子局域化. 物理学报, 2006, 55(12): 6430-6434. doi: 10.7498/aps.55.6430
    [20] 刘小良, 徐 慧, 马松山, 宋招权. 一维无序二元固体中电子局域性质的研究. 物理学报, 2006, 55(6): 2949-2954. doi: 10.7498/aps.55.2949
计量
  • 文章访问数:  1778
  • PDF下载量:  116
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-12-05
  • 修回日期:  2022-01-18
  • 上网日期:  2022-03-04
  • 刊出日期:  2022-06-05

由超辐射引起的迁移率边和重返局域化

  • 山西大学理论物理研究所, 量子光学与光量子器件国家重点实验室, 太原 030006
  • 通信作者: 徐志浩, xuzhihao@sxu.edu.cn ; 郭利平, guolp@sxu.edu.cn
    基金项目: 国家自然科学基金(批准号: 11604188, 12147215, 11904216)、山西省基础研究计划(批准号: 20210302123442)、北京凝聚态物理国家研究中心开放课题、山西省高等学校科技创新项目(批准号: 2019L0097) 和山西省“1331工程”重点学科建设计划资助的课题

摘要: 研究了束缚在梯子形光晶格中的中性原子与高精度腔耦合的系统, 发现由超辐射引起的准周期调制可以导致迁移率边和重返局域化现象的出现. 在平均场近似下, 超辐射现象可以引起两种不同频率的准周期调制, 它们可以由腔场和泵浦场有效调节. 在本文的观测范围内, 当高频调制强度小于某一临界值时, 系统随着低频调制强度的增加经历一次局域转变. 通过数值求解分形维度、密度分布、平均参与率、平均逆参与率以及标度分析, 证明了局域转变历经的临界相区存在迁移率边; 当高频调制强度大于临界值时, 随着低频调制强度的增强, 系统依次经历完全扩展相-临界相-完全局域相-临界相-完全局域相, 这是一个典型的重返局域化现象. 最后给出了局域化相图. 该研究结果为超辐射相变和重返局域化现象的研究建立了联系, 也为重返局域化的研究搭建了新的平台.

English Abstract

    • 超冷原子气与腔量子电动力学分别是凝聚态物理与量子光学中最重要的研究领域之一, 它们有着各自的发展历史与丰富成果[1,2]. 冷原子气体与光腔的耦合系统进一步结合了这两大领域, 为模拟和研究各种新奇量子物态提供了一个重要平台[3,4], 引起了人们的广泛关注. 在单模腔中的玻色爱因斯坦凝聚体加入横向泵浦会导致体系从正常相转变为超辐射相[5]. 而在超辐射相区内, 腔内原子将受到一个腔诱导的超晶格势作用[6]. 实验上已经观测到在腔诱导的公度超晶格中存在超固体相[7,8]. 当腔诱导的晶格势与静态光晶格势是非公度的, 在平均场近似下系统形成一个有效的准周期调制, 导致金属-绝缘的转变[9-11].

      另一方面, 在这种既不是周期的, 又不是完全无序的准周期系统中, 最典型的模型就是具有准周期调制的一维Aubry-André (AA)模型, 它展现了独特的局域化性质[12,13]. 由于AA模型的自对偶特征, 具有与能量无关的局域化转变, 即在某一准周期调制强度下, 所有的单粒子态要么都是局域的, 要么都是扩展的. 2008 年, Roati等[14]在冷原子体系中利用双色光晶格的方案已经实现了准周期调制的AA 模型, 并且观测了其局域化转变. 但值得注意的是, 通过破坏自对偶特性, 一维准周期调制系统中可以出现迁移率边[15-19], 即扩展态与局域态之间的能量边界. 最近, 理论上预测了长程跃迁的准周期调制系统[20]、具有特殊形式准周期调制的系统[21]、准周期调制的玻色子对系统[22] 以及二维非厄米准晶[23]等具有简洁解析形式的迁移率边. 通常来说, 无序系统经过局域相变后, 随着无序强度的进一步增强, 所有的态始终处于局域态. 然而, 最近Roy等[24]研究了一维交叉调制的二聚体准晶系统, 发现由于二聚体和准周期调制的竞争, 导致重返局域化现象的产生, 即部分态在特殊的无序强度范围内从局域态重新变成扩展态, 随着无序强度进一步增强, 所有态又重新回到局域态的现象. 这一现象在非厄米准周期系统也被发现[25]. 重返局域化现象可以利用激发动力学手段探测[26]. 通常的冷原子实验中, 实现具有重返局域化现象的无序模型存在一定的困难.

      基于以上的研究和存在的难题, 本文讨论了在原子-光腔耦合系统中由超辐射引起的迁移率边和重返局域化, 以及实验的实现方案.

    • 图1所示, 考虑一团含有$ N_{\rm{b}} $个原子的玻色爱因斯坦凝聚体被捕陷在一个梯子形的光晶格中, 并且此梯形光晶格处于x-y平面被装载进一个高精度腔中. 每条梯子腿上的原子受到x方向上的周期势$ V_{\text{static}}=V_1\cos^2{(kx)} $限制, 其中$ V_1 $正比于光格子强度和晶格波矢k. 沿着z方向梯子腿加载波矢为$ k_{\rm{p}} $、频率为$ \omega_{{\rm{p}}} $的泵浦光. 在紧束缚近似下, 系统的哈密顿量可以写为:

      图  1  实验装置示意图. 在高精度腔中, 中性原子被x-y平面的梯子形光晶格俘获, 并沿着z方向加入驱动泵浦场. 泵浦场的频率为$\omega_{\rm{p}}$远失谐于原子能级跃迁频率$\omega_{\rm{a}}$, 但接近于腔场频率$\omega_{\rm{c}}$

      Figure 1.  A schematic diagram of experimental setup. In a high-finesses optical cavity, spinless atoms are trapped by a ladder lattice in x-y plane. The atoms are driven by a pump laser beam along z direction. The frequency $ \omega_{\rm{p}} $ of the pump laser is far detuned from the atomic transition line $ \omega_{\rm{a}} $ but close to the cavity-mode frequency $ \omega_{\rm{c}} $

      $ \begin{split} \hat{H}=\;&\hat{H}_{\text{Lad}}{-\varDelta_{{\rm{c}}}}{\hat{{a}}^{†}}{\hat{a}}+{\sum\limits_{i = 1}^{N}}{\sum\limits_{\sigma = 1,2}}[{\lambda}({\hat{a}^{†}}+{\hat{a}}){\cos(2{\pi}{\beta}i)}\\ &+U{\hat{a}^{†}}{\hat{a}}{\cos^{2}(2{\pi}{\beta}i)}] \hat{c}_{i,\sigma}^{†}\hat{c}_{i,\sigma},\\[-13pt] \end{split} $

      其中

      $ \begin{split} \hat{H}_{\text{Lad}}=\;&-J{\sum\limits_{i = 1}^{N-1}}{\sum\limits_{\sigma = 1,2}}\left(\hat{c}_{i,\sigma}^{†}\hat{c}_{i+1,\sigma}+{\rm{H.c.}}\right) \\ &-K{\sum\limits_{i = 1}^{N}}\left(\hat{c}_{i,1}^{†}\hat{c}_{i,2}+{\rm{H.c.}}\right), \end{split}$

      这里, $ \hat{a} $为腔光子的湮灭算符; $ \hat{c}_{i, \sigma} $表示在第σ条梯子腿的第i个格点上原子场湮灭算符; N为每条腿上格点的数目, 即梯子系统的总格点数目为$ L = 2 N $; 腔场频率$ \omega_{{\rm{c}}} $和泵浦光频率$ \omega_{{\rm{p}}} $的失谐量定义为$ {\varDelta_{{\rm{c}}}}={\omega_{{\rm{p}}}}-{\omega_{{\rm{c}}}} $; JK分别为沿着梯子腿和梯子横档的最近格子间的跃迁强度, 它们可以由最近邻格子间的间距所控制. 选取J作为能量单位, 即$ J = 1 $, 并且全文中固定$ K = 0.8 $进行讨论. 在绝热近似下, $ {\lambda}= {g_{{\rm{c}}}}g/{\varDelta} $为泵浦光和腔场的相干耦合强度, $ {U}={g^2_{{\rm{c}}}}/{\varDelta} $为腔和原子的耦合强度. 这里, $ g_{{\rm{c}}} $是腔模的拉比频率, g是泵浦光的拉比频率, 并且$ {\varDelta}= {\omega_{{\rm{p}}}}-{\omega_{{\rm{a}}}} $为泵浦光频率$ {\omega_{{\rm{p}}}} $和原子跃迁频率$ {\omega_{{\rm{a}}}} $的失谐量. 由于腔场具有波矢$ k_{{\rm{c}}} $, 它与晶格波矢k共同决定准周期势的调制频率$ {\beta}={k_{{\rm{c}}}/(2 k)} $. 不失一般性, 选取$ \beta=(\sqrt{5}-1)/2 $, 它可以利用斐波那契数近似地得到: $\beta=\lim\limits_{m\to \infty}(F_{m-1}/F_{m})$, 这里$ F_m $可以迭代求得$ F_{m}=F_{m-1}+F_{m-2} $, 并且$ F_0=F_1 = 1 $. 在本文的计算过程中选取系统尺寸$ L = 2 F_m $, 调制频率$ \beta\approx F_{m-1}/F_{m} $, 选取周期边界条件. 由于光腔中存在光子的耗散, 系统可以由量子主方程所描述, $ \partial_{t}\rho= - {\rm{i}}[\hat{H}, \rho]+\kappa(2\hat{a}\rho \hat{a}^{†}-\hat{a}^{†}\hat{a}\rho-\rho\hat{a}^{†}\hat{a}) $, 其中κ是腔场的耗散率. 根据参考文献[27], 在$ \omega_{\rm{c}}/N_{\rm{b}}\omega_{\rm{a}}\rightarrow 0 $的条件下, 系统可以看作平均场处理. 近年来, 许多腔场超辐射相关的理论工作都用到了平均场近似处理[28-32], 所得的结果可以很好地解释相关的实验现象[5,33]. 利用主方程, 可以得到腔场$\alpha(t)= {\langle \hat{a}(t) \rangle}$的动力学方程:

      $\begin{split} \partial_{t}{\alpha}=\;&{\rm{i}}{\varDelta_{{\rm{c}}}}{\alpha}-{\rm{i}}{\sum\limits_{i = 1}^{N}}{\sum\limits_{\sigma = 1,2}}[U{\alpha}{\cos^{2}(2{\pi}{\beta}i)} \\ &+{\lambda}{\cos(2{\pi}{\beta}i)}]\langle \hat{c}^{†}_{i,\sigma}\hat{c}_{i,\sigma}\rangle -{\kappa}{\alpha}. \end{split}$

      $ \kappa > J $时, 腔场存在稳态解. 腔场的稳态解$ (\partial_{t}\alpha = 0) $可以写为:

      $ \alpha=\frac{\lambda{\displaystyle \sum\limits_{i = 1}^{N}}{\displaystyle \sum\limits_{\sigma = 1,2}}\cos(2\pi{\beta}i)\langle \hat{c}^{†}_{i,\sigma}\hat{c}_{i,\sigma}\rangle }{\varDelta_{{\rm{c}}}+{\rm{i}}\kappa-U{\displaystyle \sum\limits_{i = 1}^{N}}{\displaystyle \sum\limits_{\sigma = 1,2}}{\cos^{2}(2\pi{\beta}i)}\langle \hat{c}^{†}_{i,\sigma}\hat{c}_{i,\sigma}\rangle }, $

      其中期望值$ \langle \hat{c}^{†}_{i, \sigma}\hat{c}_{i, \sigma}\rangle $可以通过平均场哈密顿量

      $ \hat{H}_{\text{MF}}=\hat{H}_{\text{Lad}}+\hat{V}_{\text{eff}}(\alpha) $

      自洽迭代求得. 即任选一个α的初始值, 代入平均场哈密顿量(5)式, 通过对角化哈密顿量求出系统的基态, 从而得到期望值$ \langle \hat{c}^{†}_{i, \sigma}\hat{c}_{i, \sigma}\rangle $. 将期望值$ \langle \hat{c}^{†}_{i, \sigma}\hat{c}_{i, \sigma}\rangle $代入(4)式, 求一个新的α值. 若前后两个α的差值小于所允许的误差ε, 其中$ \varepsilon = 10^{-8} $, 则认为所求的$ \langle \hat{c}^{†}_{i, \sigma}\hat{c}_{i, \sigma}\rangle $α为真实值; 若差值大于ε, 则将新得到的α值作为初始值代入系统重复前面的过程, 直到符合误差条件. 定义$\lambda_{{\rm{D}}} = 2\lambda\text{Re}({\alpha})$$ U_{{\rm{D}}}=U|{\alpha}|^{2} $, 则有效势

      $\begin{split} \hat{V}_{\text{eff}}(\alpha)=\;&{\sum\limits_{i = 1}^{N}}\sum\limits_{\sigma = 1,2}[\lambda_{\rm{D}}\cos({2\pi\beta i})\\ &+U_{\rm{D}}\cos^{2}({2\pi\beta i})]\hat{c}^{†}_{i,\sigma}\hat{c}_{i,\sigma} \end{split}$

      依赖于光场序参量α. 当处于正常相($ \alpha = 0 $)时, 系统可以由一个标准的梯子模型$ \hat{H}_{\text{Lad}} $所描述, 而当系统进入超辐射相区($ \alpha\ne 0 $), 一个由超辐射引起的调制势被引入到梯子模型. 图2(a)给出了在$ L = 1974 $, $ K = 0.8 $, $ \varDelta_{\rm{c}} = -2 $$ \kappa = 1.2 $时, 对于不同U, 光场序参量$ |\alpha| $随有效耦合强度λ的变化. 对于不同的U, 在小λ情况下, $ |\alpha| = 0 $, 而随着λ的增加到超过某个临界有效耦合强度$ \lambda_{\rm{c}} $后, 光场序参量$ |\alpha| $跳变为一个有限值, 表明此时系统进入到超辐射相区. 由图2(a)可以看出, 随着U的增加, 临界耦合强度$ \lambda_{\rm{c}} $的值逐渐增加. 图2(b)给出了$ \lambda_{\rm{c}} $随着U值的变化情况, 其中黑色实线为数值拟合结果, 即$ \lambda_{\rm{c}} = -0.1149 U^2+0.4251 U+1.9875 $. 由此可见, 通过调节腔场强度可以有效地控制超辐射转变点, 从而实现对模型中调制强度$ \lambda_{\rm{D}} $$ U_{\rm{D}} $调节的控制.

      图  2  (a) 腔场耦合强度U取不同的值时, 光场$ |\alpha| $随着耦合强度λ变化的图像; (b) 临界耦合强度$ \lambda_{\rm{c}} $随着U的变化情况. 这里, $ L = 1974 $, $ K = 0.8 $, $ \beta = 610/987 $, $ \varDelta_{{\rm{c}}} = -2 $$ \kappa = 1.2 $

      Figure 2.  (a) The cavity field $ |\alpha| $ as a function of the pumping strength λ for different U; (b) the critical pumping strength $ \lambda_{\rm{c }}$ as a function of U. Here $ L = 1974 $, $ K = 0.8 $, $ \beta = 610/987 $, $ \varDelta_{{\rm{c}}} = -2 $ and $ \kappa = 1.2 $

    • 本节首先选取参数$U_{\rm{D}} = 1,\;2$作为具体的例子, 通过改变不同的调制强度$ \lambda_{\rm{D}} $, 分别讨论存在两种不同的局域化现象, 即存在迁移率边和重返局域化. 最后给出系统在$ \lambda_{\rm{D}} $-$ U_{\rm{D}} $参数平面的局域化相图.

    • 逆参与率(IPR)和参与率(NPR)可被用来表征准周期系统中本征态的局域化特性. 对于第n个本征态, 逆参与率定义为$\text{IPR}^{(n)} = \displaystyle\sum\nolimits_{i = 1}^{N} \displaystyle\sum\nolimits_{\sigma=\{1, 2\}}\times |\psi_{i, \sigma}^{(n)}|^4$, 参与率定义为$\text{NPR}^{(n)} = \displaystyle\sum\nolimits_{i = 1}^{N}\displaystyle\sum\nolimits_{\sigma=\{1, 2\}}\times \left(L|\psi_{i, \sigma}^{(n)}|^4\right)^{-1}$. 如果本征态$ \psi^{(n)} $为局域态, 则$ \text{IPR}^{(n)}\sim L^0 $为有限值, 而$ \text{NPR}^{(n)} $正比于$ 1/L $, 在$ L\to \infty $时趋近于$ 0 $. 如果$ \psi^{(n)} $为扩展态, 则在$ L\to \infty $时, $ \text{IPR}^{(n)} \sim 1/L $, 趋近于$ 0 $. 而$ \text{NPR}^{(n)} $趋近于一个有限值. 进一步, 定义平均逆参与率$\langle \text{IPR} \rangle = 1/L \displaystyle\sum\nolimits_{n = 1}^L \text{IPR}^{(n)}$和平均参与率$\langle \text{NPR} \rangle = 1/L\displaystyle\sum\nolimits_{n = 1}^L \text{NPR}^{(n)}$, 它们携带了系统所有本征态的局域化的完全信息. 如果系统处于完全的扩展相区, 则$ \langle \text{IPR} \rangle \sim L^{-1} $$ \langle \text{NPR} \rangle $为有限值; 如果系统处于完全的局域相区, 则$ \langle \text{IPR} \rangle $为有限值且$ \langle \text{NPR} \rangle \sim L^{-1} $; 而如果系统处于临界相区, 即存在扩展和局域态的混合, 则$ \langle \text{IPR} \rangle $$ \langle \text{NPR} \rangle $将都为一个有限值[34].

      图3(a)给出了$ U_{\rm{D}} = 1 $$ L = 3194 $时, $ \langle \text{IPR} \rangle $$ \langle \text{NPR} \rangle $随着准周期调制强度$ \lambda_{\rm{D}} $的变化情况. 对比AA 模型, 此时系统的$ \langle \text{IPR} \rangle $$ \langle \text{NPR} \rangle $并不是在同一点发生转变, 而是存在了一段$ \langle \text{IPR} \rangle $$ \langle \text{NPR} \rangle $同时具有有限值的临界相区, 即图3(a)中灰色标记的区域. 这一现象表明在$ \lambda_{\rm{D}}^{({\rm{c}}1)} < \lambda_{\rm{D}} < \lambda_{\rm{D}}^{({\rm{c}}2)} $的区间, 体系能谱中应该同时存在扩展和局域态. 这里, $ \lambda_{\rm{D}}^{({\rm{c}}1)}\approx 0.85 $, $ \lambda_{\rm{D}}^{({\rm{c}}2)}\approx 3.35 $. 而当$ \lambda_{\rm{D}} > \lambda_{\rm{D}}^{({\rm{c}}2)} $, 所有本征态随着$ \lambda_{\rm{D}} $的增加都变成局域态. 图3(b)图3(d)给出了在$ U_{\rm{D}} = 1 $, 对应于不同$ \lambda_{\rm{D}} $时, $ \langle \text{NPR} \rangle $随着尺寸L的标度分析. 当$ \lambda_{\rm{D}} = 0.5 $时, 如图3(b)所示, 在$ L\to \infty $的极限下, $ \langle \text{NPR} \rangle $随着尺寸L的增大趋近于一个有限值, 约为$ 0.445 $, 表明系统中所有本征态都是扩展态. 当$ \lambda_{\rm{D}} = 5 $时, 如图3(d)所示, $ \lim\limits_{L\to\infty}\langle \text{NPR} \rangle\to 0 $, 表明此时系统中所有的本征态都为局域态. 当调制强度$ \lambda_{\rm{D}} $落在中间临界相区时, 以$ \lambda_{\rm{D }}= 2 $为例, 如图3(c)所示, $ \langle \text{NPR} \rangle $$ L\to \infty $时, 随着L呈衰减趋势, 并趋近于一个有限值, 约为0.14, 它远小于扩展相区的情况, 表明此时系统中既存在扩展态又存在局域态, 即存在迁移率边.

      图  3  (a) 当$ U_{\rm{D}} = 1, K = 0.8, L = 3194 $时, $ \langle \text{IPR} \rangle $(黑色实线) 和$ \langle \text{NPR} \rangle $(红色虚线) 随$ \lambda_{\rm{D}} $变化的曲线, 灰色区域代表着具有迁移率边的临界区域; (b)$ \lambda_{\rm{D}} = 0.5 $, (c)$ \lambda_{\rm{D}} = 2 $, (d)$ \lambda_{\rm{D}} = 5 $时, $ \langle \text{NPR} \rangle $随着$ L^{-1} $的变化, 其余参数的取值是$ U_{\rm{D}} = 1 $, $ K = 0.8 $

      Figure 3.  (a) $ \langle \text{IPR} \rangle $ (the black solid line) and $ \langle \text{NPR} \rangle $ (the red dashed line) as the functions of $ \lambda_{\rm{D}} $ for $U_{\rm{D}} = 1, K = 0.8, $$ L = 3194$. The grey region denotes the critical region with mobility edges; $ \langle \text{NPR} \rangle $ as a function of $ L^{-1} $ with $ U_{\rm{D}} = 1 $ and $ K = 0.8 $ for (b) $ \lambda_{\rm{D}} = 0.5 $, (c) $ \lambda_{\rm{D}} = 2 $, and (d) $ \lambda_{\rm{D }}= 5 $

      波函数的密度分布$\rho_i=\displaystyle\sum\nolimits_{\sigma=\{1, 2\}}|\psi_{i, \sigma}|^2$可以直观地表征系统的局域化特性, 从而进一步确定体系中迁移率边的存在. 图4给出了$ L = 1220 $$ U_{\rm{D}} = 1 $时, 对于不同$ \lambda_{\rm{D}} $的取值, 不同激发态的密度分布情况. 其中图4的左栏对应于第$ 300 $个激发态, 右栏对应于第$ 1100 $个激发态, 要求本征能量按照升序排列. 图4(a1)图4(a2) 对应于$ \lambda_{\rm{D}} = 0.5 $的情况, 可以看到此时这两个本征态在整个空间中都有分布, 展示了扩展态的特征. 图4(c1)图4(c2) 展示了$ \lambda_{\rm{D}} = 5 $时两个激发态的密度分布, 可以看到两个本征态仅在某些原胞位置有分布, 对应于局域态的特征. 而在$ \lambda_{\rm{D}}^{({\rm{c}}1)} < \lambda_{\rm{D}} < \lambda_{\rm{D}}^{({\rm{c}}2)} $区域, 当能量处于迁移率边以下时, 相应的密度分布为扩展的, 见图4(b1). 而能量处于迁移率边以上时, 相应的密度分布是局域的, 如图4(b2)所示. 此时迁移率边对应的能量指标约为$ 506 $. 由此可见, 在这一区域系统中确实存在迁移率边.

      图  4  (a1) $ \lambda_{\rm{D}} = 0.5 $, (b1) $ \lambda_{\rm{D}} = 2 $和(c1) $ \lambda_{\rm{D}} = 5 $时第300个激发态的密度分布; (a2) $ \lambda_{\rm{D }}= 0.5 $, (b2) $ \lambda_{\rm{D}} = 2 $和(c2) $ \lambda_{\rm{D}} = 5 $时第 1100 个激发态的密度分布. 这里, $ K = 0.8 $, $ U_{\rm{D}} = 1 $$ L = 1220 $

      Figure 4.  Density distributions of the 300th excited eigenstates for (a1) $ \lambda_{\rm{D}} = 0.5 $, (b1) $ \lambda_{\rm{D}} = 2 $, and (c1) $ \lambda_{\rm{D}} = 5 $; density distributions of the 1100th excited eigenstates for (a2) $ \lambda_{\rm{D}} = 0.5 $, (b2) $ \lambda_{\rm{D }}= 2 $, and (c2) $ \lambda_{\rm{D}} = 5 $. Here, $ K = 0.8 $, $ U_{\rm{D}} = 1 $, and $ L = 1220 $

      计算每个本征态的分形维度$\gamma=-\lim\limits_{L\to\infty}\times \ln{(\text{IPR})}/\ln{L}$可以更加精细地表征系统的局域化特性. 当$ \gamma\to 1 $, 对应于波函数为扩展态, 当$ \gamma \to 0 $, 对应于波函数处于局域态. 图5给出了在$ L = 3194 $$ U_{\rm{D}} = 1 $时, 不同能量本征态的分形维度γ随着无序强度$ \lambda_{\rm{D}} $的变化情况. 可以看出, 系统中存在两个局域化转变点$ \lambda_{\rm{D}}^{({\rm{c}}1)} $$ \lambda_{\rm{D}}^{({\rm{c}}2)} $分别由黑色实线和黑色虚线标出. 当$ \lambda_{\rm{D}} < \lambda_{\rm{D}}^{({\rm{c}}1)} $时, 所有的本征态的分形维度γ都趋近于$ 1 $, 表明此时系统处于完全的扩展相. 当$ \lambda_{\rm{D}} > \lambda_{\rm{D}}^{({\rm{c}}2)} $时, 所有的本征态的分形维度γ都趋近于$ 0 $, 对应于系统处于完全局域相. 而在中间区域, 存在扩展与局域态的混合, 即存在迁移率边.

      图  5  在$ U_{\rm{D}} = 1 $, $ K = 0.8 $$ L = 3194 $时, 分形维度γ随着调制强度$ \lambda_{\rm{D }}$和能量本征值$ E_n $的变化. 其中, 黑色实线和黑色虚线分别对应了两个局域化转变点$ \lambda_{\rm{D}}^{({\rm{c}}1)} $$ \lambda_{\rm{D}}^{({\rm{c}}2)} $. 图中颜色代表着γ的大小

      Figure 5.  γ of the eigenstates as a function of the energy spectrum and $ \lambda_{\rm{D}} $ for $ U_{\rm{D}} = 1 $, $ K = 0.8 $, and $ L = 3194 $. Here the black solid line and black dashed line denote two transport points, $ \lambda_{\rm{D}}^{({\rm{c}}1)} $ and $ \lambda_{\rm{D}}^{({\rm{c}}2)} $. The color code represents the values of γ

    • $ U_{\rm{D}} = 2 $时, 图6(a)给出了$ \langle \text{IPR} \rangle $$ \langle \text{NPR} \rangle $随着准周期调制强度$ \lambda_{\rm{D}} $的变化情况. 在图6(a)中, 随着$ \lambda_{\rm{D}} $的增加, 存在两个$ \langle \text{IPR} \rangle $$ \langle \text{NPR} \rangle $的共存区间. 当$ \lambda_{\rm{D}} < \lambda_{\rm{D}}^{({\rm{c}}1)}\approx 0.3 $时, $ \langle \text{IPR} \rangle $趋近于$ 0 $, 而$ \langle \text{NPR} \rangle $为有限值, 表明此时系统处于完全扩展相. 当$ \lambda_{\rm{D}}^{({\rm{c}}1)} < \lambda_{\rm{D}} < \lambda_{\rm{D}}^{({\rm{c}}2)}\approx 2.14 $时, $ \langle \text{IPR} \rangle $$ \langle \text{NPR} \rangle $都保持有限值, 表明系统中存在扩展与局域态的混合, 即存在迁移率边. 当$ \lambda_{\rm{D}} > \lambda_{\rm{D}}^{({\rm{c}}2)} $时, 平均逆参与率为有限值而平均参与率趋近于零, 表明系统进入到完全局域相. 通常来说, 进一步增加调制强度$ \lambda_{\rm{D}} $, 系统将一直处于局域区, 如$ U_{\rm{D}} = 1 $时的情况. 然而通过计算发现, 在$ \lambda_{\rm{D}}\in (\lambda_{\rm{D}}^{({\rm{c}}3)}, \lambda_{\rm{D}}^{({\rm{c}}4)}) $区间, $ \langle \text{IPR} \rangle $$ \langle \text{NPR} \rangle $又重新恢复到都保持有限值的情况, 表明此时系统又重新进入到具有迁移率边的相区. 这里, $ \lambda_{\rm{D}}^{({\rm{c}}3)} \approx 2.85 $$ \lambda_{\rm{D}}^{({\rm{c}}4)}\approx 4.88 $. 而当$\lambda_{\rm{D}} > \lambda_{\rm{D}}^{({\rm{c}}4)}$, 系统又进入完全局域相. 这是典型的重返局域化的特性.

      图  6  当$ U_{\rm{D}} = 2, K = 0.8, L = 3194 $时 (a)$ \langle \text{IPR} \rangle $(黑色实线)和$ \langle \text{NPR} \rangle $ (红色虚线)随着$ \lambda_{\rm{D}} $变化的曲线, 灰色区域表示具有迁移率边的临界相区; (b)分形维度γ随着能量本征值$ E_n $和调制强度$ \lambda_{\rm{D}} $的变化. 图中颜色代表γ的大小

      Figure 6.  (a) $ \langle \text{IPR} \rangle $ (the black solid line) and $ \langle \text{NPR} \rangle $ (the red dashed line) as the functions of $ \lambda_{\rm{D}} $ for $U_{\rm{D}} = 2,\; K = 0.8, $$ L = 3194$. The grey regions denote the intermediate regimes with mobility edges; (b) fractal dimension γ of all the eigenstates as a function of energies and $ \lambda_{\rm{D}} $ for $U_{\rm{D}} = 2,$ $ K = 0.8 $, and $ L = 3194 $. Here the color code represents the values of γ

      重返局域化的特征可以由每个能量对应的本征态的分形维度γ随着$ \lambda_{\rm{D}} $的变化来精确地反映, 如图6(b)所示. 其中, 4条灰色实线从左到右分别对应于局域化转变点$ \lambda_{\rm{D}}^{({\rm{c}}1)} $, $ \lambda_{\rm{D}}^{({\rm{c}}2)} $, $ \lambda_{\rm{D}}^{({\rm{c}}3)} $$ \lambda_{\rm{D}}^{({\rm{c}}4)} $. 重返局域化现象指的是系统经过第一次局域转变, 已经处于完全局域相后, 随着调制强度的进一步增大, 部分态再次回到扩展态; 最终在更大的无序强度下, 系统重新返回到完全局域相. 从图6(b)可以直观地看出, $ \lambda_{\rm{D}} < \lambda_{\rm{D}}^{({\rm{c}}1)} $时, 所有态的γ值都趋于$ 1 $, 系统处于完全扩展相. $ \lambda_{\rm{D}}^{({\rm{c}}1)} < \lambda_{\rm{D}} < \lambda_{\rm{D}}^{({\rm{c}}2)} $时, 部分态的γ值趋于$ 0 $, 部分态的γ值趋于$ 1 $, 系统处于两种态并存的临界相区. 在$ \lambda_{\rm{D}}^{({\rm{c}}2)} < \lambda_{\rm{D}} < \lambda_{\rm{D}}^{({\rm{c}}3)} $的区间内, 所有态的γ都趋于$ 0 $, 系统经历局域转变, 已经处于完全局域相. 而在$ \lambda_{\rm{D}}^{({\rm{c}}3)} < \lambda_{\rm{D}} < \lambda_{\rm{D}}^{({\rm{c}}4)} $区间内, 部分本征态对应的γ的值趋于$ 1 $, 代表着这些态重新回到扩展态, 最终在$ \lambda_{\rm{D}} > \lambda_{\rm{D}}^{({\rm{c}}4)} $后, 所有态对应的γ都趋近于$ 0 $, 系统又一次返回到完全局域相, 该过程更加清晰地描述了重返局域化现象.

      为了进一步验证重返局域化现象是否存在, 在临界相区和局域相区, 分别研究了$ \langle \text{NPR} \rangle $的标度变化. 在图7(a)中, 展示了$ \lambda_{\rm{D}} = 1,\;4 $时, $ \langle \text{NPR} \rangle $随尺寸L变化的情况, 可以看出当$ L\to\infty $时, $ \langle \text{NPR} \rangle $均趋近于有限值, 表明此时系统处于扩展和局域混合的临界相区. 并且在热力学极限情况下, $ \lambda_{\rm{D}} = 1 $时的$ \langle \text{NPR} \rangle $要比$ \lambda_{\rm{D}} = 4 $时大很多, 这表明系统处于$ \lambda_{\rm{D}} = 1 $时的扩展态所占比例要远大于$ \lambda_{\rm{D}} = 4 $时的比例. 而图7(b)展示了$ \lambda_{\rm{D}} = 2.5,\;5.5 $$ \langle \text{NPR} \rangle $随尺寸L的变化情况. 可以看出, 当$ L\to\infty $时, $ \langle \text{NPR} \rangle $均趋于零, 展示了完全局域的特性. 由此可见, 在$ U_{\rm{D}} = 2 $时, 随着调制强度$ \lambda_{\rm{D}} $的增加, 系统出现了重返局域化现象.

      图  7  当$ U_{\rm{D}} = 2 $, $ K = 0.8 $时平均参与率$ \langle \text{NPR} \rangle $$ 1/L $变化的曲线  (a) 临界相区; (b) 局域相区

      Figure 7.  $ \langle \text{NPR} \rangle $ as a function of $ 1/L $ for $ U_{\rm{D}} = 2 $, $K = 0.8 $: (a) critical phase region; (b) localized phase region

    • 利用$ \langle \text{IPR} \rangle $$ \langle \text{NPR} \rangle $定义$ \eta = \ln{\left[ \langle \text{IPR} \rangle \langle \text{NPR} \rangle \right]} $, 可以把临界相从完全的扩展和局域相中清晰地区分出来. 这是因为在临界相时η的值大于$ -\ln{L} $, 而处于完全扩展和局域的相时, η的值小于$ -\ln{L} $. 图8(a)给出了$ \lambda_{\rm{D}}\text{-} U_{\rm{D}} $的平面内, 以η的大小为填充颜色的相图. 其中, 白色区域代表了完全局域相或完全扩展相, 红色区域代表了具有迁移率边的临界相, 黑色实线对应于临界值$ U_{\rm{D}}^{({\rm{c}})}\approx 1.97 $. 在本文所研究的调制强度$ U_{\rm{D}} $的取值范围内, 可以看到当$ U_{\rm{D}} < U_{\rm{D}}^{({\rm{c}})} $时, 系统随着$ \lambda_{\rm{D}} $的增大, 从完全扩展相经历一个扩展态和局域态共存的临界相, 最终转变到完全局域相. 而当$ U_{\rm{D}} > U_{\rm{D}}^{({\rm{c}})} $时, 随着$ \lambda_{\rm{D}} $的增大, 系统经历第一次局域转变到达完全局域相后, 又重新回到临界相区, 最终经历第二次局域转变到达完全局域相. 因此, 当$ U_{\rm{D}} > U_{\rm{D}}^{({\rm{c}})} $时, 体系中存在重返局域化现象.

      图  8  (a) $\lambda_{\rm{D}}\text{-}U_{\rm{D}}$参数平面内, 以序参量η的大小为填充颜色的相图, 其中白色区域表示完全扩展或局域相, 红色区域表示具有迁移率边的临界相. 其中, 绿色方块对应$ \lambda = 3.183 $, $ \varDelta_{{\rm{c}}} = -0.2 $, $ |\alpha| = 0.604 $, $ \lambda_{{\rm{D}}}\approx2.86 $; 黑色圆圈对应的是$ \lambda = 4.113 $, $ \varDelta_{{\rm{c}}} = -0.6 $, $ |\alpha| = 0.607 $, $ \lambda_{{\rm{D}}}\approx4.86 $; 蓝色叉号对应于$ \lambda = 5.069 $, $ \varDelta_{{\rm{c}}} = -2 $, $ |\alpha| = 0.608 $, $ \lambda_{{\rm{D}}}\approx6.08 $. 这里, $ \kappa = 1.2 $, $ U = 6 $, $ U_{\rm{D}}\approx2.21 $. 相图中的(b)绿色方块, (c) 黑色圆圈和(d) 蓝色叉号对应的参数取值下, 所有态的逆参与率随本征能变化的情况. 这里, $ K = 0.8 $$ L = 1974 $

      Figure 8.  (a) Phase diagram in the $\lambda_{\rm{D}}\text{-} U_{\rm{D}}$ plane. The color code represents the values of η, where the white regions denote the full extended or localized phase and the red region represents the critical phase. Here, the green square corresponds to $ \lambda = 3.183 $, $\varDelta_{{\rm{c}}} = -0.2$, $ |\alpha| = 0.604 $, $ \lambda_{{\rm{D}}}\approx2.86 $, the black circle corresponds to $ \lambda = 4.113 $, $ \varDelta_{{\rm{c}}} = -0.6 $, $ |\alpha| = 0.607 $, $ \lambda_{{\rm{D}}}\approx4.86 $, and the blue cross corresponds to $ \lambda = 5.069 $, $ \varDelta_{{\rm{c}}} = -2 $, $ |\alpha| = 0.608 $, $ \lambda_{{\rm{D}}}\approx6.08 $ for $ \kappa = 1.2 $, $ U = 6 $, $ U_{\rm{D}}\approx2.21 $. The IPR of different eigenstates as a function of energies for (b) the green square, (c) the black circle, and (d) the blue cross. Here, $ K = 0.8 $, and $ L = 1974 $

      为了考察实际的腔场参数与重返局域化现象的关系, 在图8(a)中选取3种不同的点去计算它们的IPR, 通过数值计算可以得到此时图8(a)中绿色方块、黑色圆圈、蓝色叉号处所对应的$ |\alpha| $分别为$ 0.604 $, $ 0.607 $$ 0.608 $, $\lambda_{{\rm{D}}}$分别为$ 2.86 $, $ 4.86 $$ 6.08 $, 且具有相同的$U_{{\rm{D}}}\approx2.21$. 为了进一步验证此时这3个点所对应的局域化性质, 图8(b)(d)分别给出了它们的IPR随着全部能量的变化情况. 可以看出$ |\alpha| = 0.604 $, $ \lambda_{{\rm{D}}}\approx2.86 $时, 系统已经处于完全局域相. 而在$ |\alpha| = 0.607 $, $ \lambda_{{\rm{D}}}\approx 4.86 $时, 部分态处于局域态, 部分态处于扩展态, 系统处于临界相. 最后, 在$|\alpha| = 0.608$, $ \lambda_{{\rm{D}}}\approx6.08 $时, 所有态又重新回到局域态. 由此可见, 通过合理地调节腔场的参数, 在超辐射相区可以实现重返局域化现象.

      之前讨论了$ K = 0.8 $时的情况, 图9分别给出了$ K = 0.5 $$ K = 1.5 $时, 超辐射相变与局域化相图的情况. 可以看出, 当K的取值在有限的变化范围内(以$ K = 0.5 $(图9(a1))和$ K = 1.5 $(图9(b1))为例), 对于不同的腔场耦合强度U, 光场$ |\alpha| $随泵浦和腔的耦合强度λ的变化具有类似的转变行为. 对于不同的K值, 如图9(a2)图9(b2)所示的$ K = 0.5 $$ K = 1.5 $情况, 展示了与前文中$ K = 0.8 $类似的局域化相图. 由此可见, 当K取值与J处于相同量级时, 系统展示了类似的超辐射和局域化相变.

      图  9  (a1)$ K = 0.5 $, (b1)$ K = 1.5 $时, 腔场耦合强度U取不同的值时, 光场$ |\alpha| $随着耦合强度λ变化的图像. 这里, $ L = 1974 $, $ \beta = 610/987 $, $\varDelta_{{\rm{c}}} = -2$$ \kappa = 1.2 $. (a2)$ K = 0.5 $, (b2)$ K = 1.5 $时, $\lambda_{\rm{D}}-U_{\rm{D}}$参数平面内, 以序参量η的大小为填充颜色的相图, 其中白色区域表示完全扩展或局域相, 红色区域表示具有迁移率边的临界相. 这里, $ L = 1974 $

      Figure 9.  The cavity field $ |\alpha| $ as a function of the pumping strength λ for different U for (a1)$ K = 0.5 $, (b1)$ K = 1.5 $. Here $ L = 1974 $, $ \beta = 610/987 $, $\varDelta_{{\rm{c}}} = -2$ and $ \kappa = 1.2 $. Phase diagram in the $\lambda_{\rm{D}}-U_{\rm{D}}$ plane with $ L = 1974 $ for (a2)$ K = 0.5 $, (b2)$ K = 1.5 $. The color code represents the values of η, where the white regions denote the full extended or localized phase and the red region represents the critical phase

    • 本文研究了在高精度光腔中由超辐射引起的准周期调制的梯子模型. 通过数值计算发现系统中存在迁移率边以及重返局域化现象, 并且得到了全参数空间的相图. 在本文研究的调制强度$U_{\rm{D}}$范围内, 系统的相图被临界值$U_{\rm{D}}^{({\rm{c}})}$分为两部分, 在小于临界值$U_{\rm{D}}^{({\rm{c}})}$的区间, 系统随着$\lambda_{\rm{D}}$的增大, 仅仅出现一次具有迁移率边的临界相; 而在大于临界值$ U_{\rm{D}}^{({\rm{c}})} $的部分, 系统随着$ \lambda_{\rm{D}} $的增大, 会出现重返局域化现象.

      重返局域化现象最初是在具有奇偶交错的准周期调制的Su-Schrieffer-Hegger (SSH) 模型中发现的[24], 并且重返局域化现象出现的区间非常狭小, 这给实验探测带来了极大的挑战. 而对于最初的理论方案中如果选取一个无奇偶交替的准周期调制时, 此一维晶格并不会出现重返局域化现象, 仅仅能观测到迁移率边的存在. 本文提出了一种梯子形冷原子体系与光腔耦合的方案, 此方案并不要求奇偶交替的调制并且具有有限的重返局域化相区. 以碱金属原子23Na为例实现该理论模型, 其中钠原子被捕陷在一个晶格长度约为$ 266\;{\rm{ nm}} $的梯子形光晶格中, 其晶格隧穿强度约为$ 2\pi\times200 $ kHz, 并且选用腔长约为$ 300\;\mu {\rm{m}}—1\; {\rm{mm}} $, 耗散率κ约为$ 2\pi\times240 $ kHz的高精度腔, 泵浦光的频率与腔场频率相差$\varDelta_{{\rm{c}}}\approx 2\pi\times140—2\pi\times600$ kHz. 通过调节腔场和泵浦场各个参数, 实现耦合强度λ和腔及原子的耦合强度U达到$ 2\pi\times800 $ kHz, 这些参数是目前实验上可实现的. 值得注意的是, 在平均场近似下, 当超辐射相变发生时, 梯子模型将出现两种不同频率的调制, 而这两种不同频率的准周期调制的共同作用是出现重返局域化的主要原因. 本文所研究的重返局域化现象的出现并不依赖于所填充的原子种类. 但是由于费米子体系存在泡利不相容原理, 造成发生超辐射相变的难易程度与原子填充数有关[28-32], 使得计算过程更为复杂. 本文为了理论计算的简便仅关注了玻色原子的情况. 该工作为原子与光腔耦合系统和重返局域化现象的研究建立了联系, 也为重返局域化的研究提供了新的思路.

参考文献 (34)

目录

    /

    返回文章
    返回