First-principles calculations of diamond/copper (silver, titanium carbide) interface properties
-
摘要: 采用第一性原理计算的方法研究金刚石/铜、金刚石/银、金刚石/碳化钛3种界面的结构、电子结构和传热。结果表明：金刚石/碳化钛的界面结构最为稳定，界面间距（1.990 Å）最小，界面黏附功（5.578 J/m2）最大，结合强度最高。电子态密度、马利肯布居分析、差分电荷密度、径向分布函数的结果表明金刚石/碳化钛存在较多的电荷转移和较强键合作用。 声子态密度的计算结果表明金刚石/碳化钛的界面热阻较低。
-
关键词:
- 金属基金刚石复合材料 /
- 第一性原理计算 /
- 界面结构 /
- 界面电子结构 /
- 界面传热
Abstract: The structure, electrical structure, and heat transmission of diamond/copper, diamond/silver, and diamond/titanium carbide surfaces have been investigated using first-principles calculations. The results show that the diamond/titanium carbide interfacial structure is the most stable, with the shortest interfacial distance (1.990 Å), the greatest interfacial adhesion effort (5.578 J/m2), and the best bond strength. The results of the electronic density of states, mulliken population analysis ,charge density difference, and radial distribution function indicate the presence of more charge transfer and stronger bonding in diamond/titanium carbide. According to the results of the phonon density calculation, the interfacial thermal resistance of diamond/titanium carbide is low. -
随着微电子技术的不断发展，电子器件不断向高功率密度、高集成度发展，但器件因严重的发热问题，使常处于高温下工作的电子器件的稳定性和寿命大大降低[1]。为解决其散热问题，各种高热导率材料成为该领域研究的热点，其中的金属基金刚石复合材料是重要的研究方向。常用的金属基体有铝、银和铜，但银和铜与金刚石之间的润湿性较差，通常采用金属基体的合金化[2-3]对金刚石表面进行处理[4-6]和采用合理的复合材料成形工艺[7-9]改善两者之间的界面结合。
目前已经有大量试验对界面结合改善情况进行了研究。LI等[10]通过对金刚石铜表面镀钛得到了高热导率的金刚石/铜复合材料。WU等[11]通过真空热压烧结法制备金刚石/铜复合材料，并引入钛改善其界面结合，与不添加钛元素的复合材料相比，其热导率提升了48%。JHONG等[12]采用无压烧结法制备金刚石/银复合材料，并添加钛元素，制备出热导率高达953 W/（m·K）的复合材料，相当于其理论热导率的98%。TANG等[13]通过在金刚石表面镀覆碳化铬来改善金刚石与银之间的界面结合，并优化了试验中的工艺参数，得到了热导率为695 W/（m·K）的金刚石/银复合材料。
另外，可以通过第一性原理揭示材料界面结合能力和热传导问题。CHEN等[14]利用密度泛函理论和原子格林函数方法研究了金刚石/铜复合材料界面的纳米尺度的热运输问题，其声子谱的相关计算表明，金刚石与铜之间的声子耦合性较差，金刚石与铜之间的界面热阻较大。XIE等[15]讨论了钼掺杂对金刚石（111）/Al（Cu）（111）界面相互作用、拉伸性能和热导率的影响，表明铝与金刚石之间的界面结合要优于铜与金刚石的，钼掺杂后对铝/金刚石界面的增强作用有限，但对于铜/金刚石界面的增强作用较大。
上述的计算和试验研究均取得了一定的成果，但对金刚石/银的界面结构和界面性质的研究尚未见报道，并且通过引入碳化钛来改善金属与金刚石之间的界面结合和热传导能力的机理尚未有相关解释。因此，采用第一性原理计算的方法研究金刚石/银、金刚石/铜、金刚石/碳化钛 3 种界面的结构、电子结构和传热行为。
1. 计算方法与建模
1.1 计算方法
采用基于密度泛函理论的第一性原理方法，对3种界面结构的性质进行研究；采用Materials Studio中的Materials Visualizer模块进行模型构建；采用CASTEP模块进行结构优化[16]、能量计算以及界面性质的相关计算；采用FORCITE模块进行径向分布函数的计算。采用平面波展开的赝势为超软赝势[17]。采用广义梯度近似（generalized gradient approximation，GGA）中的PBE（perdew-becke-ernerhof）[18]泛函对交换关联能进行处理。布里渊区中k点的选取方法为Monkhorst-Pack网格[19]。收敛标准是系统总能量的变化在1×10−5 eV以内，结构优化后每个原子上的力小于0.03 eV/Å，晶胞的剩余应力小于0.05 GPa，公差偏移小于0.001 Å。在没有原子固定的情况下，采用BFGS（broyden-fletcher-goldfarb-shanno）[20]算法进行结构优化；在有原子固定的情况下，采用两点快速下降（two-point-steepest-descent,TPSD）[21]方法进行结构优化。结构优化和能量计算均采用的截断能为500 eV，k点网络为8×8×8；表面计算采用的截断能为400 eV，k点网络为8×8×1；界面模型最佳间距测试和结构优化采用的截断能为400 eV，k点网络为7×7×1。在处理表面结构和界面结构时，选取厚度为15 Å的真空层。电子态密度计算的k点网络为14×14×2。差分电荷密度是计算当前结构中每个原子及其孤立存在时的电子密度差异。
1.2 体相性质与表面性质
人造金刚石裸露的表面为（100）面和（111）面，而铜、银和碳化钛均属于立方晶系，面心立方布拉菲点阵，密排面均为｛111｝面。选用金刚石、铜、银和碳化钛的（111）面构建金刚石（111）/铜（111）、金刚石（111）/银（111）、金刚石（111）/碳化钛（111）的界面结构，揭示3种复合材料的界面性质。
为使界面模型中各种物质的晶体结构符合现实情况，对采用的各种晶胞模型进行结构优化，将晶胞模型的晶胞参数（A）与其他人的计算结果（B）和试验数据（C）进行比对。铜、银、金刚石、碳化钛的晶胞参数见表1。由表1可得：本文计算的晶胞参数满足计算所需要的精度要求。
表 1 计算的铜、银、金刚石、碳化钛的晶胞参数以及他人试验和计算的结果Table 1. Calculated cell parameters of copper, silver, diamond and titanium carbide, as well as the results of others’ experiments and calculations确定原子层数的方法为计算其表面结构的表面能。表面能计算的结果表明，6层金刚石（111）面、6层银（111）面、6层铜（111）面的表面能趋于稳定。碳化钛原子层数的确定采用测量弛豫后原子层间距的方法，其表面原子弛豫的结果如表2所示，结果表明9层碳化钛的内部原子层间距与块体碳化钛的原子层间距相似。因此，选择9层碳化钛原子，并且以钛原子作为第1层原子的表面结构构建其界面结构。
表 2 碳化钛的表面原子弛豫Table 2. Surface atomic relaxation of titanium carbide间距
位置层间距变化 ξ / % 第3层
原子第5层
原子第7层
原子第9层
原子第11层
原子1/2 −10.49 22.61 −21.21 −21.21 −20.36 2/3 5.39 11.02 11.98 12.45 3/4 −4.93 −7.10 −7.06 4/5 0.49 1.87 5/6 −2.34 1.3 界面的构建
由表1可以看出：只有铜与金刚石的晶胞参数相近。因此，需要重新定义碳化钛和银的表面结构的向量，由u（1, 0），v（0, 1）改为m（1, −1），n（2, 1）。铜与金刚石则进行扩胞使其变为原来的
$ \sqrt{2} $ 倍。之后，用式（1）计算金刚石/铜、金刚石/银、金刚石/碳化钛的3种晶格失配度[30]，分别为1.76%，0.52%，5.20%。经过表面重构之后，界面结构的失配度均较小，可以构成稳定的界面结构。界面处的几何参数以金刚石为基准，所构建的3种界面模型和界面处的原子重叠方式如图1所示。$$ \varDelta =\left|\frac{{a}_{1}-{a}_{2}}{{a}_{2}}\right| $$ (1) 其中：
$\varDelta $ 为晶格失配度，a1、a2分别为表面1和表面2晶格常数。2. 结果与讨论
2.1 界面间距和界面黏附功
采用通用结合能曲线（universal binding energy relation，UBER）测试图1中3种界面结构的界面间距，得到最为合理和稳定界面结构的UBER曲线，如图2所示。图2中的金刚石/铜黏附功最大时的界面结构间距最小，可能是在取界面间距时间隔过大，导致其数据有所偏差。
图 2 金刚石/铜、金刚石/银、金刚石/碳化钛界面结构的UBER曲线Figure 2. UBER curves of interface structures of diamond/copper, diamond/silver and diamond/titanium carbide黏附功计算公式如式（2）所示。在UBER曲线中，黏附功从大到小依次为金刚石/碳化钛的界面、金刚石/银的界面、金刚石/铜的界面。
$$ {W_{{\rm{ad}}}} = \frac{{{E_A}^{{\rm{slab}}} + {E_B}^{{\rm{slab}}} - {E_{A/B}}}}{S} $$ (2) 其中：EA/B为界面A/B的总能量，EAslab和EBslab分别表示构成界面A/B结构的表面结构的能量，S为界面结构的面积。
在固定非界面处原子的情况下，对于界面结构进行弛豫，并计算界面结构的黏附功，数据如表3所示。界面黏附功越大，界面间距越小，界面结合性能越强，界面结构越稳定。从表3的界面间距的角度分析可知：金刚石/银（2.236 Å）的界面间距较大，而金刚石/铜（2.046 Å）和金刚石/碳化钛（1.990 Å）的界面间距相近。这表明，金刚石/银之间的界面结合性能较差。而从表3的黏附功的角度分析，金刚石/碳化钛的界面黏附功（5.578 J/m2）相比于金刚石/铜的界面黏附功（1.919 J/m2）和金刚石/银的界面黏附功（4.291 J/m2）分别增加191%和30%。这说明，金刚石/碳化钛的界面结合最强，结构最稳定。
表 3 通过UBER和完全弛豫所得到的最佳界面间距和最佳黏附功Table 3. Optimal interface spacing and adhesion work obtained by UBER and complete relaxation methods
结构名称UBER 完全弛豫 界面间距
d / Å界面黏附功
Wad / (J·m−2)界面间距
d / Å界面黏附功
Wad / (J·m−2)金刚石/铜 1.901 2.302 2.046 1.919 金刚石/银 2.128 3.920 2.236 4.291 金刚石/碳化钛 1.967 4.947 1.990 5.578 2.2 电子态密度
图3为电子态密度计算结果。在3种界面结构中，界面附近的碳原子在费米能级处的态密度均不为0，即碳原子在金属原子的诱导下呈现出一定的金属性质。并且3种金属原子均在–20 eV附近出现新的峰值，表明金属原子与碳原子的轨道存在杂化作用。所有结构的第二层原子受界面结构的影响较小。
图 3 3种界面结构的电子态密度Figure 3. Electronic density of states of the three interface structures如图3a和3b所示：在金刚石/铜界面结构中界面处2种原子在−2 eV附近存在杂化峰，金刚石/银界面结构中界面处2种原子在−1 eV附近存在杂化峰。2个峰较为尖锐，表现出较强的局域化特性。但2种结构中杂化峰数量较少，表明铜原子和银原子与碳原子之间相互作用较弱，即2种界面结构的结合能力较低。
如图3c所示：在金刚石/碳化钛的界面结构中，界面处钛原子和碳原子的波形发生了较大的变化。碳原子的态密度向费米能级移动，并且在−5～5 eV处与钛原子的s和p轨道出现了多处共振峰和多处重叠区域。表面钛原子与碳原子之间相互作用较强，这种作用增大了界面结构的结合能力。
2.3 马利肯布居分析
3种界面结构界面处原子轨道布居分析如表4、表5、表6所示。所取的原子为2种表面结构中最靠近界面处的2层原子，一共4层原子。电荷的转移数越多，形成的键合作用越强。在金刚石/铜界面处，铜原子失去0.16 eV，碳原子得到0.20 eV。在金刚石/银界面处，银原子失去0.22 eV，碳原子得到0.16 eV。这2种界面结构主要为金属原子的p轨道失去电子，碳原子的s和p轨道得到电子。综合电子得失，金刚石/银界面键合作用更强。在金刚石/碳化钛的界面结构中，碳原子得到0.27 eV，s和p轨道均得到电子。这种结构的碳原子得到的电荷数目为3种界面结构中最多，表明这种界面结构键合作用最强，界面结构最稳定。
表 4 金刚石/铜界面原子轨道布居分析Table 4. Analysis of atomic orbital population at diamond/copper interface原子种类及位置 电荷量 Q / eV s p d 总数 差值 铜第2层铜原子 0.49 0.75 9.69 10.93 0.07 铜第1层铜原子 0.54 0.60 9.70 10.84 0.16 金刚石第1层碳原子 1.20 3.00 0.00 4.20 −0.20 金刚石第2层碳原子 1.14 2.92 0.00 4.06 −0.06 表 5 金刚石/银界面原子轨道布居分析Table 5. Analysis of atomic orbital population at diamond/silver interface原子种类及位置 电荷量 Q / eV s p d 总数 差值 银第2层银原子 0.59 0.58 9.78 10.95 0.05 银第1层银原子 0.60 0.43 9.75 10.78 0.22 金刚石第1层碳原子 1.18 2.98 0.00 4.16 −0.16 金刚石第2层碳原子 1.13 2.93 0.00 4.06 −0.06 表 6 金刚石/碳化钛界面原子轨道布居分析Table 6. Atomic orbital population analysis of diamond/titanium carbide interface原子种类及位置 电荷量 Q / eV s p d 总数 差值 碳化钛第2层碳原子 1.47 3.25 0.00 4.72 −0.72 碳化钛第1层钛原子 2.18 6.36 2.64 11.19 0.81 金刚石第1层碳原子 1.20 3.08 0.00 4.27 −0.27 金刚石第2层碳原子 1.15 2.94 0.00 4.08 −0.08 2.4 差分电荷密度分析
界面结构的原子差分电荷密度，如图4a、图4b、图4c所示。在3种结构中的金刚石部分，碳原子之间存在强烈的电荷增加，并呈现球形，显示出强共价键作用的特点。在金刚石/铜的界面处仅有极少部分的黄色区域，表示金刚石表面的原子仅得到极少部分的电荷，并且其黄色区域较为分散，没有固定的取向性。而在金刚石/银的界面处，黄色区域更加集中，并且位于银原子的下方，呈现出一定的方向性，这表面银原子与碳原子之间不仅存在一定量的电荷转移，并且有一定的键合作用。这种键合作用使得金刚石/银的界面结构更加稳定。而碳化钛与金刚石界面处的电荷转移则更为强烈，其颜色为红色并且更加集中，且位于钛原子的下方。这表明：钛原子与碳原子之间具有更多的电荷转移和更强的键合作用。
为了更清晰地表现这3种界面结构的电荷转移和成键强弱，取金刚石表面正上方0.3 Å处的差分电荷密度图，如图4d、图4e、图4f所示。3种截面处的电荷区域的形状均为哑铃形的截面，说明均为碳原子的p轨道得到了电荷。图中颜色的深浅反映了电荷转移量的大小和成键的强弱，其中金刚石/铜界面结构为淡黄色，表现出极少的电荷转移和极弱的键合作用。而金刚石/碳化钛界面处则有红色部分出现，表现出强烈的电荷转移和强的键合作用。金刚石/银界面处的电荷转移情况则位于两者中间，即一般的电荷转移量和较强的键合作用。综上，金刚石/碳化钛界面结构电荷转移最多，成键作用明显，界面结构最稳定，界面结合最强。
2.5 径向分布函数
3种界面结构的径向分布函数如图5所示。所取的参考原子为金刚石表面结构的第1层中处于中心位置的碳原子。只考虑5 Å以内的径向分布函数。在金刚石表面结构中碳原子与碳原子的距离有1.511 Å，2.493 Å，2.517 Å，2.923 Å，4.365 Å等，其对应的峰均在图5出现。3种界面结构的界面间距均在2 Å左右（见表3），在图5中只有碳化钛与金刚石的界面结构在2 Å附近处出现了一个峰，表现出碳原子与钛原子存在相互作用，即两者表面的结合强度大。而在金刚石/铜与金刚石/银的图中并没有新的峰出现，表明碳原子与铜原子、碳原子与银原子间没有相互作用，即金刚石/银和金刚石/铜的结合强度较低。。
2.6 声子态密度分析
金刚石、银、铜和碳化钛的声子态密度如图6所示。从图6可以看出：金刚石声子的频率分布主要为15～40 THz，其峰值在37 THz左右。而银与铜的声子频率分布与金刚石差距较大，这表明：金刚石与铜和银之间的声子耦合性能均较差，导致了金刚石与铜和银界面处的界面热阻较大，说明只采用金属和金刚石所制备的复合材料无法完全发挥金刚石热导率高的特点。而碳化钛的声子频率主要在0～13 THz和15～20 THz。其中，较高的频率分布部分与金刚石的声子频率在15～21 THz，接近完全重合，这部分频率的声子存在使得金刚石与碳化钛之间存在一定的声子耦合性能，使其界面之间的热量传输能力有一定提高，从而提升了复合材料的热导率。
3. 结论
（1）采用第一性原理计算分析了金刚石/银、金刚石/铜、金刚石/碳化钛的界面结构。相比于金刚石/铜和金刚石/银，金刚石/碳化钛的界面间距最小（1.990 Å），界面黏附功最大（5.578 J/m2），界面黏附功分别提高191%和30%。这说明金刚石/碳化钛的界面结合最强，界面结构最稳定。
（2）通过电子态密度，马利肯布居，差分电荷密度分析界面处的键合情况。金刚石/铜界面处存在电荷转移，但键合作用极弱。金刚石/银界面处存在键合作用，但强度不高。金刚石/碳化钛界面存在较多的电荷转移和较强的键合作用。
（3）通过计算声子态密度分析界面处声子耦合情况，结果表明：金刚石/铜和金刚石/银界面的声子耦合性均较差，金刚石与铜和银之间的界面热阻较大。金刚石/碳化钛界面之间声子耦合性较好，表明碳化钛的引入降低了界面热阻，提升了复合材料的热导率。
-
图 2 金刚石/铜、金刚石/银、金刚石/碳化钛界面结构的UBER曲线
Figure 2. UBER curves of interface structures of diamond/copper, diamond/silver and diamond/titanium carbide
图 3 3种界面结构的电子态密度
Figure 3. Electronic density of states of the three interface structures
表 1 计算的铜、银、金刚石、碳化钛的晶胞参数以及他人试验和计算的结果
Table 1. Calculated cell parameters of copper, silver, diamond and titanium carbide, as well as the results of others’ experiments and calculations
下载: 导出CSV
表 2 碳化钛的表面原子弛豫
Table 2. Surface atomic relaxation of titanium carbide
间距
位置层间距变化 ξ / % 第3层
原子第5层
原子第7层
原子第9层
原子第11层
原子1/2 −10.49 22.61 −21.21 −21.21 −20.36 2/3 5.39 11.02 11.98 12.45 3/4 −4.93 −7.10 −7.06 4/5 0.49 1.87 5/6 −2.34
下载: 导出CSV
表 3 通过UBER和完全弛豫所得到的最佳界面间距和最佳黏附功
Table 3. Optimal interface spacing and adhesion work obtained by UBER and complete relaxation methods
结构名称UBER 完全弛豫 界面间距
d / Å界面黏附功
Wad / (J·m−2)界面间距
d / Å界面黏附功
Wad / (J·m−2)金刚石/铜 1.901 2.302 2.046 1.919 金刚石/银 2.128 3.920 2.236 4.291 金刚石/碳化钛 1.967 4.947 1.990 5.578
下载: 导出CSV
表 4 金刚石/铜界面原子轨道布居分析
Table 4. Analysis of atomic orbital population at diamond/copper interface
原子种类及位置 电荷量 Q / eV s p d 总数 差值 铜第2层铜原子 0.49 0.75 9.69 10.93 0.07 铜第1层铜原子 0.54 0.60 9.70 10.84 0.16 金刚石第1层碳原子 1.20 3.00 0.00 4.20 −0.20 金刚石第2层碳原子 1.14 2.92 0.00 4.06 −0.06
下载: 导出CSV
表 5 金刚石/银界面原子轨道布居分析
Table 5. Analysis of atomic orbital population at diamond/silver interface
原子种类及位置 电荷量 Q / eV s p d 总数 差值 银第2层银原子 0.59 0.58 9.78 10.95 0.05 银第1层银原子 0.60 0.43 9.75 10.78 0.22 金刚石第1层碳原子 1.18 2.98 0.00 4.16 −0.16 金刚石第2层碳原子 1.13 2.93 0.00 4.06 −0.06
下载: 导出CSV
表 6 金刚石/碳化钛界面原子轨道布居分析
Table 6. Atomic orbital population analysis of diamond/titanium carbide interface
原子种类及位置 电荷量 Q / eV s p d 总数 差值 碳化钛第2层碳原子 1.47 3.25 0.00 4.72 −0.72 碳化钛第1层钛原子 2.18 6.36 2.64 11.19 0.81 金刚石第1层碳原子 1.20 3.08 0.00 4.27 −0.27 金刚石第2层碳原子 1.15 2.94 0.00 4.08 −0.08
下载: 导出CSV
-
[1] MOORE A L, SHI L. Emerging challenges and materials for thermal management of electronics [J]. Materials Today,2014,17(4):163-174. doi: 10.1016/j.mattod.2014.04.003 [2] AZINA C, CORNU I, SILVAIN J F, et al. Effect of titanium and zirconium carbide interphases on the thermal conductivity and interfacial heat transfers in copper/diamond composite materials [J]. AIP Advances,2019,9(5):055315. doi: 10.1063/1.5052307 [3] BAI G, ZHANG Y, DAI J, et al. Tunable coefficient of thermal expansion of Cu-B/diamond composites prepared by gas pressure infiltration [J]. Journal of Alloys and Compounds,2019,794:473-481. doi: 10.1016/j.jallcom.2019.04.252 [4] JIA J, BAI S, XIONG D, et al. Effect of tungsten based coating characteristics on microstructure and thermal conductivity of diamond/Cu composites prepared by pressureless infiltration [J]. Ceramics International,2019,45(8):10810-10818. doi: 10.1016/j.ceramint.2019.02.156 [5] MOLINA-JORDÁ J M. Thermal conductivity of metal matrix composites with coated inclusions: A new modelling approach for interface engineering design in thermal management [J]. Journal of Alloys and Compounds,2018,745:849-855. doi: 10.1016/j.jallcom.2018.02.092 [6] CHE Z, LI J, WANG Q, et al. The formation of atomic-level interfacial layer and its effect on thermal conductivity of W-coated diamond particles reinforced Al matrix composites [J]. Composites Part A: Applied Science and Manufacturing,2018,107:164-170. doi: 10.1016/j.compositesa.2018.01.002 [7] CIUPIŃSKI Ł, KRUSZEWSKI M J, GRZONKA J, et al. Design of interfacial Cr3C2 carbide layer via optimization of sintering parameters used to fabricate copper/diamond composites for thermal management applications [J]. Materials & Design,2017,120:170-185. [8] CHEN H, JIA C, LI S. Interfacial characterization and thermal conductivity of diamond/Cu composites prepared by two HPHT techniques [J]. Journal of Materials Science,2012,47(7):3367-3375. doi: 10.1007/s10853-011-6180-6 [9] ARAI S, UEDA M. Fabrication of high thermal conductivity copper/diamond composites by electrodeposition under potentiostatic conditions [J]. Journal of Applied Electrochemistry,2020,50(5):631-638. doi: 10.1007/s10800-020-01414-3 [10] LI H, WANG C, WU L, et al. Optimization of process parameters, microstructure, and thermal conductivity properties of Ti-coated diamond/copper composites prepared by spark plasma sintering [J]. Journal of Materials Science: Materials in Electronics,2021,32(7):9115-9125. doi: 10.1007/s10854-021-05579-1 [11] WU X, WAN D, ZHANG W, et al. Constructing efficient heat transfer channels at the interface of diamond/Cu composites [J]. Composite Interfaces,2021,28(6):625-635. doi: 10.1080/09276440.2020.1795466 [12] JHONG Y S, TSENG H T, LIN S J. Diamond/Ag-Ti composites with high thermal conductivity and excellent thermal cycling performance fabricated by pressureless sintering [J]. Journal of Alloys and Compounds,2019,801:589-595. doi: 10.1016/j.jallcom.2019.06.167 [13] TANG Y, WANG L, ZHAO C. Enhancement of the thermal properties of silver-diamond composites with chromium carbide coating [J]. Applied Physics A,2014,115(2):379-385. doi: 10.1007/s00339-014-8254-1 [14] CHEN L, CHEN S, HOU Y. Understanding the thermal conductivity of diamond/copper composites by first-principles calculations [J]. Carbon,2019,148:249-257. doi: 10.1016/j.carbon.2019.03.051 [15] XIE H, CHEN Y, ZHANG T, et al. Adhesion, bonding and mechanical properties of Mo doped diamond/Al (Cu) interfaces: A first principles study [J]. Applied Surface Science,2020,527:146817. doi: 10.1016/j.apsusc.2020.146817 [16] SEGALL M D, LINDAN P J D, PROBERT M J, et al. First-principles simulation: Ideas, illustrations and the CASTEP code [J]. Journal of Physics: Condensed Matter,2002,14(11):2717-2744. doi: 10.1088/0953-8984/14/11/301 [17] VANDERBILT D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism [J]. Physical Review B,1990,41(11):7892. doi: 10.1103/PhysRevB.41.7892 [18] WU Z, COHEN R E. More accurate generalized gradient approximation for solids [J]. Physical Review B,2006,73(23):235116. doi: 10.1103/PhysRevB.73.235116 [19] CHADI D J. Special points for Brillouin-zone integrations [J]. Physical Review B,1977,16(4):1746. doi: 10.1103/PhysRevB.16.1746 [20] HEAD J D, ZERNER M C. A broyden-fletcher-goldfarb-shanno optimization procedure for molecular geometries [J]. Chemical Physics Letters,1985,122(3):264-270. doi: 10.1016/0009-2614(85)80574-1 [21] BARZILAI J, BORWEIN J M. Two-point step size gradient methods [J]. IMA Journal of Numerical Analysis,1988,8(1):141-148. doi: 10.1093/imanum/8.1.141 [22] KIM W J, XING S, KREMER G, et al. Electronic structure of heavy halogen atoms adsorbed on the Cu (111) surface: A combined ARPES and first principles calculations study [J]. The Journal of Physical Chemistry C,2019,123(43):26309-26314. doi: 10.1021/acs.jpcc.9b07057 [23] VILLARS P, DAAMS J L C. Atomic-environment classification of the chemical elements [J]. Journal of Alloys and Compounds,1993,197(2):177-196. doi: 10.1016/0925-8388(93)90041-K [24] QIAN H Q, MAO H Y, SONG F, et al. Study on the adsorption of fluorescein on Ag (110) substrate [J]. Applied Surface Science,2010,256(9):2686-2690. doi: 10.1016/j.apsusc.2009.10.056 [25] LIU L, Bassett W A. Compression of Ag and phase transformation of NaCl [J]. Journal of Applied Physics,1973,44(4):1475-1479. doi: 10.1063/1.1662396 [26] GU K, PANG M, ZHAN Y. Insight into interfacial structure and bonding nature of diamond (001)/Cr3C2 (001) interface [J]. Journal of Alloys and Compounds,2019,770:82-89. doi: 10.1016/j.jallcom.2018.08.112 [27] OWNBY P D, YANG X, LIU J. Calculated X‐ray diffraction data for diamond polytypes [J]. Journal of the American Ceramic Society,1992,75(7):1876-1883. doi: 10.1111/j.1151-2916.1992.tb07211.x [28] KHANZADEH M, ALAHYARIZADEH G. A DFT study on pressure dependency of TiC and ZrC properties: Interconnecting elastic constants, thermodynamic, and mechanical properties [J]. Ceramics International,2021,47(7):9990-10005. doi: 10.1016/j.ceramint.2020.12.145 [29] RAZUMOVSKIY V I, POPOV M N, DING H, et al. Formation and interaction of point defects in group IVb transition metal carbides and nitrides [J]. Computational Materials Science,2015,104:147-154. doi: 10.1016/j.commatsci.2015.03.042 [30] SUN Q, WEI Y, ZHANG J, et al. Effect of lattice mismatch on HgCdTe LPE film surface morphology [J]. Journal of Electronic Materials,2016,45(9):4674-4679. doi: 10.1007/s11664-016-4637-8 -
下载:
点击查看大图
计量
- 文章访问数: 885
- HTML全文浏览量: 355
- PDF下载量: 109
- 被引次数: 0
邮件订阅
RSS
