格点规范场与费米子物质场耦合问题研究取得进展。格点规范场和费米子物质场之间的相互作用描述了诸多基本物理规律，涉及的范围涵盖从高能物理中的夸克禁闭到凝聚态物理中的量子自旋液体等广泛的领域。其中，(2 1) 维紧致u(1) 格点规范场与费米子物质场耦合问题格外受到关注，这是因为在凝聚态物理系统中，这个问题描述了涌现规范场(emergent gauge field)和分数化任意子激发的耦合，被认为是很多二维强关联电子体系(如铜基超导体和代数量子自旋液体)的低能有效理论描述;在高能物理中，它对应于3维紧致量子电动力学与无质量费米子耦合的问题，耦合的结果能否产生去禁闭相仍然是一个尚未解决的难题。
从凝聚态物理和高能物理的角度来看，对于规范场和物质场耦合的问题，已有的解析研究都依赖大nf微扰展开技术(这里的nf是费米子味道数)，研究发现去禁闭相在大nf极限下可以存在，而在很小nf的时候可能不稳定。但是，众所周知大nf展开技术推广到小nf(比如与实际密切相关的自旋1/2 电子 nf=2)，其可靠性急剧下降，因为1/nf在小nf时不再是小量。故而该问题的解决需要找到一个在小nf下仍然能工作的方法。目前为止还没有找到这样的可靠方法，这就是为什么近十几年来，该研究方向鲜有进展的原因。
几年以来，中国科学院物理研究所/北京凝聚态物理国家研究中心凝聚态理论与材料计算实验室研究员孟子杨与合作者，包括香港科技大学许霄琰(物理所2017年毕业生)、复旦大学副教授戚扬、密歇根大学教授孙锴、维尔兹堡大学教授fakher assaad等人运用行列式量子蒙特卡洛研究费米子-玻色子耦合格点模型，不断改进数值计算方法和条件，取得了系列性成果。最近，许霄琰、戚扬、中国科学院大学卡弗里理论科学研究所助理教授张龙、加州大学圣塔芭芭拉分校教授许岑珂、fakher assaad和孟子杨组成的国际研究团队(下面简称该团队)，寻求从大规模数值计算的角度研究规范场和物质场耦合的问题，取得了突破。他们发现 (2 1) 维紧致 u(1) 格点规范场与费米子物质场耦合问题可以通过无符号问题的量子蒙特卡洛方法严格求解，并通过大规模数值计算，得到了该问题初步的相图，取得包括 nf=2 时u(1)规范场去禁闭相仍然稳定存在这样的重要结果。相关工作发表在最近一期的physics review x 上。
物理问题的模型和相图如图1所示，该团队计算得到了在不同数量的费米子味道数 nf 和不同强度的规范涨落强度j之下的完整相图。发现当规范涨落很强时(j大时)u(1)规范场禁闭，诱导出几种不同的对称性破缺相，例如反铁磁相(antiferromagnetic phase, afm)和价键固体相(valence bond solid, vbs);而当规范涨落较弱时(j小时)，系统可以存在稳定的去禁闭相(u(1) deconfined phase, u1d，即代数量子自旋液体)。从去禁闭相到禁闭相的量子相变过程，可以通过蒙特卡洛计算中得到的物理学可观测量来确定，如图2所示。u1d相中，规范场涨落传递的费米子相互作用，系统展现出奇异的性质，比如两自旋关联和四自旋关联具有同样幂律衰减，也就是说不同数目费米子算符组成的关联函数其实具有相同的标度维度，而标度维度随着费米子味道数目(nf)的变化，最终可以和 1/nf 在大 nf 极限下的解析表达式取得一致。
这些结果说明大规模量子蒙特卡洛计算揭示了即使在 nf=2(对应自旋1/2 的相互作用电子系统)这样与实际物理体系密切相关的情况下，去禁闭相仍然稳定存在这样的重要结论。这种去禁闭相，具有普遍幂律关联，就是凝聚态物理学强关联电子理论中常常出现的代数自旋液体，作为涌现规范场和物质场耦合的结果，它是量子物质科学新范式的代表性系统。同时，这样的结果回答了3维紧致量子电动力学与无质量费米子耦合的情况下，可以存在稳定的去禁闭相这个高能物理学关心的问题。这项研究开辟了很多新的方向，比如去禁闭相和禁闭相之间的量子相变，就是超越朗道-金兹堡-威耳逊理论框架所能描述的相变实例，相关的研究正在进行中。
这项工作得到科技部重点研发计划(2016yfa0300502，2018yfa0305802)、中科院先导项目(xdb28000000，xdpb08-4, xdb28040200)、国家自然科学基金委项目(11574359，11674370)以及松山湖材料实验室、香港特别行政区hkrgc 项目hkust3/crf/13g 和 c6026-16w、德国研究基金research unit for1807 和 mercator follow 等项目的支持。量子蒙特卡洛模拟所需的大规模的并行计算在物理所量子模拟科学中心和国家超算天津中心天河1号平台和国家超算广州中心天河2号平台上完成，计算过程中得到国家超算天津中心孟祥飞、工程师菅晓东等人，国家超算广州中心应用推广部部长王栋、工程师崔颖妍等人的支持。
图1. a. (2 1) 维紧致u(1) 格点规范场与费米子物质场耦合问题的相图，横轴j控制规范场涨落，纵轴nf 是费米子的味道数。u1d 代表u(1) 规范场去禁闭相，即代数量子自旋液体)，afm 代表规范场禁闭的反铁磁相，vbs代表规范场禁闭的价键固体相。b. 微观模型在2维正方晶格上实现，蓝色格点代表费米子场，黄色圆圈代表规范场，耦合的结果是规范场倾向于形成 pi-flux，一如黄色虚线的箭头方向。c. 在去禁闭相中，费米子格林函数需要考虑路径上的规范场才是规范不变的物理学可观测量。
图2. 在费米子味道数 n=2 的时候，运用关联比例(correlation ratio r_afm)确定从u(1) 规范场去禁闭相(代数量子自旋液体)到紧闭的反铁磁长程序相的量子相变点，其热力学极限外推的位置在 jc=1.6(2)。
图3. 代数量子自旋液体具有的特殊性质，其两自旋关联和四自旋关联具有相同的标度维度(scaling dimension，delta_s = delta_d)，它们只是费米子味道数目nf的函数，当味道数目无穷大时，有1/nf 展开的解析表达式可以比对。
