石墨变金刚石：蓝丝黛尔石与钻石的动力学竞争

杨金龙

(中国科学技术大学化学与材料科学学院，合肥230026)

石墨变金刚石：蓝丝黛尔石与钻石的动力学竞争

杨金龙

(中国科学技术大学化学与材料科学学院，合肥230026)

六方金刚石(hexagonal diamond，HD)，又称蓝丝黛尔石(Lonsdaleite)是自然界中神秘的物质。六方金刚石在自然界非常稀少，科学界普遍认为自然界中的六方金刚石是流星上的石墨坠入地球时形成的。流星撞击地球会产生了巨大的热量与压强，从而导致石墨发生相变形成金刚石。第一次发现六方金刚石是在美国亚利桑那州的巴林杰陨石坑(Barringermeteror crater)中的魔谷陨石(canyon diablometeorite)和印度阿萨姆邦的阿帕拉陨石(Goalparameteorite)中1。六方金刚石在地质研究中有着重要的作用。例如，在距今1.3万年的冰川世纪(Ice Age)，又称为更新世(Pleistocene)，地球上发生大量的大型动物灭绝。许多学者认为，这与流星撞击地球有关，而目前科学家们判定是否有流星撞击事件的重要证据就是在相应的地质层中，是否可以找到六方金刚石2。相比于六方金刚石，立方金刚石(cubic diamond，CD)则经常出现在日常生活，通常被称为钻石。人类早就知道石墨和钻石是同素异构体，利用高温高压固态相变虽然可以将石墨转变成金刚石，但是一般需要巨大的压强和非常高的温度3,4(见图1a，相图)；即使如此，得到的钻石颗粒度也极小。

图1 (a)温度－压力(p－T)碳相图(G：石墨；D：金刚石)4。(b)通过采用SSW方法在15GPa下进行势能面搜索得到的5种亚稳的双相异质节，图中虚线为金刚石密排面，即HD(001)，CD(111);gh2中界面与密排面垂直；橙色标记界面原子。(c)三条最低反应路径的能量势能面图：它们分别经过gh1，gh2和gc1等界面结构，其中TS表示过渡态，MS表示亚稳态10

寻找石墨合成金刚石的最优合成条件是科学家过去60多年来不断努力的研究方向。早在1967年，在对石墨逐步增加压强时发现，能合成稳定存在的金刚石的起始压强大约为15 GPa，起始温度大约1000°C。然而，经过原子结构分析发现，在这种条件下合成的金刚石竟然不是常见的立方金刚石，而是与来自宇宙的六方金刚石有着同样的结构1,3。因此，理解石墨到六方金刚石的反应机理，特别是从原子层次揭示其机理，不仅对合成钻石有着基础性的指导意义，而且对深入了解陨石中钻石的具体来源也有着重要意义。更为广泛的是，在高压物理中，在不同的压强中物质将表现成完全不一样的性质，深入理解石墨到金刚石的相变机理有普适的基础物理意义。尽管过去几十年来，人们采用了许多现代实验和理论方法对石墨到金刚石的相变进行了广泛研究，但是由于难于了解固体相变的原子级别反应路径，对于为何六方金刚石比起立方金刚石更加容易在固固相变中被合成的这一关键问题，一直没有答案。

2013年复旦大学刘智攀教授课题组开发了新一代势能面全局搜索方法，即随机势能面行走(stochastic surfacewalking，SSW)方法5,6。SSW方法与常见结构搜索算法的区别在于，SSW方法能在高维度势能面上给定初始结构的周围系统地搜寻亚稳态，并能通过高通量计算搜索相变发生的最优反应路径。该方法已成功应用于TiO2、ZrO2和MnO2等体系的固态相变研究7－9。最近，他们采用SSW方法系统研究了石墨到金刚石的固体相变过程，通过寻找所有可能的异相结中间态，为石墨到金刚石的成核相变路径的研究打开了大门。研究确定了七种异相结的原子结构(见图1b)，并发现石墨到六方金刚石的相变过程中会形成三维共格成核核心，而到立方金刚石的相变过程中必然伴随非共格界面的存在，这就导致了石墨到六方金刚石的相变动力学更为容易(见图1c)。同时，该研究还指出，相变过程中的应变力大小对相变势垒起到至关重要的影响。该研究工作最近发表在Journal of the American Chem ical Society杂志上10。这一研究工作预测了石墨到金刚石固固相变的成核动力学，为人工合成金刚石提供了新的思路。

(1)Hanneman,R.E.;Strong,H.M.;Bundy,F.P.Science 1967,

(2)Kennett,D.J.;West,A.;West,G.J.;Bunch,T.E.;Culleton,B. J.;Erlandson,J.M.;Hee,S.S.Q.;Johnson,J.R.;Mercer,C.; Shen,F.;Sellers,M.;Staffor,T.W.;Stich,A.;Weaver,J.C.; Wittke,J.H.;Wolbach,W.S.Proc.Natl.Acad.Sci.U.S.A. 2009,106,12623.doi:10.1073/pnas.0906374106

(3)Bundy,F.P.;Kasper,J.S.J.Chem.Phys.1967,46,3437.doi: 10.1063/1.1841236

(4)Bundy,F.P.Phys.A 1989,156,169.doi:10.1016/0378-4371(89) 90115-5

(5)Shang,C.;Liu,Z.P.J.Chem.Theory Comput.2013,9,1838. doi:10.1021/ct301010b

(6)Zhang,X.J.;Liu,Z.P.J.Chem.Theory Comput.2015,11, 4885.doi：10.1021/acs.jctc.5b00641

(7)Guan,S.H.;Zhang,X.J.;Liu,Z.P.J.Am.Chem.Soc.2015, 137,8010.doi:10.1021/jacs.5b04528

(8)Zhu,S.C.;Xie,S.H.;Liu,Z.P.J.Am.Chem.Soc.2015,137, 11532.doi:10.1021/jacs.5b07734

(9)Li,Y.F.;Zhu,S.C.;Liu,Z.P.J.Am.Chem.Soc.2016,138, 5371.doi:10.1021/jacs.6b01768

(10)Xie,Y.P.;Zhang,X.J.;Liu,Z.P.J.Am.Chem.Soc.2017, doi:10.1021/jacs.6b11193

Graphite to Diam ond:Kinetics Selec tivity Lead ing to Hexagonal Diam ond o r Cubic Diam ond

YANG Jin-Long
(SchoolofChemistry and Materials Science,University ofScience and Technology ofChina,Hefei230026,P.R.China)

995.

10.1126/science.155.3765.995

doi:10.3866/PKU.WHXB201702172

www.whxb.pku.edu.cn