← 返回 PaperDaily
大模型与智能体
清华DMRG新发现:手性自旋液体一场就倒?
手性自旋液体向来是强关联物理的“圣杯”级物态,但外界条件一变化它会不会垮?清华大学姚宏团队用大规模DMRG数值实验,把位移场这个“旋钮”拧开:手性自旋液体只撑到φ≈0.02π,随后跳出个带长程标量手性的反铁磁姐妹相,中间还可能藏着一个连续相变。一篇paper把相图玩出三种新花样,值得一读。
龙哥读论文
发布于 2026-08-19 00:20:09
阅读 3
查看原文
原论文信息如下:
手性量子物态的新舞台:moiré体系中的位移场调控
手性量子物态的新舞台:moiré体系中的位移场调控
先给不熟悉moiré体系的朋友补个背景:把两层过渡金属硫族化合物(TMD,transition metal dichalcogenide)叠在一起,稍微转一个角度,就会形成周期性的莫尔条纹——这套“转魔方”的技术,如今已经是凝聚态物理学家变出各种新奇量子相的标配舞台。在这个舞台上,转动角度、载流子浓度、垂直电场(即位移场)都是随时可调的“旋钮”。过去几年里,科学家们在这类moiré体系中陆续找到了关联绝缘体、分数量子反常霍尔态、复合费米液体甚至超导,简直像在玩一个“量子相的宝可梦收集游戏”。
本篇要聊的论文,来自清华大学高等研究院的姚宏团队,主角就是上面提到的“位移场”这个旋钮。论文用大规模密度矩阵重整化群(DMRG)数值模拟,系统刻画了半填充moiré Hubbard模型在位移场影响下的基态相图。结果相当有意思:本来被寄予厚望的手性自旋液体(CSL,chiral spin liquid) 在位移场面前相当“脆”,稍微一加场就撑不住了;但它倒下去之后,又扶起来一个带着长程标量手性的反铁磁姐妹相;更妙的是,这两个相之间还可能出现连续量子相变——这在强关联物理里可是个香饽饽。
一句话概括:加个电场,手性自旋液体就变身了? 没错,这篇论文说的就是这么回事。下面咱们一层层拆开看。
从Hubbard模型到手性自旋液体:DMRG揭示的丰富相图
先看模型。研究者考虑的是一个定义在三角晶格上的最小moiré Hubbard模型,只保留最近邻跳跃和在位Hubbard排斥项,哈密顿量写出来长这样:
公式中t是跳跃振幅,U是在位库仑排斥强度,σ=↑,↓标记自旋-谷锁定指标(spin-valley locked index),c†和n分别是费米子产生算符和密度算符。最特别的是跳跃项里的相位因子e^{iσφ_ij}——这正是位移场的“指纹”。
这个相位因子的物理起源值得展开:微观上位移场本身跟自旋无关,它只是打破了两层之间的镜面对称。但在投影到自旋-谷锁定能带之后,层指标跟自旋指标耦合在一起,于是位移场就变成了一种有效的自旋轨道耦合(SOC,spin-orbit coupling) 。更具体地说,它把体系点群对称性从D6h降到C3v,自旋对称性从SU(2)降到U(1)。从模型上看,自旋向上的电子感受到一组交错通量,自旋向下的电子则因为时间反演对称而感受到相反的通量——也就是说,每个原胞的总磁通为零,但两个不等价三角格点上的通量分别为±3φ。
图1给出了本文最核心的结果之一:(a)是DMRG在YC-48×4圆柱上得到的强关联区基态相图,黑点是模拟数据点;(b)是自旋向上电子的交错跳跃通量示意图;(c)是自旋向上和自旋向下在φ=0和φ=0.04π时的费米面。可以看到,随位移场从零增大,相图发生了显著变化,这块后面细说。
数值方法方面,研究者用的是密度矩阵重整化群(DMRG,density matrix renormalization group)——这在强关联费米子体系中堪称“数值显微镜”。众所周知,三角晶格Hubbard模型是出了名的棘手:量子涨落强、阻挫严重,传统平均场几乎没法信任。DMRG通过在有限宽圆柱上做高精度变分,能可靠地刻画这类强关联基态。本文用了YC圆柱几何,主系统圆周L_y=4、长度L_x最大到48,保留U(1)_charge × U(1)_spin对称性,最大键维D=10000,截断误差在10⁻⁵到10⁻⁷之间。这精度放在几年以前是想都不敢想的。
先复习一下什么叫手性自旋液体(CSL)。传统的自旋液体,简单说就是自旋系统即便在零温下也不形成任何磁序,自旋就像液体一样无序地“流动”着。而“手性自旋液体”更进一步:它虽然不形成常规磁序,却自发破坏了时间反演对称性,在微观三重态上形成有规则的“涡旋”,也就是长程标量手性序。更重要的是,CSL具有拓扑序,体内有分数化激发,边界上有手性边缘态——换句话说,这是一种真正“非平凡”的量子物态。
在零位移场下,DMRG早就确认了CSL存在于金属和120°反铁磁之间的一个窄窗口里(大约9 ≲ U/t ≲ 10.5),这也是近年来三角格点Hubbard模型研究最受关注的结论之一。关键问题是:加上位移场之后,这个CSL还活得下去吗?
论文给出的答案相当干脆:CSL只撑到φ≈0.02π 。再大的位移场就凉了。要刻画CSL,研究者测量了标量手性关联、单粒子格林函数和自旋关联三种量。标量手性的定义如下:
其中(i, i+e₁, i+e₂)是同一个三角格点上的三个格点,S_i是格点i的自旋算符。这个量刻画的是三个自旋在某个原胞上形成的“手性”程度——如果自旋构型存在某种有序的“拧转”,这个关联函数在长距离下就不衰减。
图2展示了沿某条参数线的CSL证据:(a)标量手性关联函数在长距离下保持有限,说明存在长程手性序;(b)单粒子格林函数指数衰减,说明体系是单粒子有能隙的绝缘体;(c)自旋关联函数快速衰减到零,排除磁有序。三者合在一起,就是CSL的典型“体检报告”。不过也能清楚看到,位移场增大时,手性关联的抑制非常明显——CSL确实在变弱。
位移场如何重塑量子磁性:从CSL到手性反铁磁
CSL倒下了,那总该有个“接盘侠”。论文的重头戏之一,就是在强耦合区发现了一个全新的手性反铁磁(chiral AFM)相。
先把模型做强耦合展开。当U/t比较大时,通过标准的二阶微扰理论,可以把Hubbard模型映射成一个自旋模型,得到位移场可调的XXZ各向异性加上Dzyaloshinskii-Moriya(DM)相互作用。这就是本文分析的起点。
值得注意的是,φ=0时这个模型还原为各向同性海森堡模型,基态就是120°反铁磁,自旋取向平面可任意选择,矢量手性κ也有两个简并方向。一旦φ≠0,SU(2)旋转对称性和空间反演对称性都被打破,DM项开始主导,它自然会挑选面内的120°反铁磁序,并且固定矢量手性的方向。按这个思路,位移场本质上只是给原来的反铁磁相“挑了挑衣服”——确定自旋取向平面和手性方向而已。如果仅此而已,那这篇文章的戏剧性就得大打折扣了。
真正的惊喜出现在标量手性关联上。图3(d)给出了U/t=12时不同φ下的标量手性关联函数:φ=0时关联快速衰减,跟普通共面反铁磁的预期一致;但φ≠0时,关联在长距离下保持有限,出现长程标量手性序!这就意味着,有限位移场下的反铁磁相不仅仅是“选了个平面”,它还自发打破了复合对称性T~=e^{-iπS^z_tot}T(即全局自旋z旋转与时间反演的复合操作),成为一个与常规120°反铁磁截然不同的新物态——手性反铁磁相。
图3还对磁结构本身做了全面“体检”。(a)和(b)分别展示自旋结构因子S^{+-}(q)和S^{zz}(q),(c)是两者峰值随φ的变化。可以看到,随着φ增加,面内通道S^{+-}的峰越来越强、越来越尖锐,而面外通道S^{zz}的峰被抹平、削减。这正好和强耦合展开的预期一致:DM项在加强面内120°反铁磁序。所以最终图像是:位移场不仅强化了面内磁序,还额外地给整个磁结构“拧上”了一个标量手性维度。
这里有两个可能的来源:要么是微弱的面外非共面磁结构在半经典层面贡献了标量手性;要么是即使自旋还在面内共面排列,纯粹的量子涨落也能产生标量手性关联。论文没有强行区分两种机制,但从对称性角度看,不管哪种来源,这个相都确实打破了T~对称性,是货真价实的新相。
连续相变的蛛丝马迹:分数化量子临界性的可能证据
既然CSL和手性反铁磁在手性这个维度上是“同一个阵营”的,那它们之间的相变会是什么性质?这个问题恐怕是全文最诱人的部分。
为什么说它诱人?因为在φ=0的常规场景里,CSL到120°反铁磁的相变在理论上非常难处理:CSL有拓扑序,反铁磁没有拓扑序,两者之间的连续相变往往需要精细调节,一般不容易实现。但如今有限位移场下的场景不一样了——CSL和手性反铁磁两侧都保持标量手性序,这等于说两个相在“手性”这个序参量上是连续的,拓扑序的消失和磁序的出现可以在同一个连续相变里完成,这在理论上顺理成章得多。
图S1展示了φ=0.005π下穿过相变点的数值证据:(a)单粒子格林函数在不同U/t下从代数衰减变为指数衰减;(b)标量手性关联在相变两侧都保持长程;(c)自旋结构因子S^{+-}(q)在强耦合侧出现尖锐磁峰;(d)磁峰峰值随U/t连续增大。综合来看,在数值分辨率范围内没有看到任何一级相变的跳变特征,更像是一个连续相变。
论文提出了两个可能的连续相变场景。第一个是自旋子凝聚机制 (spinon condensation):CSL里的分数化激发——自旋子——在临界点处凝聚,同时产生磁序并摧毁拓扑序,但手性序保持不变。这个场景和ν=1/2玻色子Laughlin态到超流态的相变颇为相似,可能由N_f=2的QED₃-Chern-Simons理论描述。第二个是常规磁性相变 ,即拓扑序一路保持到手性反铁磁相中,相变由常规磁有序驱动,归入3D-XY普适类。两种场景都很难得,未来需要更仔细的数值分析来区分它们。
如果第一种场景成立,意味着我们在moiré TMD体系里就能看到分数化量子临界性——这是强关联物理里公认的“圣杯”之一。虽然目前还只是理论预测,但至少地图已经画出来了,接下来就看实验物理学家怎么去对号入座了。
实验启示与未来展望:通向moiré TMD中的新奇量子相
相图里除了CSL和手性反铁磁这两个“主角”,还有一位“配角”也值得关注——非公度自旋密度波(IC-SDW,incommensurate spin-density wave)。在中等相互作用强度、有限位移场下,电子系统会形成一种调制周期与晶格不匹配的自旋密度波态。
图4展示了φ=0.04π时不同U/t下的行为。(a)是单粒子格林函数,U/t≤7时缓慢衰减对应金属态,U/t≥7.5时快速指数衰减对应绝缘体;(b)是自旋结构因子S(q),在U/t≥8.5时只有一个尖锐的K点附近峰,而在U/t≤8时峰劈裂成两个非公度峰;(c)显示了峰值随U/t的系统尺寸依赖,据此估计金属到IC-SDW金属的临界U/t约4.5。
IC-SDW的形成从弱耦合视角很好理解:有限位移场下自旋向上和向下的费米面之间存在近似K点的嵌套波矢,这会增强自旋密度波的不稳定性,但嵌套不完美,所以需要有限相互作用才能发展出长程序。在L_y=4圆柱上,IC-SDW绝缘相看起来插在金属和强耦合反铁磁之间;但论文自己也承认,L_y=6的初步结果显示强耦合公度反铁磁可能更占优——IC-SDW绝缘相在L_y=4上的稳定性更像是一个“激烈竞争的候选者”,稍微加一点次近邻跳跃就可能改变胜负。这个点处理得很诚实。
那实验上该怎么探测这些相?论文给了几个很具体的建议。标量手性序会破坏复合时间反演对称T~,因此可以做磁光Kerr效应 (magneto-optic Kerr effect)和热霍尔输运 (thermal Hall transport)实验。测温依赖也会有用,因为标量手性序在温度升高时会熔化。理论上估算一下:CSL能存活到的φ≈0.02π,对应到4°扭转双层WSe₂大约3 meV的层间能量差,或者约10 mV·nm⁻¹的垂直位移场——这些都在现代器件能做到的范围内,并不是遥不可及的数字。
论文还跟最近的实验对了一下账。近年来在扭转双层WSe₂里确实观察到了半填充附近的超导,而一个更近的实验在Mott绝缘区内发现了一个展宽的非磁绝缘区,随着位移场增大,它插在超导和反铁磁绝缘体之间。论文指出这个实验相序和他们在U/t≈8-9看到的金属→CSL→手性反铁磁序列惊人地相似。如果后续实验能在候选CSL区和反铁磁区都探测到标量手性,而且它们之间是连续转变,那本论文的“剧本”就真的被验证了。
最后说点个人感受。这篇论文最让人舒爽的地方在于它不飘:每个结论都给出具体数值、具体几何、具体的有限尺寸检查。CSL只活到φ≈0.02π、IC-SDW在L_y=6上的敏感性、连续相变的数值分辨率限制——这些“丑话”都说在前面。在拓扑序、分数化、量子临界性这些动不动就“突破”“颠覆”的领域里,这种诚实尤为可贵。
龙迷三问
这篇论文到底在解决什么问题? 清华大学姚宏团队用大规模DMRG研究半填充莫尔Hubbard模型,发现位移场显著重塑相图:手性自旋液体仅存活至φ≈0.02π,强耦合区出现长程标量手性的反铁磁相,中间耦合区稳定出非公度自旋密度波,CSL到手性AFM或为连续相变。
这篇工作最值得看的点是什么? 论文通过DMRG模拟成功绘制了位移场依赖的基态相图,发现位移场显著丰富三角晶格Hubbard模型的相图:快速抑制CSL相、诱导强耦合手性反铁磁相中的长程标量手性序、稳定中等耦合下的非公度自旋密度波相,并发现CSL到手性AFM的连续相变迹象。
这篇工作的边界或风险在哪里? 优点:(1) 使用大规模DMRG模拟,键维数高达10000,结果可靠;(2) 系统研究了位移场对三角晶格Hubbard模型相图的影响,填补了该领域空白;(3) 发现CSL到手性AFM的连续相变迹象,为分数化量子临界性提供了可能平台;(4) 在Ly=4和Ly=6两种圆柱几何上进行了交叉验证。缺点:(1) 有限周长效应仍需进一步评估,特别是IC-SDW绝缘体相在Ly=6上未被观察到;(2) CSL-AFM转变的连续性受数值分辨率限制,不能排除弱一级相变或窄中间相的可能性;(3) 未直接研究超导性,对实验观测的指导有限。
如果你还有哪些想要了解的,欢迎在评论区留言或者讨论~
龙哥点评 论文创新性分数: ★★★★☆
利用大规模密度矩阵重正化群(DMRG)模拟,系统研究半填充三角晶格moiré Hubbard模型在位移场作用下的基态相图,揭示手性自旋液体、手性反铁磁及非公度自旋密度波等新奇量子相。
实验合理度: ★★★☆☆
现有材料未完整覆盖数据划分、基线公平性和统计显著性,因此按中性评价处理。
学术研究价值: ★★★★☆
利用大规模密度矩阵重正化群(DMRG)模拟,系统研究半填充三角晶格moiré Hubbard模型在位移场作用下的基态相图,揭示手性自旋液体、手性反铁磁及非公度自旋密度波等新奇量子相;更关键的是问题定义是否可复用到同类任务。
稳定性: ★★★☆☆
现有材料未提供充分的极端条件、重复运行或扰动测试,稳定性暂按中性评价。
适应性以及泛化能力: ★★★☆☆
现有材料未完整展示跨数据集、跨场景或分布外实验,泛化能力仍需进一步验证。
硬件需求及成本: ★★★☆☆
DMRG模拟,键维数最高达D=10000,截断误差约10^-5至10^-7
复现难度: ★★★☆☆
现有材料未确认完整代码、配置、数据处理脚本和权重是否齐备,复现难度暂按中性评价。
产品化成熟度: ★★★☆☆
论文验证以研究实验为主,真实部署中的时延、成本、维护和异常场景仍需补充验证。
可能的问题: 现有材料尚未充分覆盖分布外泛化、部署成本、长期稳定性和失败案例。
[1] Tuo, C., & Yao, H. Chiral spin liquid and chiral antiferromagnetism in half-filled moiré Hubbard model: possible applications to twisted bilayer TMDs. arXiv:2608.14503 (2026).
[2] Pan, H., Wu, F., & Das Sarma, S. Band topology, Hubbard model, Heisenberg model, and Dzyaloshinskii-Moriya interaction in twisted bilayer WSe₂. Phys. Rev. Research 2, 033087 (2020).
[3] Zang, J., Wang, J., Cano, J., & Millis, A. J. Hartree-Fock study of the moiré Hubbard model for twisted bilayer transition metal dichalcogenides. Phys. Rev. B 104, 075150 (2021).
[4] Xia, Y., Han, Z., Watanabe, K., Taniguchi, T., Shan, J., & Mak, K. F. Superconductivity in twisted bilayer WSe₂. Nature 637, 833 (2025).
[5] Guo, Y., et al. Superconductivity in 5.0° twisted bilayer WSe₂. Nature 637, 839 (2025).
*本文仅代表个人理解及观点,不构成任何论文审核或者项目落地推荐意见,具体以相关组织评审结果为准。欢迎就论文内容交流探讨,理性发言哦~ 想了解更多原文细节的小伙伴,可以点击 "阅读原文", 查看更多原论文细节哦!
*本文仅代表个人理解及观点,不构成任何论文审核或者项目落地推荐意见,具体以相关组织评审结果为准。欢迎就论文内容交流探讨,理性发言哦~ 想了解更多原文细节的小伙伴,可以点击 "阅读原文", 查看更多原论文细节哦!