具有三體相互作用的一維硬核玻色鏈的粒子輸運

任　杰，吳銀忠，2，許雪芬

（1.常熟理工學院物理與電子工程學院，江蘇常熟215500；2.蘇州科技學院物理系，江蘇蘇州215009；3.無錫職業技術學院數理部，江蘇無錫214121）

具有三體相互作用的一維硬核玻色鏈的粒子輸運

任杰1，吳銀忠1，2，許雪芬3

（1.常熟理工學院物理與電子工程學院，江蘇常熟215500；2.蘇州科技學院物理系，江蘇蘇州215009；3.無錫職業技術學院數理部，江蘇無錫214121）

摘要：利用密度矩陣重正化群的方法，研究了具有三體相互作用的一維硬核玻色鏈的粒子輸運.研究表明，在弱相互作用下，系統是彈道輸運，其彈道傳輸速度依賴于系統的粒子數，并且彈道傳輸速度隨三體相互作用增大而減小.另研究了系統處于2/3填充時基態的淬火行為.研究表明系統在無能隙相時，其輸運仍然是彈道輸運.這一有趣現象可以幫助我們探測系統的量子相變.

關鍵詞：輸運；玻色子；量子相變

近年來，對強關聯體系中非平衡態動力學行為的研究一直是理論物理和實驗物理的重要方向之一.科學家主要關注強關聯系統在非平衡態時的熱化及粒子輸運等問題.眾所周知，在經典物理中粒子的輸運主要分為彈道輸運和擴散輸運，量子系統比經典物理有著更豐富的物理現象.最近，基于實驗技術的發展，科學家已實現用冷原子系統來模擬低維量子體系.因此研究這些低維量子體系中的粒子輸運非常有意義［1］.本課題組研究了具有三體相互作用的一維硬核玻色子模型中的粒子輸運問題，該系統的哈密頓量可以寫成，其中是格點i上波色子產生（湮滅）算符，n是局域密度算符，t為交換相互作用，這里取t=1. W為三體相互作用強度. L為鏈長.

眾所周知，數值對角化系統哈密頓量是研究量子多體系統的重要方法，其困難之處在于系統的希爾伯特空間維度隨著系統的尺寸呈指數增長，鏈長為L的一維硬核玻色系統其希爾伯特空間維度大小為2L，系統的哈密頓量是2L*2L矩陣.而即便利用超級計算機數值對角化哈密頓量，能夠處理的系統鏈長上限L也不會超過40.為了避免系統邊界的影響，我們所需研究的系統鏈長必須達到百的量級.因而我們采用含時密度矩陣重正化群數值方法［2-3］，它可使我們數值模擬鏈長L達到200的系統，是研究強關聯系統的重要方法之一.

圖1　不同三體相互作用強度下，密度算符的預期值（ni）隨時間的變化關系（a）W=0，（b）W=1，（c）W=3，白色虛線是系統的最大幾率速度2t［4］

我們數值模擬了系統的局域密度算符的預期值在不同三體相互作用W下隨時間的變化關系（見圖1）.粒子的擴散行為表明在弱相互作用下，系統中的波色子傳輸行為是一種彈道行為.同時我們利用波色子的填充情況得到其彈道傳輸速度.研究表明彈道傳輸速度隨著三體相互作用增強而減小，同時還依賴于系統中波色子數目（見圖2）.

圖2　不同玻色子數目條件下的輸運速度與三體相互作用之間的關系圖（插圖：粒子輸運速度V隨玻色子數目N在W=1，2下的關系變化圖［4］）

此外，我們研究了2/3填充（N/L=2/3）情況下基態的淬火行為，得到當系統初態處于無能隙相時，系統中波色子的傳輸行為仍然是一種彈道行為.同時研究了粒子輸運速度與三體相互作用的變化關系，發現其隨W先增大后減小（見圖3）.這一有趣現象可以幫助我們探測系統中的量子相變.系統的尺度效應如圖3所示，從圖3可以得到系統的超流相與固態相的相變點發生在WC= 2.74±0.1.這與文獻［5］的結論吻合得非常好.

圖3　粒子輸運速度V隨三體作用W變化關系圖（插圖：粒子輸運速度V的最大值隨系統長度L的關系變化圖［4］）

參考文獻：

［1］POLKOVNIKOV A, SENGUPTA K, SILVA A, et al. Colloquium：Nonequilibrium dynamics of closed interacting quantumsystems ［J］. Rev Mod Phys, 2011, 83：863.

［2］SCHOLLW?CK U. The density-matrix renormalization group［J］. Rev Mod Phys, 2005, 77：259.

［3］REN J, WU Y Z, ZHU S Q. Study of dynamics in spin-1/2 chain using adaptive time-dependent density-matrix renormalizationgroup method［J］. Int J Theor Phys, 2013, 3167：52 .

［4］REN J, WU Y Z, XU X F. Expansion dynamics in a one-dimensional hard-core boson model with three-body interactions［J］. Scientific reports 5, 2015：14743.

［5］CAPOGROSSO-SANSONE B, WESSEL S, BüCHLER H P, et al. Phase diagram of one-dimensional hard-core bosons with three-body interactions［J］. Phys Rev B, 2009, 79：020503（R）.

Expansion Dynamics in a One-dimensional Hard-core Boson Model with Three-body Interactions

REN Jie1, WU Yinzhong1,2, XU Xuefen3
1.School of Physics and Electronic Engineering, Changshu Institute of Technology, Changshu 215500, China；2. Department of Physics, Suzhou University of Science and Technology, Suzhou 215009, China；3. Department of Mathematics and Physics, Wuxi Institute of Technology, Wuxi 214121, China）

Abstract：Using the adaptive time-dependent density matrix renormalization group method, the authors of this paper numerically investigate the expansion dynamics of bosons in a one-dimensional hard-core boson model with three-body interactions. It is found that the bosons expand ballistically with weak interaction. It is also shown that the expansion velocity V is dependent on the number of bosons. As a prominent result, the expansion velocity decreases with the enhancement of three-body interaction. Besides, a further study is made of the dynamics of the system, which quenches from the ground state with two-thirds filling, and the results indicate the expansion is also ballistic in the gapless phase regime, which can help us detect the phase transition in the system.

Key words：expansion；boson；quantum phase transition

中圖分類號：O469

文獻標識碼：A

文章編號：1008-2794（2016）02-0123-02

收稿日期：2015-12-26

基金項目：國家自然科學基金“應用重整化群研究強關聯系統中的量子糾纏和調控”（11104021）

通信作者：任杰，副教授，研究方向：凝聚態物理、量子信息，E-mail：jren@cslg.cn.