金刚石氮-空位色心的原子自旋声子耦合机理*

2021-03-26沈翔赵立业黄璞孔熙季鲁敏

物理学报 2021年6期

沈翔赵立业† 黄璞孔熙季鲁敏

1) (东南大学仪器科学与工程学院，微惯性仪表与先进导航技术教育部重点实验室，南京 210096)

2) (南京大学物理学院，南京 210093)

1 引言

氮-空位(Nitrogen-Vacancy， NV)色心是金刚石中的一种原子点缺陷.自从实现NV 色心的光探测磁共振以来， NV 色心受到了人们的广泛关注[1].室温下， NV 色心这一固态单自旋体系具有毫秒量级的相干时间，可用光学共聚焦系统进行初始化和读出，并且利用交变磁场实现单个自旋量子态的调控.这些优良性质使之被用于量子计算[2]、量子模拟[3]及量子物理基本问题[4-7]等研究.近年来，基于体系设计的纳米分辨率量子探测器，在物理学、化学、生物学等领域具有广泛而重要的应用前景[8].

将声子场耦合到量子位引起了量子声学的研究热潮.声子是晶体中离子或原子振动的量子化表现[9].声子特征与晶体材料的输运性质及多体相互作用下的元激发性质、元激发寿命、谱线宽度、磁化率、热性质等都密切相关[10，11].因声子波长小、无辐射等优点已经应用到量子计算[12-16]和超导量子位计算[17-20]中.通过将声子机械或电磁耦合到人造含掺物(如NV 色心)，研究者已开展自旋纠缠、自旋压缩和声子换能的研究[12，21-24].学者通过研究发现NV 色心的自旋态与金刚石的声子模式相关[25-40].NV 色心声子-自旋耦合强度的形变势与能隙相关[41]，其自旋相干性导致基态自旋的声子耦合非常弱[42，43]，激发态下NV 色心自旋比基态自旋的声子耦合程度大6 个量级[44，45].因此可通过声子辅助多物理场对NV 色心进行耦合控制[46].

目前，学者工作主要集中在研究声子作用下NV 色心的自旋跃迁性质，即“声子耦合作用于NV 色心，其自旋跃迁性质会发生什么变化? ”这一命题.目前对金刚石NV 色心的原子自旋声子耦合机理研究很少，特别是缺乏NV 色心自旋跃迁的声子耦合介导机理和声子模式下含NV 色心金刚石热力学性质的文献报道.本文重点研究金刚石氮-空位色心的原子自旋声子耦合机理，即“声子如何耦合作用于NV 色心? ”这一关键命题.本文引入了耦合声子场对NV 色心自旋进行共振调控，重点研究声子模式下NV 色心的结构特性、声子特性和热力学特性，为后续研究声子辅助操控NV 色心原子自旋，实现自旋跃迁效率的提高具有一定的指导意义.本文各部分研究内容阐述如下.

第一部分: 介绍了NV 色心量子传感的应用背景，具体阐述了目前声子耦合操控NV 色心自旋跃迁的研究现状.

第二部分: 阐述了含NV 色心金刚石的晶体结构、自旋跃迁性质和量子敏感机理.

第三部分: 首先，基于波函数和晶格的点阵位移矢量关系，分析了声子与晶格能量交互作用，研究了基于声子共振调控的NV 色心的自旋跃迁机理，建立了基于应变诱导的能量转移声子-自旋交互耦合激发模型.其次，基于NV 色心晶格振动理论，引入满足布洛赫定理的系数矩阵，建立了不同轴向NV 色心第一布里渊区特征区域的声子谱特征.同时，基于德拜(Debye)模型，考虑热膨胀效应，解析该声子共振系统的声子热平衡性质，并对其比热性质进行研究.

第四部分: 基于分子动力学仿真软件CASTEP 和密度泛函理论进行第一性原理研究，构建了声子模式下不同轴向NV 色心的结构优化模型，并分析了其结构特性、声子特性和热力学特性.

第五部分: 对本文工作进行总结.

2 NV 色心的性质及量子敏感原理

2.1 NV 色心的结构和自旋跃迁性质

金刚石NV 色心的几何结构如图1(a)所示，当金刚石晶格中有一个N 原子取代了C 原子，而相邻位置又缺少一个C 原子而留下空位时，就形成了NV 色心.大多数NV 色心带有一个负电荷或处于电中性状态，分别记为NV—和NV0.本文中提到的NV 色心均指NV—，它有6 个价电子，占据4 个分子轨道.NV 色心的能级结构及自旋跃迁性质如图1(b)所示.第一激发态与基态的能级差为1.945 eV，对应的波长是637 nm，采用532 nm 波长的激光可激发这一跃迁.处于激发态的NV 色心有两条主要的退激发途径，其中之一是发射荧光回到基态，另一条途径是发生系间窜越，经由单态1A1和1E 回到基态，这一过程不辐射荧光[8].

由于金刚石具有 C3V对称性，在金刚石样品中NV 色心存在4 个不同的量子化轴，分别为[1， 1， 1]轴向、[1， —1， —1]轴向、[—1， 1， —1]轴向和[—1， —1， 1]轴向.NV 色心量子化轴如图2 所示.金刚石中4 个不同轴向NV 色心分布及NV 坐标系如图3 所示.

图1 NV 色心几何结构和自旋跃迁性质 (a)几何结构; (b)能级结构及自旋跃迁性质Fig.1.Structures and spin transition properties of a negatively charged NV center: (a) Geometric structure; (b) energy level structure and spin transition properties.

图2 NV 色心量子化轴示意图Fig.2.Schematic diagram of quantization axis for NV center.

图3 金刚石中4 个轴向NV 色心分布及NV 坐标系Fig.3.Four axial NV center distributions and their NV coordinate systems in diamond.

2.2 NV 色心的量子敏感原理

NV 色心基态自旋能量的哈密顿方程为

式中D 为NV 色心的零场分裂参量， D = 2.87 GHz;B 为外部磁场矢量; g 为波谱分裂因子， g =2 ;μB为玻尔磁子， μB=1.4 MHz/G; ℏ =h/(2π) 是约化普朗克常数; S 为不确定度，其表达式为

其中， SX， SY和SZ分别为量子化轴X， Y 和Z 方向的不确定度.

(1)式等号右侧第一项为NV 色心零场分裂哈密顿量; 第二项为加外部磁场后NV 色心自旋作用哈密顿量.

因此，由两个退激发频率(通过荧光检测)就可获得磁场强度B 和相对于该NV 轴的极化角 θ :

3 NV 色心的声子共振调控机理及热力学解析

3.1 NV 色心自旋跃迁的声子共振调控机理

声子本质是声子场在金刚石内部传输时，介质的晶格振动能量的量子.在基于声子共振调控的NV 色心原子自旋量子态操控系统中，声子与自旋的交互作用实为晶格振动与自旋的相互作用.自旋在晶格的周期性势场中，由于晶格原子不停地在其平衡位置附近发生小振动，从而破坏了晶格的周期性势场，产生了附加形变势.自旋因受到附加势场的作用，激发势能发生改变[47].

绝热近似下，金刚石晶格的运动简化为两个相对较小的子系统(电子系统和原子实系统)的运动.对给定的电子系统状态n，原子实系统受到有效势场 VN(R) 为原子实间的库仑相互作用 VLL(R) 与电子能 En(R) 之和，即:

图4 (a)声子场共振结构示意图; (b)声子场共振调控机理示意图[48]Fig.4.(a) Schematic diagram of phonon field resonance structure; (b) mechanism diagram of phonon field resonance control.

则原子实系统运动的哈密顿方程表达如下:

其中， X (R) 是原子实系统运动轨迹向量.

采用周期性边界条件，在金刚石第n 个原胞中，第l 个原子的平衡位置为 Rnα=Rn+Rl，Rn和 Rl分别为原胞位置和原子l 在 l′原胞中的相对位置，则金刚石面心立方晶体的原子振动方程转化为如下形式:

由于金刚石是对周期排布的原子晶体，在正则坐标系下对(9)式进行简谐近似，转化为如下形式的格波

由上述分析可知，在声子场作用下，可实现NV色心晶格原子的定向振动，声子场共振结构及调控机理示意图分别如图4(a)和图4(b)[48]所示.

在无声子场作用时，由拉比频率 Ω1光场驱动的NV 色心自旋跃迁系统在大偶极子 Δ 失谐极限内，等效于的两个较低状态的光驱动之间的自旋跃迁.在含声子场作用时，由声子场驱动系统从 |g〉到|e〉的声子共振调控自旋跃迁，其有效拉莫尔频率为ΩP.在大偶极子 Δ 失谐极限下，该系统等效于2 个较低状态的声子阶梯之间的自旋跃迁，其中n 表示声子数.

NV 色心的激发态与波长晶格应变强烈耦合.这种电子-声子耦合可以导致应变诱导的能量转移，也可以导致相关激发态的状态混合.对于 Ey态，描述应变诱导的能量转移的电子-声子相互作用哈密顿量模型为

其中， b 是声子的湮灭算符; b+是声子的产生算符;是电子-声子有效耦合率;εc表示晶格振动时的带边能量; δ εc是形变势，且是晶格振动的纵向声学模引起的晶格体积V 的局域变化; km是声子模的波数;m 是有效的力学质量振荡器; ωm是声子振荡频率.

对于从 ms=0 到 Ey态的声子共振调控的自旋跃迁，与位移成线性关系的有效哈密顿量 HR为[29]

其中， Ω0是耦合到 ms=0 到 Ey态跃迁的激光场的拉莫尔频率.因此，声子共振调控的自旋跃迁(红色带边跃迁)的有效拉比频率 ΩP由下式给出:

其中，〈 n〉是平均声子数.当激光从红色带边过渡到蓝色带边时，其哈密顿量类似.

将声子共振调控的自旋跃迁结合到三能级系统中，并使用激发态电子-声子耦合来介导NV 色心自旋和机械自由度之间的相互作用.电子被困在两个较低的态中，形成由两个较低能级的相干叠加导致的暗态，该暗态可以介导和控制NV 色心自旋状态与相关声子模式之间的相互作用，由下式给出:

其中， ΩP和 Ω±是分别耦合到 ms=0 和 ms=±1 状态的两个跃迁的拉比频率.结合(1)式，对于声子共振调控下NV 色心的量子测量模型，其哈密顿方程为

从(15)式可以看出，以声子场辅助射频场和磁场共振耦合调控，可实现NV 色心原子自旋的共振耦合调控.

3.2 NV 色心的声子谱解析

根据(11)式和(12)式，声子场辅助射频场和磁场共振耦合调控取决于电子-声子相互作用He-phonon和声子共振调控的自旋跃迁作用 HR.He-phonon和 HR与声子模式下晶格的形变势 δ εc密切相关，关系方程如下:

其中， fNV为NV 色心晶格结构特征函数; ωσ(q) 为频率为 ωσ，波矢为 q 的声子谱; VLA为晶格的纵向声学模; qz为波矢为 q 的z 向分量.因此对NV 色心声子谱 ωσ(q) 的解析至关重要.

由于(10)式所示的格波色散关系在波矢空间是周期为 π /a 的周期性函数( E (k)=E(k+π/a) )，因此只需用第一布里渊区中的波矢来描述能带电子、点阵振动和自旋波的状态.纯净无缺陷的金刚石为面心立方晶体，其第一布里渊区模型如图5(a)所示.该模型为十四面体结构，包含 Γ 点、L 点、U点、X 点、W 点和K 点6 个布里渊区特征点，还包括 Λ 线、 Σ 线和 Δ 线3 个布里渊区特征线.对于含NV 色心的金刚石，根据晶格动力学理论[49]，本文通过在单个NV 色心晶格的交界面上对称扩展方式建立第一布里渊区模型，如图5(b)所示，该模型为六面体结构，包含 Γ 点、Z 点、Q 点、K 点、B 点和F 点6 个布里渊区特征点，还包括 Λ 线、 Σ 线、Δ线和M 线4 个布里渊区特征线.

图5 金刚石第一布里渊区特征 (a)不含NV 色心; (b)含NV 色心Fig.5.Characteristics of first Brillouin zone of diamond:(a) Without NV center; (b) contain NV center.

本文基于晶体原子振动方程((9)式)对金刚石NV 色心的声子谱进行解析，由于(9)式是3Nm(Nm为第一布里渊区内的原子数)个二阶微分方程联立耦合的方程组，为了对其进行声子谱解析，引入满足布洛赫(Bloch)定理的系数矩阵[49]，结合(10)式，得到晶格动力学矩阵元 Dαβ，表达如下:

晶格动力学矩阵元 Dαβ的非零解的条件是

其中， ωσ(q) 为频率为 ωσ，波矢 q =(qx，qy，qz) 的声子谱; δαβ表示原子沿 β 方向位移单位距离时原子作用力沿 α 方向分量的声子谱系数.

将(18)式展开，忽略三次及以上近邻力高阶项，可得到3 个声子谱 ωσ(q)(σ=1，2，3) 及对应的3 个极化向量 eqσ=(eqxσ， eqyσ， eqzσ) (表示金刚石原子的振动方向)，其动力学矩阵元分解式为

表1 不同轴向NV 色心的晶格动力学矩阵元的不对称关系Table 1.Asymmetrical relations of lattice dynamics matrix elements for NV centers of different axes.

其中， f1是最近邻力常数， f2是次近邻力常数，η =f2/f1是次近邻力常数与最近邻力常数之比，α是金刚石近邻原子间平衡间距， a =1.544 Å，波矢q =(qx，qy，qz).

对于四种轴向的NV 色心金刚石晶体，晶格动力学矩阵元的不对称关系如表1 所列.其中， k[1，1，1]，k[1，-1，-1]， k[-1，1，-1]和 k[-1，-1，1]分别是[1， 1， 1]轴向、[1， —1， —1]轴向、[—1， 1， —1]轴向和[—1， —1， 1]轴向NV 色心金刚石晶体的氮-空位不对称系数.

由于四种轴向NV 色心的金刚石晶体具有相似性，在进行声子谱解析时仅晶格动力学矩阵元不对称关系不同，所以本节仅以[1， 1， 1]轴向NV 色心的金刚石晶体进行详细解析.在[1， 1， 1]轴向NV 色心金刚石晶体中，令:

金刚石NV 色心的第一布里渊区一般点(非简并点)上的声子能量为

式中， X1σ和 X2σ分别为原子实系统的一阶本征态振动量和二阶本征态振动量，表达式如下:

其中， λσ为第 σ 本征态( σ =1，2，3 )的本征值.

由(19)—(23)式可以解析出[1， 1， 1]轴向NV 色心金刚石晶体特征线及特征面的声子谱，其中， Λ 特征线及 qz=0 特征面的解析结果如下.

1)第一布里渊区 Λ 特征线上( qx=qy=qz=q )的声子谱:

当 qx=qy=qz=q 时，令

(24)式化简为

根据(22)式， ω1， ω2， ω3分别对应的极化向量为

因此， ω1是纵向声子， ω2与 ω3是横向声子，而且横向声子在第一布里渊区 Λ 特征线上有二重简并.

2)第一布里渊区 qz=0 平面的声子谱.

第一布里渊区 qz=0 特征面的声子谱同(21)式，但其条件满足 qz=0 ，令:

可得

根据(22)式， ω1， ω2， ω3分别对应的极化向量为

其中，

因此， ω2是纵向声子， ω1与 ω3是横向声子，而且横向声子在第一布里渊区 qz=0 平面上有二重简并.

[1， 1， 1]轴向NV 色心金刚石第一布里渊区内的其余特征线或特征面的声子谱解析过程与 Λ 特征线和 qz=0 特征面类似，均基于(19)—(23)式仅变换波矢q 的声子谱条件，详细解析过程不再赘述.相关特征线的声子谱分析结果如表2 所列.

表2 [1， 1， 1]轴向NV 色心金刚石布里渊区特征线的声子谱解析结果Table 2.Phonon spectrum analysis results at the characteristic line of the Brillouin zone in the diamond with the NV center of [1， 1， 1] axis.

3.3 NV 色心的声子热力学解析

3.3.1 声子热平衡解析

在NV 色心自旋跃迁过程是在振动态声子热平衡条件下进行的，且达到热平衡的过程比自旋跃迁的速率快得多，所以需要对NV 色心声子共振系统的热力学性质进行解析.

在热平衡条件下，一个频率为 ωσ的振动模 ωσ(q) ，处于本征态 σ ，模中有n 个声子的概率 Pn，正比于玻尔兹曼(Boltzmann)因子: e xp[-nℏωσ(q)/(kBT)] ，kB是Boltzmann 常数， T 为声子热平衡温度.

令γ =e-nℏωσ(q)/(kBT)， (31)式变换为

频率为 ωσ(q) 的振动模中的热平均声子数〈n〉可以表示为

化简(33)式得到

代入 γ ，可得到声子热平衡温度函数，即:

为了定量计算不同轴向结构参数金刚石NV色心的声子热力学性质，结合(21)式和(35)式NV色心金刚石晶体的不同轴向结构参数声子热平衡温度为

其中，〈 n〉表示平均声子数.(36)式中，参数u， b，n 均与NV 色心金刚石晶体轴向结构有关.特别地，取 σ =1 时的声子热平衡温度，所以(36)式可转化为

对于[1， 1， 1]轴向NV 色心金刚石，声子热平衡温度解析结果如表3 所列，其中， TΛ， TΔ， TΣ和TM分别为声子极化方向为 Λ 线、 Δ 线、 Σ 线和M 线方向的声子热平衡温度.

3.3.2 比热解析

比热是反应晶体的重要物理量之一，本节基于Debye 模型对声子模式下NV 色心金刚石晶格三维模型进行比热解析，探究声子作用对NV 色心金刚石晶格的比热影响机理.

在金刚石NV 色心三维波矢空间内，格波的等频面为球面，在 q →q+dq 方向上波速为 vq的格波数目 d zq为

其中， V 是含NV 色心金刚石的晶格体积，且3/(vq)3=1/(vH)3+2/(vT)3; vH是横向声子波速;vT是纵向声子波速.

由(38)式可得到声子模式下的NV 色心金刚石的晶格振动密度 ρ [ωσ(q)] :

根据(39)式可得到声子模式下的NV 色心金刚石的晶格自由能:

根据(40)式可以得到声子模式下的NV 色心金刚石的晶格热容量:

为了定量计算不同轴向结构参数金刚石NV色心的声子热力学性质，根据(21)式，令ωσ(q)=参数u， b， n 均与NV 色心金刚石晶体轴向结构有关，将(21)式代入(40)式和(41)式可分别解析出NV 色心金刚石晶体的不同轴向结构晶格自由能和晶格热容量.

表3 [1， 1， 1]轴向NV 色心金刚石的声子热平衡温度解析结果Table 3.Phonon thermal equilibrium temperature analysis results of the diamond with the NV center of [1， 1， 1] axis.

4 NV 色心声子谱及热力学性质第一性原理计算

4.1 计算方法

本文基于应变诱导能量转移的电子-声子相互作用哈密顿量模型和密度泛函理论[50]，采用分子动力学仿真软件CASTEP 构建不同轴向的NV 色心金刚石晶体的结构模型，在声子共振调控条件下的氮-空位色心的结构优化模型，分析该优化模型不同NV 色心占位时的结构特性、声子特性和热力学特性.

本文的第一性原理计算是基于密度泛函理论，采用分子动力学软件CASTEP 通过平面波赝势方法实现.为了准确描述氮-空位原子域电子间的关联作用，计算基于广义密度近似的Ceperley Alder-Perdew Zunger 模型，晶胞中价电子与离子核的相互作用通过超软赝势Ultrasoft 描述.第一布里渊区内采用8 × 8 × 8 的Monkhorst-Pack 网格.自洽计算的收敛精度为总能量的变化低于10-8eV.晶胞内各原子的受力不超过0.01 eV/Å.

声子谱的计算是基于密度泛函微扰理论的局域密度近似下的Perdew-Burke-Ernzerh 模型实现.计算过程中，布里渊区内格波的波矢q 采用6 ×6 × 6 的无偏移网格.考虑材料泊松系数随体积的变化，声子热力学性质的计算采用晶格振动谐振近似下的Debye 模型.在采用有限粒子数来模拟实际体系中粒子的运动时，通过施加三维周期性边界条件使处于模拟体系中的粒子的运动空间成为赝无限来减少界面对体系中粒子振动的影响.在模拟体系设定的温度、压力和粒子的初始速度下，基于等温等压系综，根据埃瓦尔德(Ewald)求和技术[51]进行解析，每个平衡态计算10000 步，步长为1 fs.

4.2 计算结果

4.2.1 结构特性

基于分子动力学软件CASTEP 对不同轴向NV色心的晶格结构进行优化，使晶格能最低.不同轴向NV 色心金刚石的晶格能优化特征如图6 所示.从图6 可以看出，经过结构优化迭代后，晶格能降到平衡状态.结构优化后纯净无缺陷的金刚石晶格总能量为—309.802 eV， [1， 1， 1]轴向的NV 色心金刚石晶格能为—1198.548 eV，其余三种轴向的NV 色心金刚石晶格能相差不大，约为—1197.3 eV.所以，含NV 色心缺陷的金刚石晶格能小于纯净无缺陷的金刚石.

图6 不同轴向NV 色心金刚石的晶格能优化特征Fig.6.Lattice energy optimization characteristics for the diamond with NV centers of different axes.

晶格能优化后的不同轴向NV 色心金刚石的带隙特征如图7 所示.从图7 可以看出，无NV 色心金刚石的带隙特征是沿着第一布里渊区特征点W 点、L 点、 Γ 点、X 点、W 点和K 点连线展开的.四种轴向NV 色心金刚石的带隙特征是沿着第一布里渊区特征点 Γ 点、F 点、Q 点、Z 点和 Γ 点连线展开的.无NV 色心金刚石的带隙为4.325 eV.[1， 1， 1]轴向、[1， —1， —1]轴向、[—1， 1， —1]轴向和[—1， —1， 1]轴向的NV 色心金刚石的带隙分别为0.6257， 1.6012， 1.2295， 1.2864 eV.因此，相比于纯净无缺陷金刚石，含有NV 色心的金刚石减小了导带的最低点和价带的最高点的能量之差.电子由价带被激发到导带更容易，本征载流子浓度更高，电导率更高.对比四种不同轴向NV 色心金刚石的带隙特征可知， [1， 1， 1]轴向NV 色心金刚石的带隙最低，具有最高的本征载流子浓度和电导率.反之，[1， —1， —1]轴向NV 色心金刚石的带隙最高，其本征载流子浓度和电导率最低.

图7 不同轴向NV 色心金刚石的带隙特征 (a)无NV 色心; (b) [1， 1， 1]轴向; (c) [1， —1， —1]轴向; (d) [—1， 1， —1]轴向;(e) [—1， —1， 1]轴向Fig.7.Band gap characteristics for the diamond with NV centers of different axes: (a) Without NV center; (b) axis direction of[1， 1， 1]; (c) axis direction of [1， —1， —1]; (d) axis direction of [—1， 1， —1]; (e) axis direction of [—1， —1， 1].

不同轴向NV 色心金刚石的态密度曲线如图8 所示，可以看出，含NV 色心的金刚石态密度(density of states， DOS)分布较纯净无缺陷金刚石的局域尖峰更加明显.[1， 1， 1]轴向、[1， —1， —1]轴向、[—1， 1， —1]轴向和[—1， —1， 1]轴向的NV色心金刚石的最大DOS 分别为3.3149， 3.7512，3.5727， 3.5695 electrons/eV，对应的晶格能位分别为—5.78， —5.96， —5.83， —5.82 eV.纯净无缺陷的金刚石的最大DOS 仅为1.2516 electrons/eV，对应的晶格能位为8.29 eV.因此，含有NV 色心的金刚石相比于纯净无缺陷金刚石具有更强的电子局域化性质.从图8 还能看出不同轴向NV 色心金刚石的能带特性.由于赝能隙直接反映了该体系成键的共价性的强弱.赝能隙越宽，共价性越强.含NV 色心的金刚石的赝能隙较纯净无缺陷金刚石的赝能隙更窄.因此，其共价键较纯净无缺陷金刚石更弱.从DOS 曲线还可以看出不同轴向NV 色心金刚石在一定的晶格能区间均具有零带隙，与图7 结果相吻合.

4.2.2 声子特性

图9 为不同轴向NV 色心金刚石沿着第一布里渊区特征线的声子谱.从图9 可以看出，纯净无缺陷金刚石的原胞具有2 个原子，产生6 个色散关系(3 个声分支和3 个光学分支).含NV 色心金刚石的第一布里渊区模型是通过对称扩展交界面的方式形成的，其原胞具有7 个原子，产生21 个色散关系(6 个声分支和15 个光学分支).纯净无缺陷金刚石和四种不同轴向NV 色心的金刚石的禁带宽度最窄的地方均出现在 Γ 点处，且价带的最高点与导带最低点在同一点，故属于直接带隙.四种不同轴向NV 色心的金刚石声子谱略有差异，这是不同轴向NV 色心金刚石晶体的晶格动力学矩阵元不对称关系导致的.

图8 不同轴向NV 色心金刚石的态密度曲线Fig.8.State density curves of the diamond with NV centers of different axes.

图9 不同轴向NV 色心金刚石的声子谱 (a)无NV 色心; (b) [1， 1， 1]轴向; (c) [1， —1， —1]轴向; (d) [—1， 1， —1]轴向; (e) [—1，—1， 1]轴向Fig.9.Phonon spectrum curves of the diamond with NV centers of different axes: (a) Without NV center; (b) axis direction of[1， 1， 1]; (c) axis direction of [1， —1， —1]; (d) axis direction of [—1， 1， —1]; (e) axis direction of [—1， —1， 1].

分析不同轴向NV 色心金刚石声子态密度(图10)可发现，纯净无缺陷金刚石晶体中声子振动模式主要出现在10—38 THz 频段，特别在35—37 THz 频段最大声子态密度为0.13395 (unit cell)/THz，在此频段下出现声子振动模式的概率较高.四种不同轴向的声子态密度曲线略有差异，[1， 1， 1]轴向NV 色心金刚石具有较多的大概率声子振动模式频段，其最大声子态密度为0.10904(unit cell)/THz.[1， —1， —1]轴向、[—1， 1， —1]轴向和[—1， —1， 1]轴向的NV 色心金刚石的大概率声子振动模式频段基本一致.[1， —1， —1]轴向、[—1， 1，—1]轴向和[—1， —1， 1]轴向的NV 色心金刚石的最大声子态密度分别为0.12681， 0.1308， 0.13211(unit cell)/THz.由图10 可知，四种不同轴向NV色心金刚石的声子主共振频段均处于THz 量级，声子耦合激励至主共振频率技术较难实现.NV 色心金刚石的次共振频率约为[800， 1200] MHz.根据该频段设计等效的声表面波共振机构，其叉指的宽度w 为1.5 μm.由于声表面波在金刚石中传播速度vs约为5600 m/s，且声表面波中心频率ωm=vs/4w.推算设计的声表面波叉指电极中心频率约为930 MHz，处于次共振频率范围.在该声子共振调控参数条件下，声子共振调控方法可有效增大氮-空位色心的自旋跃迁概率.

图10 不同轴向NV 色心金刚石的声子态密度曲线Fig.10.Phonon state density curves of the diamond with NV centers of different axes.

4.2.3 热力学特性

Debye 温度不仅反映晶体点阵的动畸变程度，还是物质原子间结合力的表征，物质的弹性、硬度、熔点和比热等物理量均与它有关.根据不同轴向NV 色心金刚石的Debye 温度曲线(图11)可知，当温度为25 K 时， Debye 温度均达到C 点，纯净无缺陷金刚石、[1， 1， 1]轴向、[1， —1， —1]轴向、[—1，1， —1]轴向和[—1， —1， 1]轴向的NV 色心金刚石的C 点Debye 温度值分别为2614.50934， 3379.62378，1233.70678， 1145.47014， 1165.59051 K.当温度为500 K 以上时，不同轴向NV 色心金刚石的Debye温度均达到稳定，纯净无缺陷金刚石、[1， 1， 1]轴向、[1， —1， —1]轴向、[—1， 1， —1]轴向和[—1， —1， 1]轴向的NV 色心金刚石的稳定Debye温度值分别为1489.99201， 1575.97472， 1542.95619， 1547.77852，1547.94141 K.从稳定Debye 温度可知，纯净无缺陷金刚石的稳定Debye 温度高于含NV 色心的金刚石，因此，纯净无缺陷金刚石的原子健结合力更强，熔点更高.同理，对于四种不同轴向NV 色心的金刚石来说， [1， 1， 1]轴向NV 色心金刚石的原子健结合力和熔点略大于其余三种轴向.

图11 不同轴向NV 色心金刚石的Debye 温度特征 (a)特征曲线; (b)特征值Fig.11.Debye temperture characteristics of the diamond with NV centers of different axes: (a) Characteristic curves;(b) characteristic values.

不同轴向NV 色心金刚石的声子热力学曲线(图12)表明，热焓、热熵随着温度的升高而增加，且纯净无缺陷金刚石的热焓随着温度升高，其增长幅度远小于含NV 色心的金刚石(图12(a)).同时，热熵的增长幅度也远小于含NV 色心的金刚石(图12(b)).另外，自由能随着温度的升高而下降，且纯净无缺陷金刚石自由能的下降幅度不及含NV色心的金刚石(图12(c)).四种不同轴向NV 色心的金刚石的声子热力学性质差异不明显.纯净无缺陷金刚石与含NV 色心的金刚石相比，具有更加稳定的热力学性质.图12(a)、图12(b)和图12(c)还分别基于三次多项式方程( y =A+Bx2+Cx3)对不同轴向NV 色心金刚石的热焓、热熵和自由能曲线进行拟合，其拟合方程分别列入图12 中.基于上述拟合方程可描述和预测不同轴向NV 色心金刚石的热焓、热熵和自由能等声子热力学性质.研究结果还表明，系统声子模式的演化依赖于NV 色心的占位，声子模式强化伴随着热力学熵的降低.

不同轴向NV 色心金刚石的热容曲线(图13)表明，热容随着温度的升高而增加，且逐渐趋向于平衡.高温下纯净无缺陷金刚石的热容远小于含NV色心的金刚石.四种不同轴向NV 色心的金刚石的热容性质差异不明显.特别的，当温度为250 K 时，纯净无缺陷金刚石、[1， 1， 1]轴向、[1， —1， —1]轴向、[—1， 1， —1]轴向和[—1， —1， 1]轴向的热容分别为1.999， 11.874， 12.751， 12.473， 12.469 cal/(cell·K).当温度为1000 K 时，对应的热容分别为10.108，36.975， 37.157， 37.131， 37.130 cal/(cell·K).

图13 不同轴向NV 色心金刚石的热容特性 (a)热容曲线; (b)热容值Fig.13.Heat capacity characteristics of the diamond with NV centers of different axes: (a) Heat capacity curves; (b)heat capacity values.

5 结论

本文引入耦合声子场对NV 色心自旋进行共振调控，重点研究声子模式下NV 色心的自旋操控机理和热力学性质，研究内容及相关结论如下:

1)基于波函数和晶格的点阵位移矢量关系，分析了声子与晶格能量交互作用，研究了基于声子共振调控的NV 色心的自旋跃迁机理，建立了基于应变诱导的能量转移声子-自旋交互耦合激发模型.

2)基于NV 色心晶格振动理论，引入满足布洛赫定理的系数矩阵，建立了[1， 1， 1]轴向NV 色心第一布里渊区特征区域的声子谱特征函数，并对相应的极化向量进行了解析.解析结果表明: 横向声子 ω2与 ω3在NV 色心金刚石晶格的第一布里渊区 Λ 特征线产生二重简并; 横向声子 ω1与 ω3在其第一布里渊区 qz=0 平面上产生二重简并.

3)基于Debye 模型，考虑热膨胀效应，对声子共振系统的声子热平衡性质进行解析，建立[1， 1， 1]轴向NV 色心沿着声子极化方向的声子热平衡温度模型.同时，对声子模式下不同轴向NV 色心结构的晶格自由能和晶格热容量进行了解析.解析结果表明: NV 色心金刚石晶体的声子热平衡温度，晶格自由能和晶格热容量均与其晶体结构参数密切相关.

4)基于分子动力学仿真软件CASTEP 和密度泛函理论进行第一性原理研究，构建了声子模式下不同轴向NV 色心的结构优化模型，并分析了其结构特性、声子特性和热力学特性.研究结果表明，系统声子模式的演化依赖于氮-空位的占位，声子模式强化伴随着热力学熵的降低.含NV 色心的金刚石共价键较纯净无缺陷金刚石更弱，热力学性质更不稳定.含NV 色心金刚石的声子主共振频段均处于THz 量级，次共振频率约为[800， 1200] MHz.声子耦合激励至主共振频率技术较难实现.根据次共振频段设计叉指宽度为1.5 μm 的声表面波共振机构，其中心频率约为930 MHz.在该声子共振调控参数条件下，声子共振调控方法可有效增大氮-空位色心的自旋跃迁概率，实现氮-空位色心原子自旋操控效率的提高.

本文工作为后续研究声子辅助操控NV 色心原子自旋，实现自旋跃迁效率的提高具有一定的指导意义.