冷原子体系的量子波动性、宏观量子相干性和人工可调控性，使其成为了一个全新的量子体系，其新颖的量子态和奇异物性的研究是国际上具有前瞻性和挑战性的前沿领域。自年实现稀薄气体玻色—爱因斯坦凝聚以来，从单组分、简单相互作用的研究逐渐过渡到多组分、复杂多体效应以及自旋—轨道耦合、非厄米、强关联、无序效应等新物理的研究。文章介绍了近几年冷原子方面的研究进展，包括冷原子的相关技术，冷原子在量子精密测量、量子模拟和量子计算方面的重要工作，期望给未来的研究以新启迪。
撰文
邢益辉、李文通、刘伍明（中国科学院物理研究所)
来源
本文选自《物理》年第2期
01
引言
在微观物理学中，温度是原子做无规则热运动剧烈程度的度量标准：温度越高，原子的运动速率越快，热运动越剧烈；温度越低，原子的运动速率越慢，热运动越微弱。这种无规则热运动是难以预测和操控的，往往会在实验过程中引入误差。可想而知，如果能够尽量减弱原子热运动的剧烈程度，也就是尽可能地降低温度，就会使实验得到极大优化，这一研究领域统称为冷原子物理。
冷原子具备以下两个特征：第一，冷原子运动速率慢，相互之间碰撞少，使其易于操控的同时也降低了能谱的展宽，有助于提高测量精度；第二，冷原子的德布罗意波长很长，表现出明显的量子特征，具有很强的相干性，可以产生宏观量子效应，一个最重要的代表就是玻色—爱因斯坦凝聚态(BEC)，其中所有的玻色子都处于能量最低的基态。基于这两个优异的特性，冷原子物理在量子精密测量、量子模拟、量子计算等诸多领域具有广泛的应用。基于原子干涉的高精度绝对重力仪可以通过测量重力场的变化进行地震监测和地下水探测，更准确的原子钟为全球定位和通信系统提供了强有力的保障，使用冷原子制造的量子存储器、量子芯片将是量子计算机的基本部件，高度可控、纯净的冷原子系统为Hubbard模型、Su—Schrieffer—Heeger(SSH)模型等量子体系的模拟提供了理想的实验平台。
02
冷原子技术
年，Hnsch和Schawlow提出了激光冷却的方法[1]，基本原理是利用光子与原子之间的散射降低原子的运动速率，从而达到冷却的效果。当一个原子与相向运动的光子碰撞时，原子会将光子吸收并跃迁到激发态，由于动量守恒，原子将会以一个比之前慢的速率继续向前运动，之后在某个时刻自发跃迁回基态，随机向一个方向放出一个光子。如果大量重复这一过程，向各个方向放出的光子对原子的反冲作用会相互抵消，最终表现为原子沿原来的方向减速运动。
激光冷却可以获得极低的温度，为人们研究冷原子物理打开了一扇大门。年，Phillips等人首次实现了中性钠原子的激光冷却[2]，将钠原子的速度降低到原来的4%，相当于把温度降低到70mK。仅仅在三年后，朱棣文等人就将钠原子进一步冷却到μK[3]，达到了多普勒冷却的极限。之后，偏振梯度激光冷却、速度选择相干粒子数囚禁冷却等方法更是突破了多普勒冷却极限。时至今日，人们已经可以在实验室中获得pK级别的温度，正逐渐向绝对零度逼近，并且实现了几乎所有碱金属原子以及部分碱土和稀土原子的冷却。除了原子外，对分子[4]和离子的冷却也逐步被实现。Langin等人通过对超冷中性原子气体进行光电离得到了超冷中性等离子体[5]，温度达到50mK，打破了传统的高温等离子体的限制，使中性等离子体的研究深入到强耦合区域。甚至人们还对反物质进行了激光冷却的研究。年，Baker等人首次对一个反质子和一个反电子组成的反氢原子进行了一维的激光冷却[6]，观察到了比未冷却的反氢原子更窄的1s—2s跃迁光谱(图1)，这一结果会促进正在进行的对反氢原子的光谱和引力的研究，并为未来的反物质实验开创了新方法。
图1冷却后的反氢原子具有比未冷却的反氢原子更窄的1s—2s跃迁光谱[6]
与激光冷却同时发展起来的另一种方法是蒸发冷却，具体过程是在原子团中速率分布达到平衡时，将其中速率较大的原子剔除，剩余的原子将在碰撞中达到新的平衡，这时原子的速率总体上会减小，能够使温度进一步降低，这对于BEC的获得具有重要作用。此外，Gisbert等人在研究一维冷原子链时发现了一种由原子间交换光子导致的协同冷却机制[7]，当链足够长时能够消减原子的自发辐射，使链更稳定，分析还表明这种作用在二维情况下可能会更显著。
原子无时无刻不在运动，为了防止原子到处乱跑，需要将它俘获在势阱中。利用反向激光束之间的干涉，在空间上可以形成稳定的周期性光学势阱阵列，被俘获的原子在这些光学势阱中有序排列，类似于晶体结构，因此称为光晶格。通过设计激光束的方向、波长和势阱深度等参数可以得到各种晶格结构，并控制原子在不同格点间的跃迁。年，Endres等人利用个光镊在不到ms的时间内实现了将50多个原子逐个组装成无缺陷的一维阵列[8]。李在勋(LeeJaeHoon)等人设计的一种光学补偿变焦透镜[9]，可以产生一个光学偶极子阱用来输运冷原子，具有精度高、势阱深度恒定、移动距离远、配置容易等优点，可以对光晶格的组装方法进行补充。三维光晶格的组装比较困难，尤其是每层之间形状存在差异或堆叠方式不同的三维晶格，如扭曲双层石墨烯。因此郝磊(HaoLei)提出了一种逐层组装三维光晶格的方案[10]，首先用激光分束器获得多层二维光晶格，再用两个圆柱形薄透镜组成的光束整形器压缩z方向的激光束，使层与层之间的距离变小，就可以得到三维光晶格，并可以通过控制相位参数实现不同层之间的相对滑移，和通过旋转x—y平面上的反射镜实现不同层之间的相对扭转，这一方案在复杂晶格结构的研究中具有很好的应用前景。
通过对冷原子系综折射或吸收的光图像进行分析，就能够提取出有关原子系综的物理信息。常用的成像技术有吸收、荧光和相位对比，还有离共振离焦成像(ORDI)[11]。对于冷原子气中的离子，Gross等人提出了一种利用离子—里德伯原子相互作用诱导的吸收进行成像的方法[12]，可以以时间分辨的方式成像离子的动力学演化过程。
03
量子精密测量
量子精密测量是冷原子的重要应用之一，它的主要任务是不断提高测量的精度。通常我们有两种手段来提高测量精度，最直接的方法就是寻找最小刻度更小的“尺子”，例如最小刻度一米的尺子无法测量几厘米的长度，但最小刻度是一厘米的尺子就可以测量几厘米的长度；另一种方法是进行多次测量，利用统计规律来减小每次测量产生的误差，数学上的中心极限定理告诉我们，对同一个量进行N次独立重复的测量，得到的所有结果服从正态分布，每次测量的误差为测量值的
，称为散粒噪声极限，是经典测量方法在理论上所能达到的最高精度。
随着量子技术的不断发展，人们发现如果在测量中利用量子力学的独特性质，如相干性、纠缠性等，就可以打破经典的散粒噪声极限的限制，进一步提高测量的精度。如果让N个探测粒子的量子态相互纠缠，外界对这N个粒子的作用就会相干叠加，最终得到的误差为测量值的
，比经典的散粒噪声极限允许的精度提高了
倍，称为海森伯极限。这是由海森伯不确定原理导致的测量误差，来源于真空中的量子涨落，是无法被消除的，换言之，海森伯极限是量子测量方法在理论上所能达到的最高精度。
尽管理想的量子测量方法可以得到很高的精度，但将N个粒子纠缠起来是极其困难的，实际上难以实现。量子精密测量的研究任务就是要突破散粒噪声极限，不断逼近海森伯极限。由于冷原子的德布罗意波长很长，非常容易表现出量子性质，因此适合用来作为量子精密测量的工具。下面将从三个方面介绍近年来在逼近海森伯极限过程中的进展。
3.1使用压缩态进行量子测量
海森伯不确定关系限制了我们在同时测量多个物理量时得到的整体精度，如果放宽对其中一个物理量的精度要求，就能使另一个物理量得到更高的精度，这种量子态称为压缩态。Hosten等人基于光学腔测量时钟态铷原子的实验对噪声水平进行了计算[13]，在没有纠缠的情况下，误差可以表示为
N为原子数，通过对量子态在自旋自由度上进行“压缩”，即产生了纠缠，如图2所示，未压缩时量子态在y和z方向的不确定度相同，压缩后y方向的不确定度增大，而z方向的不确定度减小了，使用压缩态再次进行测量得到z方向的测量误差，通过计算可以得到
其中α是由原子云和激光束的形状确定的参数，由于α1，可知在有纠缠的情况下测量误差比没有纠缠的情况下减小了很多，意味着超出了经典的散粒噪声极限。
图2(a)未压缩的自旋量子态在布洛赫球上的表示，斑点代表不确定度，未压缩时在y和z方向的不确定度相同；(b)未压缩的量子态在z方向的概率分布；(c)使用未压缩态时的误差期望值；(d)压缩的自旋量子态在布洛赫球上的表示，y方向的不确定度增加，z方向的不确定度减小了；(e)蓝色为未压缩态的概率分布，红色和*色为(d)图中两个压缩态的概率分布[13]
年，Szigeti等人基于对冷原子重力仪的研究，提出BEC中固有的原子间相互作用会使得测量的精度提高[14]。体系的哈密顿量为
3.2与弱测量相结合
在量子力学中，测量过程实际上是使仪器与量子系统耦合的过程，测量的结果是将系统的量子态投影到某个本征态上的本征值，测量之后系统的波函数就坍缩到这个本征态上，使得系统的状态回不到原来的状态，也就不能进行重复测量，这是我们所熟知的强测量。年，Aharonov提出了弱测量的概念，使仪器与量子系统的耦合作用很弱，只从系统中提取
达到了海森伯极限。实验中使用了大约十万个光子，相较于经典测量方法精度提高了两个数量级。这是世界上首次在实际测量中达到海森伯极限，并且不需要使用纠缠就能实现。
图3测量精度与光子数的关系，图中点是选取了一组特定参数得到的结果，蓝线是对点的拟合，紫线是混态测量精度的边界[15]
图4三个参数的测量误差。图中实线和虚线代表最小理论误差，点表示实验测量结果，可以看出实验中三个参数同时达到了最小误差。其中曲线①为经典独立测量，②为纠缠独立测量，③为纠缠联合测量，④为控制增强顺序测量[16]
3.3多参数测量的最优方案
在实际测量时，往往需要同时对多个参数进行测量，根据海森伯不确定原理，如果两个参数是不对易的，就无法同时对它们进行精确测量，之前提到的压缩态就是舍弃其中一个参数的精度来换取另一个参数的更高精度。那么，如何同时兼顾多个参数来进行最佳的测量呢？
图5对复数α的实部和虚部使用相干态进行参数估计的均方误差，在曲线左下方的区域是被不确定原理所禁止的。黑色实线代表制约关系；红色虚线代表基于右对数导数的几何平均量子测量极限；蓝色虚线代表基于右对数导数的算术平均量子测量极限；绿色点线代表基于对称对数导数的调和平均量子测量极限[17]
04
量子模拟
得益于原子控制技术的进步，冷原子成为量子模拟的一个理想化的平台[18,19]。我们选择如下4个主题来展示冷原子的模拟能力。一是冷原子模拟自旋—轨道耦合(SOC)，也包括最近提出的质心轨道角动量和自旋间的耦合，即自旋—轨道—角动量耦合(SOAMC)。二是近几年凝聚态物理领域被广泛研究的物相，如超固态和拓扑相，体系由包括了SOC，Hubbard项的哈密顿量描述，在这些相互作用的竞争下，物理体系可以表现出奇特的相变行为。三是冷原子模拟非厄米量子系统，着重分析其不同于厄米体系的行为。四是不同于之前在平坦黎曼流形上的分析，而是在双曲流形上对BEC现象给予新的物理图像。
当然，冷原子的模拟能力远比这强大，许多有趣的物理系统并未包含在这些课题中，如模拟非线性系统[20]、耗散系统[21]等。
4.1自旋—轨道耦合
自旋—轨道耦合(SOC)即原子内部的超精细自旋通过双光子拉曼过程与原子质心动量耦合[22—24]，等价于受到SU(2)非阿贝尔规范势的作用[25]，可以产生很多有趣的物理效应[26，27]。近来，原子质心轨道角动量和超精细自旋耦合产生的新现象被发现[28]，通过两束Laguerre—Gaussian型激光，可以在旋量BEC中实现这种SOAMC[29，30]。
图6自旋为1/2的费米原子与拉曼激光作用(a)两束同向传播的携有不同轨道角动量(-l1和-l2)的高斯型激光与原子作用诱导的SOAMC；(b)能级跃迁示意图；(c)Ω0—δ平面内二维SOAMC的费米超流的多体相图，其他参数设置为l=3，EB/EF=0.5。图中包括超流态(SF)、正常态(N)、有能隙的涡旋态(V1)和无能隙的涡旋态(V2)[31]
考虑约束在x—y平面的二组分自旋1/2费米气体，与两束不同轨道角动量的拉曼激光耦合(图6(a)，(b))。通过一个幺正变换，得到极坐标系下的单粒子有效哈密顿量[31]：
由此可见，利用高自旋的SOC或者SOAMC，可以发现许多其他新奇的物理相，这使其成为一种“技术手段”，与各种模型结合，得到更丰富的物理图像，如超固体(SS)、超辐射[33]、拓扑物态等。
4.2量子相
超固体既是具有非对角长程序的超流态，也是具有对角序的固态，概念最早在固体4He中被提出[34]。作为一种新奇的物态，其实现方法一直是理论和实验追寻的目标。理论中预测可以实现超固态的方法有很多，例如在一维或二维旋量BEC系统中施加SOC[35，36]、偶极费米系统等。
文献[35]中考虑二维有SOC和软核长程相互作用的玻色气体，在Gross—Pitaevskii平均场近似下的哈密顿量可以写成：
通过虚时算法数值最小化哈密顿量得到多体基态，如图8所示，参数固定时，涡旋环流方向相同，呈现周期性的涡旋排列，此即超固态。同时相互作用强度较弱的组分构成涡旋的中心部分，被相互作用强的组分包围。当交换相互作用强度比，涡旋环流方向改变，不同的组分交换分布方式，故而此超固态称为手征超固态。当继续改变相互作用强度比，还会出现其他超固态相，如平面波超固态、驻波超固态、自旋极化的超流相等。
周期性光学势阱中的冷原子可以只占据少量的低能布洛赫带，且同一势阱中的冷原子间相互作用很大，可以很好地模拟强关联的晶格模型，如Hubbard模型[38，39]。同时Hubbard模型可以产生许多新奇的拓扑物态，使得Hubbard模型成为研究的热潮。文献[40]指出可以用光晶格中的冷原子系统模拟二维Hubbard—Hofstadter模型，对其中Hubbard项作类似于BCS平均场近似的处理。发现Majorana—Kramers对(MKPs)的存在，为高阶拓扑现象的发现提供支持。同样的文献[41]也证明在具有Rashba和Dresselhaus两种SOC的二维拓展Hubbard模型中得到了一阶(存在零能边缘态)和二阶(存在零能角态)拓扑超导相，改变温度或化学势可以有效地调控这种相变。
4.3非厄米
随着研究的逐步深入，对于量子系统的研究从封闭系统走向开放系统，即非厄米体系，展现出其独有的物理效应，如：奇异点(EPs)[42]、非厄米趋肤效应[43]等。而光学系统中的冷原子由于其可操纵性，可以很好地模拟非厄米效应[44—46]。
文献[47]中提出可以利用冷原子双光子拉曼辅助跃迁实现哈密顿量中的跃迁项，通过射频脉冲将原子共振转移得到损耗项，从而获得非厄米的三维连续哈密顿量：
计算非厄米的能谱，发现高维体系中也存在拓扑相，同时与厄米情形不同，非厄米系统的零能表面模受晶格尺寸调控：由于表面模局域化程度不高，当格点数N太小时，零能模并不严格局域在边缘，所以表面模消失，但是增强非厄米强度，即增大γz，可以使表面模重新出现。另外，改变动量可以调控表面模局域中心的位置。
非厄米体系的体边对应关系由非布洛赫理论描述，由此计算拓扑不变量的积分区间并非是厄米情形的布里渊区，而是由广义布里渊哈密顿量得到的广义布里渊区。
厄米情形的许多模型在引入非厄米的作用时会涌现出许多新奇的行为，如文献[48]中研究了非厄米费米超导的多体行为，其中非厄米项是二体非弹性碰撞产生的损耗γ。体系s波配对的非厄米的Bardeen—Cooper—Schriefferv(BCS)有效哈密顿量为
05
量子计算
量子计算近年来发展迅速。冷原子的高度可操控性和可拓展性，使得其有很大可能实现高比特的量子纠缠，文献[52]中首次提出可以用冷中性原子构建量子计算机。中国科学技术大学潘建伟团队自年开始，对超冷原子光晶格展开研究，实现了多对高保真度的超冷87Rb原子纠缠[53]，模拟了拓扑量子计算的任意子激发模型[54]。在年，他们首次提出新的原子冷却方案[55]：在晶格中调整局域化学式，使绝缘态冷原子样品与大态密度的超流相交替出现(图11(a))，两相中的原子通过隧穿交换熵，热量以低能激发的形式储存在超流相中。通过将超流相中的原子移除并调控晶格势，可以得到超过点位的低温均匀填充态。通过三能级的超交换作用，可以得到对两原子比特纠缠门。由图11(b)中自旋关联的测量值，可以得到实验生成密度矩阵的保真度为99.3±0.1%，远高于之前实验的中性原子二量子比特门的保真度。
06
结语
冷原子体系由于其宏观的量子特性和高度可调控性提供了一个全新的研究平台，其新颖量子态和奇异物性的研究是国际上具有前瞻性和挑战性的前沿领域。冷原子在量子精密测量、量子模拟和量子计算等领域有着重要的应用价值。超冷原子作为一种微观尺度上可调控的多体系统，在研究和调控宏观新奇物态上具有独特优势，并成为物理学的重要分支。近十年里，超冷原子呈现出诸多重要的新发展，为非常规物态的实现和研究提供了全新的视角和可靠的实验平台。其中，基于超冷原子的拓扑物态研究从无相互作用或弱关联的拓扑绝缘态开始延伸至相互作用强关联体系；基于少内禀自由度、短程作用超冷气体的研究推进到多自由度、长程相互作用的新型超冷气体平台，有效提升对新型多体物态的模拟能力；基于平衡态理论的物态表征和传统研究推进到常规凝聚态系统难以调控的远离平衡的多体动力学系统，催生表征动力学物态的新概念和基本理论，拓宽对非平衡物理的认知。这些超越常规的多体物态研究正在为凝聚态物理、统计物理、非平衡场论等多个物理学科打开一扇全新的通向新物理基础规律和新科学发现的大门，同时也将启迪全新的功能器件技术和高科技产业。
在以冷原子系综为基础的精密测量领域取得了丰硕的成果：冷原子喷泉频标、干涉仪测量引力常数、魔术波长光晶格频标、冷离子频标、冷原子重力仪、冷原子陀螺仪等新型冷原子精密测量仪器正在成为计量标准或是多次刷新了精密测量的记录，有些应用也已经在勘探、*事国防等方面发挥着重要作用。极低的热噪声、良好的相干性、灵活的操控性以及极高的信噪比使得冷原子系统可以用于制造量子信息的新型量子器件，也使其在量子精密测量、量子信息存储、量子信息传输等重要量子信息领域有着重大科学意义与技术应用。这些研究方向代表了未来信息技术发展的重要战略趋势，是世界各国展开激烈竞争的下一代量子信息体系的焦点，并极有可能对人类社会的经济发展产生巨大的影响。
参考文献
[1]HnschTW，SchawlowAL.OpticsCommunications，，13(1)：68
[2]ProdanJV，PhillipsWD，MetcalfH.Phys.Rev.Lett.，，49(16)：
[3]ChuS，HollbergL，BjorkholmJEetal.Phys.Rev.Lett.，，55(1)：48
[4]BohnJL，ReyAM，YeJ.Science，，()：2
[5]LanginTK，GormanGM，KillianTC.Science，，()：61
[6]BakerC，BertscheW，CapraAetal.Nature，，()：35
[7]GisbertAT，PiovellaN，BachelardR.PhysicalReviewA，，(1)：
[8]EndresM，BernienH，KeeslingAetal.Science，，()：4
[9]LeeJH，JungH，ChoiJYetal.PhysicalReviewA，，(6)：
[10]HaoL.TheEuropeanPhysicalJournalD，，75(5)：1
[11]PerryA，SugawaS，Salces-CarcobaFetal.OpticsExpress，，29(11)：
[12]GrossC，VogtT，LiW.Phys.Rev.Lett.，，(5)：
[13]HostenO，EngelsenNJ，KrishnakumarRetal.Nature，，()：
[14]SzigetiSS，NolanSP，CloseJDetal.Phys.Rev.Lett.，，(10)：402
[15]ChenG，AharonN，SunYNetal.NatureCommunications，，9(1)：1
[16]HouZ，TangJF，ChenHetal.ScienceAdvances，，7(1)：eabd
[17]LuXM，WangXG.Phys.Rev.Lett.，，(12)：7
[18]BlochI，DalibardJ，NascimbeneS.NaturePhysics，，8(4)：
[19]NohC，AngelakisDG.Rep.Prog.Phys.，，80(1)：37
[20]KengneE，LiuWM，MalomedBA.PhysicsReports，，：1
[21]ZhuHJ，ZhangGF，ZhuangLetal.Phys.Rev.Lett.，，(22)：
[22]StringariS.Phys.Rev.Lett.，，(14)：
[23]KroezeRM，GuoY，LevBL.Phys.Rev.Lett.，，(16)：
[24]SunW，WangBZ，XuXTetal.Phys.Rev.Lett.，，(15)：
[25]ZhangW.SyntheticSpin-OrbitCouplinginColdAtoms.WorldScientific，
[26]KolkowitzS，BromleySL，BothwellTetal.Nature，，()：66
[27]LaiW，MaYQ，ZhuangLetal.Phys.Rev.Lett.，，(22)：
[28]ChenL，PuH，ZhangY.PhysicalReviewA，，93(1)：
[29]ChenHR，LinKY，ChenPKetal.Phys.Rev.Lett.，，(11)：
[30]ZhangD，GaoT，ZouPetal.Phys.Rev.Lett.，，(11)：
[31]ChenKJ，WuF，PengSGetal.Phys.Rev.Lett.，，(26)：
[32]ChenL，ZhangY，PuH.Phys.Rev.Lett.，，(19)：
[33]ZhuHJ，XuK，ZhangGFetal.Phys.Rev.Lett.，，(5)：
[34]LeggettAJ.Phys.Rev.Lett.，，25(22)：
[35]HanW，ZhangXF，WangDSetal.Phys.Rev.Lett.，，(3)：7
[36]LiaoR.PhysicalReviewLetters，，(14)：
[37]HenkelN，NathR，PohlT.Phys.Rev.Lett.，，(19)：
[38]MazurenkoA，ChiuCS，JiGetal.Nature，，()：
[39]BrownPT，MitraD，Guardado-SanchezEetal.Science，，()：
[40]ZengC，StanescuT，ZhangCetal.Phys.Rev.Lett.，，(6)：
[41]KheirkhahM，YanZ，NagaiYetal.Phys.Rev.Lett.，，(1)：
[42]RotterI.JournalofPhysicsA：MathematicalandTheoretical，，42(15)：
[43]KawabataK，SatoM，ShiozakiK.PhysicalReviewB，，(20)：205
[44]LiJ，HarterAK，LiuJetal.NatureCommunications，，10(1)：1
[45]LiL，LeeCH，GongJ.Phys.Rev.Lett.，，(25)：
[46]GouW，ChenT，XieDetal.Phys.Rev.Lett.，，(7)：
[47]HeP，FuJH，ZhangDWetal.PhysicalReviewA，，(2)：
[48]YamamotoK，NakagawaM，AdachiKetal.Phys.Rev.Lett.，，(12)：601
[49]WuB，NiuQ.PhysicalReviewA，，64(6)：
[50]NguyenJH，LuoD，HuletRG.Science，，()：
[51]LyuC，LvC，ZhouQ.Phys.Rev.Lett.，，(25)：
[52]JakschD，CiracJ，ZollerPetal.Phys.Rev.Lett.，，85(10)：
[53]DaiHN，YangB，ReingruberAetal.NaturePhysics，，12(8)：
[54]DaiHN，YangB，ReingruberAetal.NaturePhysics，，13(12)：
[55]YangB,SunH,HuangCJetal.Science，，()：
本文经授权转载自