-
研究了一维复相互作用调制的非厄米玻色子模型. 通过数值计算能谱的实-复转变、Shannon熵、标准参与比率与拓扑缠绕数发现, 当相互作用强度低于临界相互作用强度时, 系统的能谱全为实数, 处于扩展相, 且系统是拓扑平庸的; 而当相互作用强度超过临界相互作用强度时, 系统开始出现复能谱, 处于扩展态与局域态混合相, 且此时系统是拓扑非平庸的. 计算结果表明, 能谱的实-复转变点、扩展-局域的转变点与拓扑转变点相一致. 动力学演化结果可以验证系统的实-复转变与局域化转变. 最后, 提出利用二维光子波导阵列可以模拟这一复相互作用调制的一维玻色子模型. 此项工作将为非厄米两体系统的局域性质提供很好的参考.
In this work, we investigate a one-dimensional two-boson system with complex interaction modulation, described by the Hamiltonian: $\hat{H}=-J\displaystyle\sum\nolimits_{j}\left(\hat{c}_j^\dagger\hat{c}_{j+1}+{\rm h.c}\right)+\sum\nolimits_{j}\frac{U}{2}{\rm e}^{2{\rm i}\pi\alpha j}\hat{n}_j\left(\hat{n}_j-1\right), $ where U is the interaction amplitude, and the modulation frequency $\alpha=(\sqrt{5}-1)$ is an irrational number. The interaction satisfies $U_{-j}=U^*_j$, which ensures that the system possesses party-time (PT) reversal symmetry. Using the exact diagonalization method, we numerically calculate the real-to-complex transition of the energy spectrum, Shannon entropy, the normalized participation ration, and the topological winding number. For small U, all eigenvalues are real. However, as U increases, eigenvalues corresponding to two particles occupying the same site become complex, marking a PT symmetry-breaking transition at $U=2$. This point signifies a real-to-complex transition in the spectrum. To characterize the localization properties of the system, we employ the Shannon entropy and the normalized participation ration (NPR). When $U<2$, all the eigenstates are extended, exhibiting high Shannon entropy and NPR values. Conversely, for $U>2$, states with complex eigenvalues show low Shannon entropy and significantly reduced NPR, indicating localization. Meanwhile, states with real eigenvalues remain extended in this regime. We further analyze the topological aspects of the system by using the winding number. A topological phase transition occurs at $U=2$, where the winding number changes from 0 to 1. This transition coincides with the onset of PT symmetry breaking and the localization transition. The dynamical evolution can be used to detect the localization properties and the real-to-complex transition, with the initial state being two bosons occupying the center site of the chain simultaneously. Finally, we propose an experimental realization by using a two-dimensional linear photonic waveguide array. The modulated interaction can be controlled by adjusting the real part and imaginary part of the refractive index of diagonal waveguide. To simulate this non-Hermitian two-body problem, we numerically calculate the density distribution of the wave packet in a two-dimensional plane, which indirectly reflects the propagation of light in a two-dimensional waveguide array. We hope that our work can deepen the understanding of the relation between interaction and disorder while arousing further interest in two-body systems and non-Hermitian localization. -
Keywords:
- non-Hermitian /
- localization /
- topology /
- party-time symmetry
[1] Anderson P W 1958 Phys. Rev. 109 1492
Google Scholar
[2] Billy J, Josse V, Zuo Z, Bernard A, Hambrecht B, Lugan P, Clément D, Sanchez-Palencia L, Bouyer P, Aspect A 2008 Nature 453 891
Google Scholar
[3] Roati G, D'Errico C, Fallani L, Fattori M, Fort C, Zaccanti M, Modugno G, Modugno M, Inguscio M 2008 Nature 453 895
Google Scholar
[4] Chabanov A A, Stoytchev M, Genack A Z 2000 Nature 404 850
Google Scholar
[5] Pradhan P, Sridhar S 2000 Phys. Rev. Lett. 85 2360
Google Scholar
[6] Lahini Y, Pugatch R, Pozzi F, Sorel M, Morandotti R, Davidson N, Silberberg Y 2009 Phys. Rev. Lett. 103 013901
Google Scholar
[7] Liu J H, Xu Z H 2023 Phys. Rev. B 108 184205
Google Scholar
[8] 徐志浩, 皇甫宏丽, 张云波 2019 68 087201
Google Scholar
Xu Z H, Huang F H L, Zhang Y B 2019 Acta Phys. Sin. 68 087201
Google Scholar
[9] Bender Carl M, Stefan Boettcher 1998 Phys. Rev. Lett. 80 5243
Google Scholar
[10] Heiss W D 2012 Phys. Rev. A 45 444016
Google Scholar
[11] Miri M A, Alù A 2019 Science 363 086803
Google Scholar
[12] Yao S, Song F, Wang Z 2018 Phys. Rev. Lett. 121 136802
Google Scholar
[13] Ou Z, Wang Y, Li L 2023 Phys. Rev. B 107 L161404
Google Scholar
[14] Yao S, Wang Z 2018 Phys. Rev. Lett. 121 086803
Google Scholar
[15] Borgnia D S, Kruchkov A J, Slager R J 2020 Phys. Rev. Lett. 124 056802
Google Scholar
[16] Feinberg J, Zee A 1999 Phys. Rev. E 59 6433
Google Scholar
[17] Kawabata K, Shiozaki K, Ueda M, Sato M 2019 Phys. Rev. X 9 041015
Google Scholar
[18] Hatano N, Nelson D R 1996 Phys. Rev. Lett. 77 570
Google Scholar
[19] Hatano N, Nelson D R 1997 Phys. Rev. B 56 8651
Google Scholar
[20] Hatano N, Nelson D R 1998 Phys. Rev. B 58 8384
Google Scholar
[21] Longhi, Stefano 2019 Phys. Rev. B 100 125157
Google Scholar
[22] Hamazaki R, Kawabata K, Ueda M 2019 Phys. Rev. Lett. 123 090603
Google Scholar
[23] Zhai L J, Yin S, Huang G Y 2020 Phys. Rev. B 102 064206
Google Scholar
[24] Tomita T, Nakajima S, Danshita I, Takasu Y, Takahashi Y 2017 Sci. Adv. 3 e1701513
Google Scholar
[25] Sponselee K, Freystatzky L, Abeln B, et al. 2018 Quantum Sci. Technol 4 0140027
Google Scholar
[26] Lee C H 2021 Phys. Rev. B 104 195102
Google Scholar
[27] Shen R, Lee C H 2022 Commun. Phys. 5 238
Google Scholar
[28] Gao M, Sheng C, Zhao Y, et al. 2024 Phys. Rev. B 110 094308
Google Scholar
[29] Corrielli G, Crespi A, Della Valle G, Longhi S, Osellame R 2013 Nat. Commun. 4 1555
Google Scholar
[30] XingY, Zhao X D, Lü Z, et al. 2021 Optics Express 29 40428
Google Scholar
[31] Wang Y C, Hu H P, Chen S 2016 Eur. Phys. J B 89 1
Google Scholar
[32] Zhou L W, Han W Q 2021 Chin. Phys. B 30 100308
Google Scholar
[33] Tian Q, Gu Y J, Zhou L W 2024 Phys. Rev. B 109 054204
Google Scholar
-
图 1 (a) $ U=1 $和(b) $ U=3 $时系统的能谱; (c)不同尺寸下能量虚部的绝对值的最大值随U的变化; (d) $ U=1 $ 时, 能量虚部最大值所对应本征态的两玻色子的坐标关联; (e) $ U=3 $时, 复能量环上两玻色子的坐标关联
Fig. 1. Complex energy spectrum for (a) $ U=1 $ and (b) $ U=3 $, respectively; (c) the maximum value of the absolute value of the imaginary part of energy as a function of U for different L; (d) the coordinate correlation of two bosons corresponding to the eigenstate with the maximum value of the imaginary part of energy for $ U=1 $; (e) the coordinate correlation of two bosons on the complex energy ring for $ U=3 $.
图 2 (a) $ U=1 $和(b) $ U=3 $时系统的Shannon熵随能量本征值的变化; (c) $ \eta(E) $在不同的相互作用U下随能量本征值$ E_n $的实部的变化; (d) $ \eta(E_n) $在不同的相互作用U下随能量本征值$ E_n $的虚部的变化; 这里, $ L=89 $
Fig. 2. The Shannon entropy as the function of $ \mathrm{Re}(E_n) $ for (a) $ U=1 $ and (b) $ U=3 $, respectively; (c) $ \eta(E_n) $ as the function of $ \mathrm{Re}(E_n) $; (d) $ \eta(E_n) $ as the function of $ \mathrm{Im}(E_n) $. Here, $ L=89 $.
图 5 (a) 二维光子波导阵列的示意图(每个圆圈代表一个波导. 橘黄色空圆圈和绿色圆圈分别标示具有不同折射率的波导, 从而有效地实现复调制的在位相互作用. 这个二维几何结构关于对角线对称); (b) $ U=1 $和(c) $ U=3 $时$ t=200 $时刻粒子在二维平面内的密度扩散
Fig. 5. (a) Schematic diagram of the two-dimensional linear photonic waveguide array. Each circle represents a waveguide. Orange-yellow hollow and green-colored circles label waveguides with different refractive indices, which can effectively control the complex modulated on-site interactions. This two-dimensional geometric structure is symmetric with respect to the diagonal. The density distributions of the wave packets in the 2D plane at time $ t = 200 $ for (b) $ U=1 $ and (c) $ U=3 $, respectively.
-
[1] Anderson P W 1958 Phys. Rev. 109 1492
Google Scholar
[2] Billy J, Josse V, Zuo Z, Bernard A, Hambrecht B, Lugan P, Clément D, Sanchez-Palencia L, Bouyer P, Aspect A 2008 Nature 453 891
Google Scholar
[3] Roati G, D'Errico C, Fallani L, Fattori M, Fort C, Zaccanti M, Modugno G, Modugno M, Inguscio M 2008 Nature 453 895
Google Scholar
[4] Chabanov A A, Stoytchev M, Genack A Z 2000 Nature 404 850
Google Scholar
[5] Pradhan P, Sridhar S 2000 Phys. Rev. Lett. 85 2360
Google Scholar
[6] Lahini Y, Pugatch R, Pozzi F, Sorel M, Morandotti R, Davidson N, Silberberg Y 2009 Phys. Rev. Lett. 103 013901
Google Scholar
[7] Liu J H, Xu Z H 2023 Phys. Rev. B 108 184205
Google Scholar
[8] 徐志浩, 皇甫宏丽, 张云波 2019 68 087201
Google Scholar
Xu Z H, Huang F H L, Zhang Y B 2019 Acta Phys. Sin. 68 087201
Google Scholar
[9] Bender Carl M, Stefan Boettcher 1998 Phys. Rev. Lett. 80 5243
Google Scholar
[10] Heiss W D 2012 Phys. Rev. A 45 444016
Google Scholar
[11] Miri M A, Alù A 2019 Science 363 086803
Google Scholar
[12] Yao S, Song F, Wang Z 2018 Phys. Rev. Lett. 121 136802
Google Scholar
[13] Ou Z, Wang Y, Li L 2023 Phys. Rev. B 107 L161404
Google Scholar
[14] Yao S, Wang Z 2018 Phys. Rev. Lett. 121 086803
Google Scholar
[15] Borgnia D S, Kruchkov A J, Slager R J 2020 Phys. Rev. Lett. 124 056802
Google Scholar
[16] Feinberg J, Zee A 1999 Phys. Rev. E 59 6433
Google Scholar
[17] Kawabata K, Shiozaki K, Ueda M, Sato M 2019 Phys. Rev. X 9 041015
Google Scholar
[18] Hatano N, Nelson D R 1996 Phys. Rev. Lett. 77 570
Google Scholar
[19] Hatano N, Nelson D R 1997 Phys. Rev. B 56 8651
Google Scholar
[20] Hatano N, Nelson D R 1998 Phys. Rev. B 58 8384
Google Scholar
[21] Longhi, Stefano 2019 Phys. Rev. B 100 125157
Google Scholar
[22] Hamazaki R, Kawabata K, Ueda M 2019 Phys. Rev. Lett. 123 090603
Google Scholar
[23] Zhai L J, Yin S, Huang G Y 2020 Phys. Rev. B 102 064206
Google Scholar
[24] Tomita T, Nakajima S, Danshita I, Takasu Y, Takahashi Y 2017 Sci. Adv. 3 e1701513
Google Scholar
[25] Sponselee K, Freystatzky L, Abeln B, et al. 2018 Quantum Sci. Technol 4 0140027
Google Scholar
[26] Lee C H 2021 Phys. Rev. B 104 195102
Google Scholar
[27] Shen R, Lee C H 2022 Commun. Phys. 5 238
Google Scholar
[28] Gao M, Sheng C, Zhao Y, et al. 2024 Phys. Rev. B 110 094308
Google Scholar
[29] Corrielli G, Crespi A, Della Valle G, Longhi S, Osellame R 2013 Nat. Commun. 4 1555
Google Scholar
[30] XingY, Zhao X D, Lü Z, et al. 2021 Optics Express 29 40428
Google Scholar
[31] Wang Y C, Hu H P, Chen S 2016 Eur. Phys. J B 89 1
Google Scholar
[32] Zhou L W, Han W Q 2021 Chin. Phys. B 30 100308
Google Scholar
[33] Tian Q, Gu Y J, Zhou L W 2024 Phys. Rev. B 109 054204
Google Scholar
计量
- 文章访问数: 450
- PDF下载量: 15
- 被引次数: 0