搜索

x

留言板

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

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

基于U(1)对称的无限矩阵乘积态张量网络算法提取Luttinger液体参数K

王秀娟 李生好

引用本文:
Citation:

基于U(1)对称的无限矩阵乘积态张量网络算法提取Luttinger液体参数K

王秀娟, 李生好

Extracting Luttinger liquid parameter K based on U(1) symmetric infinite matrix product states

Wang Xiu-Juan, Li Sheng-Hao
PDF
HTML
导出引用
  • 本文数值研究了自旋$ S=1/2,1,2 $的各向异性量子XXZD模型的Luttinger液体参数K. 首先, 利用$ U(1) $对称的无限矩阵乘积态算法(iMPS)得到在Luttinger液体相中的基态波函数. 通过二分量子涨落F和有限纠缠标度指数$ \kappa $的关系可以提取出Luttinger液体参数K. 对于自旋$ S=1/2, D=0 $的量子XXZD模型, 本文利用$ U(1) $对称的iMPS的算法得到的数值结果与精确解符合得很好. 在参数$ D\leqslant -2 $的区域, 自旋$ S=1 $的XXZD模型的哈密顿量可以被映射到一个自旋$ S=1/2 $的有效XXZ模型, 本文计算了在这个区域内的Luttinger液体参数K与精确解基本是一致的, 相对误差小于$ 1\% $. 此外, 在参数$ \varDelta=-0.5, D=0 $处, 本文数值计算的Luttinger液体参数与密度矩阵重整化群(DMRG)的结果也是一致的. 这些研究结果表明: 当系统具有$ U(1) $对称性时, 利用$ U(1) $对称的iMPS的方法可以提取无能隙相中的Luttinger液体参数. 本文利用此方法还研究了自旋$ S=1 $的XXZD模型在其他参数下的Luttinger液体参数, 以及自旋$ S=2 $的XXZD模型的Luttinger液体参数.
    We numerically calculate Luttinger liquid parameter K in the anisotropic spin XXZD models with spin $s = 1/2$, 1, and 2. In order to obtain groundstate wavefunctions in Luttinger liquid phases, we employ the $U(1)$ symmetric infinite matrix product states algorithm (iMPS). By using relation between the bipartite quantum fluctuations F and the so-called finite-entanglement scaling exponents $\kappa$, the Luttinger liquid parameter K can be extracted. For $s = 1/2$ and $D=0$, the numerically extracted Luttinger liquid parameter K is shown to be good agreement with the exact value. On using the fact that the spin-1 XXZD Hamiltonian with $ D \leqslant - 2$ can be mapped to an effective spin-1/2 XXZ model, we calculate the Luttinger liquid parameter for the region of $ D \leqslant - 2$. It is shown that our numerical value of the Luttinger liquid parameter agree well with the exact values, here, the relative error less than $1\%$. Also, our Luttinger liquid parameter at $\Delta = - 0.5$ and $ D = 0$ is shown to be consistent with the result form the density matrix renormalization group (DMRG) method. These results suggest that the $U(1)$ symmetric iMPS method can be applicable to calculate Luttinger liquid parameters if any system has a $U(1)$ symmetry for gapless phases. For instance, we present our Luttinger liquid parameters for the first time for the spin-1 XXZD model under the other parameters and the spin-2 XXZD model with $D = 1.5$.
      通信作者: 李生好, shanshui510@163.com
    • 基金项目: 国家自然科学基金(批准号: 11104362)、重庆市基础科学与前沿技术研究专项(批准号: cstc2018jcyjAX0812)、重庆市教委科学技术研究项目(批准号: KJQN201801212)和陕西省自然科学基金(批准号: 2019JM-017)资助的课题.
      Corresponding author: Li Sheng-Hao, shanshui510@163.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11104362), the Research Program of Basic Research and Frontier Technology of Chongqing, China (Grant No. cstc2018jcyjAX0812), the Science and Technology Research Program of Chongqing Municipal Education Commission, China (Grant No. KJQN201801212), and the Natural Science Foundation of Shanxi Province, China (Grant No. 2019JM-017).
    [1]

    Sachdev S 1999 Quantum Phase Transitions (Cambridge: Cambridge University)pp3—5

    [2]

    White S R 1992 Phys. Rev. Lett. 69 2863Google Scholar

    [3]

    Landau D P, Binder K 2011 A Guide to Monte-Carlo Simulatios in Statistical Physics (Cambridge: Cambridge University)pp70-73

    [4]

    Evenbly G, Vidal G 2009 Phys. Rev. B 79 144108.

    [5]

    Jordan J, Orus R, Vidal G, Verstraete F, Cirac J I 2008 Phys. Rev. Lett. 101 250602Google Scholar

    [6]

    Singh S, Zhou H Q, Vidal G 2010 New J. Phys. 12 033029Google Scholar

    [7]

    Jiang H C, Weng Z Y, Xiang T 2008 Phys. Rev. Lett. 101 090603Google Scholar

    [8]

    Czarnik P, Dziarmaga J 2015 Phys. Rev. B 92 035120.

    [9]

    Li W, Ran S J, Gong S S, Zhao Y, Xi B, Ye F, Su G 2011 Phys. Rev. Lett. 106 127202Google Scholar

    [10]

    Chen B B, Chen L, Chen Z Y, Li W, Weichselbaum A 2018 Phys. Rev. X 8 031082

    [11]

    Corboz P, Czarnik P, Kapteijns G, Tagliacozzo L 2018 Phys. Rev. X 8 031031

    [12]

    Singh S, Pfeifer R N C, Vidal G 2010 Phys. Rev. A 82 050301Google Scholar

    [13]

    Haghshenas R, Sheng D N 2017 arXiv: 1711.07584v1 [cond-mat.str-el]

    [14]

    Singh S, Pfeifer R N C, Vidal G 2011 Phys. Rev. B 83 115125Google Scholar

    [15]

    Song H F, Rachel S, Hur K Le 2010 Phys. Rev. B 82 012405Google Scholar

    [16]

    Song H F, Rachel S, Flindt C,Klich I, Laflorencie N, Hur K Le 2012 Phys. Rev. B 85 035409Google Scholar

    [17]

    Yang M F 2007 Phys. Rev. B 76 180403(R).

    [18]

    Boschi C D E, Ercolessi E, Ortolani F, Roncaglia M 2003 Eur. Phys. J. B 35 465Google Scholar

    [19]

    Pollmann F, Mukerjee S, Turner A, Moore J E 2009 Phys. Rev. Lett. 102 255701Google Scholar

    [20]

    苏耀恒, 陈爱民, 王洪雷, 相春环 2017 物理学报 66 120301Google Scholar

    Su Y H, Chen A M, Wang H L, Xiang C H 2017 Acta Phys. Sin. 66 120301Google Scholar

    [21]

    Yang C N, Yang C P 1966 Phys. Rev. 150 321Google Scholar

    [22]

    Chen W, Hida K, Sanctuary B C 2003 Phys. Rev. B 67 104401Google Scholar

    [23]

    Kjall J A, Zaletel M P, Mong R S K, Bardarson J H, Pollmann F 2013 Phys. Rev. B 87 235106Google Scholar

  • 图 1  (i)三指标张量${ {\varGamma}}$, 奇异值矩阵${ {\lambda}}$和粒子数n; (ii)具有$U(1)$对称的iMPS表示

    Fig. 1.  (a) Three index tensor ${ {\varGamma}}$, singular value matrix ${ {\lambda}}$ and particle number n; (b) An U(1) symmetric iMPS representation.

    图 2  更新具有$U(1)$对称的MPS的过程(a)把U门作用在具有$U(1)$对称的MPS上; (b)吸收U门缩并(a)中的张量使之成为一个两指标张量${ {\varTheta}}$;(c)对张量${ {\varTheta}}$进行奇异值分解(SVD), 得到新的张量X, Y$\tilde{{ {\lambda}}_{\rm B}}$, 同时得到新的粒子数$\tilde{n_r}$;(d)插入逆矩阵, 还原原来的原胞结构; (e)得到更新的张量$\tilde{ {\varGamma}}^{\rm A}$, $\tilde{ {\varGamma}}^{\rm B}$$\tilde{ {\lambda}}^{\rm B}$及粒子数$\tilde{n_r}$

    Fig. 2.  The process of update the U(1) symmetric MPS (a) applied gate U on the U(1) symmetric MPS, then contract the tensor network (a) into a single tensor ${ {\varTheta}}$. We compute the singular value decomposition of tensor ${ {\varTheta}}$, and get the new tensor X, Y$\tilde{{ {\lambda}}_{\rm B}}$ and particle number $\tilde{n_r}$ as in (c). (d) Insert inverse matrix and restore the original tensor structure, we obtain the new tensor $\tilde{ {\varGamma}}^{\rm A}$, $\tilde{ {\varGamma}}^{\rm B}$, $\tilde{ {\lambda}}^{\rm B}$ and particle number $\tilde{n_r}$ as in (e).

    图 3  在不同控制参量条件下, 自旋$S = 1/2$的XXZD模型的关联长度$\xi$和涨落F是截断维数$\chi$的函数. 其中, 参数$D = 0$

    Fig. 3.  Correlation length $\xi$ and fluctuation F of spin $S = 1/2$ XXZD model as a function of the truncation dimension $\chi$ for various parameters $\varDelta$. Here, fixed parameter $D = 0$

    图 4  在不同控制参量条件下, 自旋$S = 1$的XXZD模型的关联长度$\xi$和涨落F是截断维数$\chi$的函数. 其中, 各向异性参数$\varDelta = -0.5$

    Fig. 4.  Correlation length $\xi$ and fluctuation F of spin $S = 1$ XXZD model as a function of the truncation dimension $\chi$ for various parameters D. Here, fixed anisotropic parameter $\varDelta = -0.5$

    图 5  在不同控制参量条件下, 自旋$S = 2$的XXZD模型的关联长度$\xi$和涨落F是截断维数$\chi$的函数. 其中, 参数$D = 1.5$

    Fig. 5.  Correlation length $\xi$ and fluctuation F of spin $S = 2$ XXZD model as a function of the truncation dimension $\chi$ for various parameters $\varDelta$. Here, fixed parameter $D = 1.5$

    表 1  自旋S = 1/2的XXZD模型在临界区的Luttinger液体参数K, 其中参数D = 0

    Table 1.  Estimates for Luttinger liquid parameter K in the critical phase of spin S = 1/2 XXZD model with the parameter D = 0.

    Δ 0 0.25 0.5 0.75 1
    K[精确] 1.0 0.8614... 0.75 0.6493... 0.5
    K[数值] 0.999 0.856 0.7529 0.6457 0.5198
    相对误差 0.001 0.0063 0.0039 0.0053 0.0396
    下载: 导出CSV

    表 2  自旋$S = 1$的XXZD模型在临界区的Luttinger液体参数K, 固定参数$\varDelta = -0.5$

    Table 2.  Estimates for Luttinger liquid parameter K in the critical phase of spin $S = 1$ XXZD model with the parameter $\varDelta = -0.5$.

    D –0.3 0 0.3 0.5 0.6
    K[数值] 3.3679 3.1275 2.6834 2.4126 2.2745
    下载: 导出CSV

    表 3  自旋S = 2的XXZD模型在临界区的Luttinger液体参数K, 固定参数$D = 1.5$

    Table 3.  Estimates for Luttinger liquid parameter K in the critical phase of spin S = 2 XXZD model with the parameter $D = 1.5$.

    Δ 0.4 0.8 1 1.2 1.6
    K[数值] 1.652 2.5227 2.4096 2.3546 2.1107
    下载: 导出CSV
  • [1]

    Sachdev S 1999 Quantum Phase Transitions (Cambridge: Cambridge University)pp3—5

    [2]

    White S R 1992 Phys. Rev. Lett. 69 2863Google Scholar

    [3]

    Landau D P, Binder K 2011 A Guide to Monte-Carlo Simulatios in Statistical Physics (Cambridge: Cambridge University)pp70-73

    [4]

    Evenbly G, Vidal G 2009 Phys. Rev. B 79 144108.

    [5]

    Jordan J, Orus R, Vidal G, Verstraete F, Cirac J I 2008 Phys. Rev. Lett. 101 250602Google Scholar

    [6]

    Singh S, Zhou H Q, Vidal G 2010 New J. Phys. 12 033029Google Scholar

    [7]

    Jiang H C, Weng Z Y, Xiang T 2008 Phys. Rev. Lett. 101 090603Google Scholar

    [8]

    Czarnik P, Dziarmaga J 2015 Phys. Rev. B 92 035120.

    [9]

    Li W, Ran S J, Gong S S, Zhao Y, Xi B, Ye F, Su G 2011 Phys. Rev. Lett. 106 127202Google Scholar

    [10]

    Chen B B, Chen L, Chen Z Y, Li W, Weichselbaum A 2018 Phys. Rev. X 8 031082

    [11]

    Corboz P, Czarnik P, Kapteijns G, Tagliacozzo L 2018 Phys. Rev. X 8 031031

    [12]

    Singh S, Pfeifer R N C, Vidal G 2010 Phys. Rev. A 82 050301Google Scholar

    [13]

    Haghshenas R, Sheng D N 2017 arXiv: 1711.07584v1 [cond-mat.str-el]

    [14]

    Singh S, Pfeifer R N C, Vidal G 2011 Phys. Rev. B 83 115125Google Scholar

    [15]

    Song H F, Rachel S, Hur K Le 2010 Phys. Rev. B 82 012405Google Scholar

    [16]

    Song H F, Rachel S, Flindt C,Klich I, Laflorencie N, Hur K Le 2012 Phys. Rev. B 85 035409Google Scholar

    [17]

    Yang M F 2007 Phys. Rev. B 76 180403(R).

    [18]

    Boschi C D E, Ercolessi E, Ortolani F, Roncaglia M 2003 Eur. Phys. J. B 35 465Google Scholar

    [19]

    Pollmann F, Mukerjee S, Turner A, Moore J E 2009 Phys. Rev. Lett. 102 255701Google Scholar

    [20]

    苏耀恒, 陈爱民, 王洪雷, 相春环 2017 物理学报 66 120301Google Scholar

    Su Y H, Chen A M, Wang H L, Xiang C H 2017 Acta Phys. Sin. 66 120301Google Scholar

    [21]

    Yang C N, Yang C P 1966 Phys. Rev. 150 321Google Scholar

    [22]

    Chen W, Hida K, Sanctuary B C 2003 Phys. Rev. B 67 104401Google Scholar

    [23]

    Kjall J A, Zaletel M P, Mong R S K, Bardarson J H, Pollmann F 2013 Phys. Rev. B 87 235106Google Scholar

  • [1] 赖红, 万林春. 基于矩阵乘积压缩态的动态可扩展秘密共享方案. 物理学报, 2024, 73(18): 180302. doi: 10.7498/aps.73.20240191
    [2] 李新月, 祁娟娟, 赵敦, 刘伍明. 自旋-轨道耦合二分量玻色-爱因斯坦凝聚系统的孤子解. 物理学报, 2023, 72(10): 106701. doi: 10.7498/aps.72.20222319
    [3] 陈若凡. 时间演化矩阵乘积算符方法及其在量子开放系统中的应用. 物理学报, 2023, 72(12): 120201. doi: 10.7498/aps.72.20222267
    [4] 焦荣珍, 唐少杰, 张弨. 诱惑态量子密钥分配系统中统计涨落的研究. 物理学报, 2012, 61(5): 050302. doi: 10.7498/aps.61.050302
    [5] 卢道明. 三参数双模压缩粒子数态的量子特性. 物理学报, 2012, 61(21): 210302. doi: 10.7498/aps.61.210302
    [6] 田立新, 贺莹环, 黄益. 一种新型二分网络类局域世界演化模型. 物理学报, 2012, 61(22): 228903. doi: 10.7498/aps.61.228903
    [7] 万宝惠, 张鹏, 张晶, 狄增如, 樊瑛. 二分网上的靴襻渗流. 物理学报, 2012, 61(16): 166402. doi: 10.7498/aps.61.166402
    [8] 肖海林, 欧阳缮, 谢武. 量子Turbo乘积码. 物理学报, 2011, 60(2): 020301. doi: 10.7498/aps.60.020301
    [9] 周磊, 支蓉, 冯爱霞, 龚志强. 基于二分图的温度网络拓扑性质研究. 物理学报, 2010, 59(9): 6689-6696. doi: 10.7498/aps.59.6689
    [10] 阮 文, 雷敏生, 嵇英华, 谢安东. 热克尔态下介观LC电路的量子涨落. 物理学报, 2005, 54(5): 2291-2295. doi: 10.7498/aps.54.2291
    [11] 张志勇, 王太宏. 用散粒噪声测量碳纳米管中Luttinger参数. 物理学报, 2004, 53(3): 942-946. doi: 10.7498/aps.53.942
    [12] 张海燕, GNgele, 马红孺. 二分量带电胶体悬浮系统的等效硬球模型. 物理学报, 2002, 51(8): 1892-1896. doi: 10.7498/aps.51.1892
    [13] 曹天德, 黄清龙. 二分量高温超导机理. 物理学报, 2002, 51(7): 1600-1603. doi: 10.7498/aps.51.1600
    [14] 二分量胶体悬浮系统的短时间动力学. 物理学报, 2001, 50(9): 1810-1817. doi: 10.7498/aps.50.1810
    [15] 王继锁, 孙长勇. 压缩真空态下介观电路的量子涨落. 物理学报, 1997, 46(10): 2007-2009. doi: 10.7498/aps.46.2007
    [16] 谭维翰, 刘仁红. 二分岔理论的抛物线近似. 物理学报, 1990, 39(7): 35-39. doi: 10.7498/aps.39.35-2
    [17] 崔世民, 蔡建华. 二维费密液体理论(Ⅱ)——3He朗道参数的计算. 物理学报, 1990, 39(4): 572-579. doi: 10.7498/aps.39.572
    [18] 周玉魁, 云国宏. 最一般的超矩阵量子非线性Schr?dinger模型的本征态研究. 物理学报, 1989, 38(4): 648-652. doi: 10.7498/aps.38.648
    [19] 文振翼. 计算酉群生成元乘积矩阵元的一个新途径. 物理学报, 1983, 32(9): 1149-1158. doi: 10.7498/aps.32.1149
    [20] 郑兆勃. 无限次微扰理论的分块矩阵法证明. 物理学报, 1981, 30(7): 866-877. doi: 10.7498/aps.30.866
计量
  • 文章访问数:  7517
  • PDF下载量:  76
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-17
  • 修回日期:  2019-05-31
  • 上网日期:  2019-08-01
  • 刊出日期:  2019-08-20

/

返回文章
返回