Abstract
本文的阵列设计Ç84 -嵌入式Si衬底捏造使用受控的自组装方法中的超高真空室的报告。的C- 84的特性-嵌入式Si表面,如原子分辨率地形,态本地电子密度,带隙能量,场发射性能,纳米机械刚度,以及表面磁性,使用各种表面分析技术下超进行了检查,在高真空(UHV)条件下,以及在一个大气压的系统。实验结果表明的C- 84的高度均匀性-嵌入式Si表面捏造使用受控的自组装纳米机制,表示场发射显示器(FED),光电器件制造,微机电系统的应用的一个重要发展切削工具,并在努力找到硬质合金半导体合适的替代者。分子动力学(MD)方法,半经验潜力b用于研究的ç84纳米压痕ê -嵌入式硅衬底。执行MD模拟详细描述这里提出。对MD模拟的力学分析,如压痕力,杨氏模量,表面硬度,原子应力,应变原子的综合研究的细节都包括在内。压痕模型的基本应力和冯·米塞斯应变分布可以计算出监测与原子级别时间评估变形机制。
Introduction
富勒烯分子和它们包含有纳米材料之间独特的,由于其优异的结构特点,电子传导性,机械强度,化学性能1-4复合材料。这些材料已被证明在一个范围内的领域,如电子,电脑,燃料电池技术,太阳能电池,以及场致发射技术5,6非常有益的。
在这些材料中,碳化硅(SiC)纳米粒子的复合材料已经受到了特别的关注由于它们的宽带隙,高导热性和稳定性,高电击穿能力,和化学惰性。这些好处是在光电器件特别明显,金属氧化物半导体场效应晶体管(MOSFET),发光二极管(LED),和高功率,高频率,和高温应用。然而,高密度的缺陷通常conventi的表面上观察到的Onal地区碳化硅可以对电子结构的不利影响,甚至导致器件失效7,8。尽管碳化硅的应用程序已自1960年以来研究的事实,这个特殊的未解决的问题仍然存在。
本研究的目的是一个C 84的制造-嵌入式Si衬底异质结和随后的分析,得到所得到的材料的电子,光电,机械,磁性,和场发射性能的全面理解。我们还讨论通过数值模拟来预测纳米材料的特性,通过分子动力学计算得到新的应用的问题。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
注:本文概述在一个半导体衬底的表面上的自组装富勒烯阵列的形成中使用的方法。具体地讲,我们提出了用作微电子机械系统(MEMS),以及在高温,高功率光电器件,应用的场致发射体或基底富勒烯嵌入硅衬底的制备的新方法,以及在高-频率器件9-13。
在Si衬底1的制作六角形封闭封装(HCP)加铺结构Ç84
- 准备清洁的Si(111)基板
- 主体的Si衬底的RCA清洗(美国无线电公司),涉及随后加热在超高真空系统中除去从基板的表面上的氧化物层和杂质的溶剂的应用(见辅助材料)。
注意:在这里,术语“UHV-超高真空系统”是指到低于1×10的真空-8帕中的Si(111)的制备中使用。
- 主体的Si衬底的RCA清洗(美国无线电公司),涉及随后加热在超高真空系统中除去从基板的表面上的氧化物层和杂质的溶剂的应用(见辅助材料)。
- 存款ç84硅表面以特高压系统使用热蒸发
- 通过加热丝至500℃的预加热用外部电源的K细胞蒸发器,以促进杂质的放气。
- 载荷C 84纳米到K-电池容器。电阻加热的K细胞至650℃。汽化ç84纳米粒子为C 84纳米粒子在容器组成蒸汽。蒸发下在直线84纳米至纳米颗粒通过一个控制阀在低于5×10 -8帕压力取得一个Si衬底。
- 嵌入Ç84通过自组装机理Si表面分子内
- 在一个超高真空系统的预退火的Si(111)基片,在900℃,得到(1×1)结构。降低温度至650℃用于该depos 30分钟在衬底的表面的C 84纳米的银行足球比赛。
- 在〜750退火在Si基板 ℃下12小时,在此期间,粉状-C的84纳米自组装成在Si(111)基板的表面上的高度均匀的富勒烯阵列。
注意:在这里,术语“高度均匀的富勒烯阵列”是指富勒烯的基板上的均匀分布,其中大部分的纳米颗粒在垂直于衬底的表面取向在一个紧凑的布置。这种配置有助于确保富勒烯阵列的垂直高度为所有样本中基本上是相同的。
2.Ç84电子性能的测量-嵌入式硅衬底
- 衡量各国的本地电子密度使用超高真空扫描隧道显微镜
- 测量使用UHV-SPM特定原子的IV曲线
- 将C 84 子>在SPM样品架 - 嵌入式硅衬底。引进的持有人到UHV-STM(扫描隧道显微镜)扫描头系统。扫描应用样本偏差-5 V至5 V
- 点击“IV”测量项目测量隧道电流我在原子分辨率。选择进行测量的C 84 -嵌入式Si衬底上至少20的特定位置。计算隧道电流I 20的特定位置的平均值。派生我作为电压的一个函数。情节IV曲线。
- 对于V.计算的I(V)的衍生物转换IV曲线到的dI / DV作为电压的函数,以便确定的C- 84的本地电子状态-嵌入式Si衬底。
- 测量带隙能量
- 根据从以下2.1.2和2.1.3的程序获得IV曲线:SI(111)-7x7表面的Si(111)-1x1表面,单一个体Ç84纳米粒子在Si,7-19Ç84在Si集群,20-50Ç84簇在Si和嵌入式硅表面内Ç84单层。
- 对于V.计算的I(V)的衍生物转换IV曲线到的dI / DV曲线来测量的HOMO-LUMO能级差异(称为带隙能量)在每个测量位置, 如图2a中 。
- 获取场发射(FE)属性
- 将C 84在FE样品架-嵌入式硅衬底。将支架插入有限元分析室。抽空腔室,以约5×10 -5 Pa的是FE测量的压力。
注意:C 84 -嵌入式硅衬底充当阴极和与〜0.71毫米2充当阳极的横截面面积的铜探针。在阴极和阳极之间的距离为大约590微米。 - 在衬底手动增加施加电压从100伏到1,100 V.测量对应ING场发射电流作为使用高电压源测量单元与电流放大器一个施加电压的功能。
- 根据如图2b所示的功函数〜5 eV的计算福勒-诺德海姆场发射的相关性。
- 具有约4,383一个β值F(场)=β(V / D):得到的几何场增强因子(β)如下。
- 获取基于所述自然对数的斜率(J / E 2)与(1 / E),这给了我们的〜4.0×10 6 V / cm的为C 84 -嵌入式Si衬底的值在真空下的击穿电场如图2c所示。
- 将C 84在FE样品架-嵌入式硅衬底。将支架插入有限元分析室。抽空腔室,以约5×10 -5 Pa的是FE测量的压力。
- 光电性能
- 传输测试衬底的光学发射测量系统。注重与位于样品室的中心在基板上325纳米的排放He-Cd激光器光源。设置在适当的位置上的光谱仪。使用SPEctrometer通过收集和分析发射光子以获得光致发光光谱。光电结果示于图2d。
3.表面磁性测量
- 获得MFM(磁力显微镜)地形。
- 磁化的Ç84样品-嵌入式硅前MFM测量通过施加磁铁用约2千奥斯特的磁场强度。
- 广场上的MFM样品台的磁化样本。点击“获取MFM地志”项目。观察在嵌入升降模式中使用MFM与磁化垂直于样品的表面上的应用在硅衬底内的磁畴的富勒烯的微观结构。
- 使用MFM测量纳米级PPP-MFMR悬臂( 图3a)。确定表面磁性如果MFM地形显得较暗(亮)时尖的磁矩是在SAM基材时刻E(对面)的方向。
- SQUID(超导量子干涉器件)测量
- 准备Ç84单层-嵌入式基于C 84嵌入硅衬底硅衬底和C 84簇。
- 磁化的ç84样品嵌入硅衬底之前SQUID实验通过应用磁铁约2千奥斯特的磁场强度-嵌入式Si和C 84簇基于C 84。
- 将在SQUID样品。在一系列〜2千奥斯特应用扫描磁场。获得作图在室温下于SQUID测量的外部磁场的磁化循环。
注:如如图3b中所示为铁磁性材料的典型MH曲线可以得到。
4. AFM纳米力学性能的测量
注:原子力显微镜(AFM)提供在微观和材料和机械性能在空气中的纳米尺度的特征,以及在一个UHV环境强大的工具
- 衡量Ç84嵌入式硅衬底在大气条件下的刚度
- 上放置的AFM样品台的基板。将一个锋利的尖端在使用扫描仪的基材。监视尖端的位移作为针尖 - 样品的相互作用力的量度。点击“力测量”项记录在某一位置沿垂直方向的许多针尖样的距离运动。
- 获得使用来自天然氧化物的2-3纳米的层一个RCA清洗Si基板以及从C 84大气压条件下的AFM测力-嵌入式Si衬底和涂覆有碳化硅的薄膜的Si基片。
- 使用AFM软件,绘制大气条件下力 - 距离曲线。
注:AFM悬臂是硅探针的尖端半径〜〜的40牛顿/米5-20纳米和弹簧常数。
- 衡量Ç84的刚度嵌入式硅基特高压室
- 根据4.1.1在从一个RCA清洗Si基板,一个干净的Si(111)-7x7表面的超高真空系统中使用的AFM指导获得力的测量,一个C 84 -嵌入式Si衬底,衬底和一个硅衬底涂覆的与碳化硅的薄膜。
- 剧情以特高压系统力 - 距离曲线。注:AFM悬臂是硅探头的〜5-20 nm和弹簧常数〜40牛顿/米的尖端半径图4给出无序硅表面的力-距离分析,7×7表面的,独立的自章C 84嵌入Si表面和Si表面组装内层,如使用UHV-AFM决定。
通过分子动力学模拟纳米力学性能的5.测量
注:在模拟部分,OVITO 16(开源visualizati上的软件),并且OSSD 17(开放面结构数据库)用于创建仿真模型和结果可视化。 LAMMPS 14(一个开源的分子动力学(MD)模拟包)被用来进行纳米压痕模拟和分析模拟结果15。所有仿真作业是用在先进的大型并行超星系的NCHC(ALPS),并行计算执行。
注意:要使用MD模拟研究的C 84单层/ Si衬底异质结,应准备几个步骤,一个模拟模型来获得嵌入到Si衬底严格的C 84单层。需要注意的是它难以产生一个恰好从实验数据相同的结构,因为该间结构的Ç84单层和Si(111)衬底异质结之间的复杂的。其结果是,我们使用一个人工的方式与程序的若干步骤,生成模拟模型,其在图5中示出。细节在以下协议说明。我们描述的MD在LAMMPS参数如何设置,建立嵌入到衬底的严格的C 84富勒烯单层,执行一个缩进过程,并分析模拟结果。
- 在LAMMPS输入文件参数设置
- 使用边界命令设置在x周期性边界条件和y方向。
- 用“固定速度”命令给初始速度与系统的各原子的高斯分布分配,随机。
- 使用“修复pair_style”命令来指定Tersoff 18和19 AIREBO潜力来形容的Si-Si和Si-C相互作用和CC的互动,分别为。
- 使用“修正NVT”和“固定的npt”命令来采取鼻胡佛方法20,以确保系统保持在所需的温度和压力下,以基因评价规范和等温等压系综20,在该系统中,速度- verlet的算法20被用来预测原子的轨迹。同时使用“修复NVT”和“运行”命令来设置3 K /皮秒退火过程的冷却速度。
- 使用“时间步”命令来设置0.2 FSEC的时间积分的时间步长。
- 使用“修复墙/体现”命令采用反射墙来限制自由(5.3.2)的程度。
- 使用“区”和“基团”的基片成不同的控制层(5.4.3)分成:牛顿原子层,热控制层和底部固定层,其可以通过使用“修正NVE”成立,“修复NVT“和”修复setforce“命令,分别。
- 使用“区域”和“create_atoms”命令来创建一个球形探头。
- 使用“修复移动”命令来嵌入C84单分子层到基板(5.4。2)与模拟(5.5.2)中移动探针。
- 使用“运行”命令来执行分子动力学模拟。
- 使用“计算力”(5.6.1)和“计算应力/原子”(5.6.4)命令来评价原子应力和压痕力。
注意:在下文中,除了结构建立,所有的步骤被LAMMPS脚本完成的。
- 使用OSSD和OVITO到硅(111)7×7表面的制备。
- 打开OSSD软件。点击“搜索”按钮。 “搜索条件”面板呈现。之所以选择Si衬底,元素类型,重建的结构,半导体ELEC,菱形格,111脸和7×7的模式。点击“搜索”和“接受”按钮。 “结构列表”面板呈现。单击所需的结构( 即 ,硅(111)7×7)。点击“文件”按钮。保存协调文件名为.xyz文件。
注:我们指出的结构从OSSD抽取的数据库是不是我们的压痕模拟足够大。其结果是,我们通过以下步骤重建较大和较厚衬底。 - 打开OVITO软件。加载名为.xyz文件到OVITO。使用“切片”命令来捕获所述Si的(111)7×7表面在x和y方向26.878点¯x46.554埃2的尺寸的超晶胞。导出的数据文件。使用“切片”命令来捕捉底部的Si(111)衬底与3的26.878点¯x46.554点¯x9.7大小的超级。使用“查看定期图像”命令来复制超级12倍z方向。导出的数据文件。
- 结合Si的数据文件(111)7×7表面和Si(111)衬底模型通过记事本+ +(一个免费的源代码编辑器)。最后,将合并的数据加载到OVITO。使用“查看周期性图像”复制在x和y方向上具有5×3超元放大基板的大小。
- 使用LAMMPS执行20微秒MD模拟时间放宽仿真模型。在下文中,1550ķ执行淬火处理至室温以500微秒仿真时间。最后,执行一个10皮秒的模拟时间进行最后的放松过程。
- 打开OSSD软件。点击“搜索”按钮。 “搜索条件”面板呈现。之所以选择Si衬底,元素类型,重建的结构,半导体ELEC,菱形格,111脸和7×7的模式。点击“搜索”和“接受”按钮。 “结构列表”面板呈现。单击所需的结构( 即 ,硅(111)7×7)。点击“文件”按钮。保存协调文件名为.xyz文件。
- 章C 84富勒烯单分子膜制备方法
- 从网络21下载C 84富勒烯的结构优化的协调文件,并写入FORTRAN程序复制安排在一个蜂窝结构49℃84富勒烯。
- 使用LAMMPS设置反映壁在与下面的C 84单层以确保该分子留在一个计划。执行MD模拟时间为200皮秒放松仿真模型。在下文中,执行从700 K A淬火过程到室温,以获得500微秒仿真时间水珠最小状态。最后,执行一个10皮秒的模拟时间进行最后的放松过程。
- 建立Indenta在硅ç84富勒烯单分子膜化模型(111)7×7表面。
- 写FORTRAN代码躺在硅的C 84单层(111)7×7表面3的距离,建立压痕模型。
- 使用LAMMPS嵌入的C 84单层与2〜3埃的深度的基板。在下文中,运行一个40微秒的模拟时间为系统的放松。最后,退火系统至室温。
- 划分硅衬底成上层牛顿原子层,热控制层和底部固定层,分别是0.7,2,和厚度5.3纳米。的C 84单层也建模为牛顿原子。
- MD压痕工艺
- 使用LAMMPS与基于C 84 / SI(111)7×7表面的模式直径5纳米( 图5)创建一个球形探头。探头设定为刚体。指定的PR 10米/秒的恒定速度OBE向下朝向缩进过程试样移动。
- 向下移动探针与样品以恒定的速度,直到特定装载深度( 即,包括1.5的情况下,2.5,4.5,10,15,20,和30埃,以便探索的C 84富勒烯单层的效果在装载过程中的Si衬底,其中C 84富勒烯的尺寸为11埃)上。持探针在夹持过程基板以允许原子的松弛。最后,在缩回过程的恒定速度提取从该衬底的探针。
- 计算与分析
- 总结根据下列公式在探针原子的垂直力计算压痕力:
(1) - 萃取从压痕的力 - 距离曲线的还原模量和硬度。基于奥利弗和法尔&#39氏方法22,线性关系可以的杨氏模量和卸载刚度之间的。刚度( 即 ,初始部分的斜率)的卸载曲线的定义为
(2)
其中,P,H,A和E r为压痕负荷,探针的弹性位移,压痕投影面积,降低了模量。 β(= 1为圆形压头)是形状修正因子。降低模量和杨氏模量之间的关系可以写成
(3)
其中E和v是杨氏模量和泊松比的样品和E I和V i是杨氏模量和泊松比的压头。 - 的H定义计算硬度= P 马X / A,其中P 最大和A是最大压痕力和探头的投影面积。
- 计算由在所述n -方向基板的m面的维里原子应力22
(4)
其中m i为原子i的质量; 和是在第m原子i的速度分量-和n -directions分别; V I是围绕原子我分配体积; N s是包含区域S,其中S是指原子的相互作用的区域内的粒子数;Φ(R IJ)是潜在的功能; R ij是原子i和之间的距离和就在M -从原子到我的Atom j中的向量的n个晶向组件。 - 使用OVITO表示根据下列公式每个原子不变的冯米塞斯应力:
(5)
- 总结根据下列公式在探针原子的垂直力计算压痕力:
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
Ç84分子的无序的Si(111)面A单层使用受控的自组装过程中一个超高真空腔室制造图1示出了一系列通过UHV-STM不同程度的覆盖率的测量地形图像:(一) 0.01毫升,(b)0.2毫升,(C)0.7毫升,和(d)0.9毫升。的C 84嵌入Si衬底的电子和光学性质使用各种表面分析技术,如STM和PL( 图2)还研究。得到的样品的优异材料性质证明纳米技术可以如何用于物质在atomic-和纳米尺度上的控制。 图3中的MFM和SQUID结果表明Ç84嵌入基片的表面的磁性。 图4给出了参考提出衬底的纳米机械的UHV-AFM结果。我们的实验结果表明,中C 84的潜力图4中)。
在模拟部分,所有的程序都用LAMMPS进行MD模拟完成。富勒烯埋置基板的机械性能(压痕力和接触应力)的计算,并在图6中示出。在不同的时间步原子的冯米塞斯应变分析用于表征局部变形。相应的快照作为压痕深度的函数可以在图6的插入件,其分别计算与由OVITO可视看到。的压痕力作为压痕深度的函数的结果被用来计算硬度H( 图7a),减少的模量E <子> R( 图7b),并装载刚度S的C 84单层( 图8)。结果可以与由实验确定被比较和提供的观点更详细地点来解释的机械性能的变化。
图1:C 84 嵌入式Si衬底具有不同的覆盖的STM地形图(40×40纳米2)表示ç吸附在Si 84分子的系列(111)面以2伏的负样本偏置,如通过UHV-STM测量。不同程度的覆盖范围:(1)0.01毫升,(b)0.2毫升,(C)0.7毫升,和(d)0.9毫升。
图2: 电子属性第C 84 测量 嵌入式Si衬底 ( 一个 )IV曲线和差分衍生物电导(DI / DV)与章C 84单自组装层的电压曲线,如通过UHV-STM确定。 ; ( 二 )场发射电流密度VS电场曲线; ( 三 )对应的表面的FN情节与嵌入式C 84,如使用源测量单元测量; ( 四 )Ç84单自组装层的光致发光光谱。重新打印与基准12的许可。 请点击此处查看该图的放大版本。
图3: 第C 84 表面磁性 嵌入式Si衬底 ( 一个 )嵌入用C 84 Si基的MFM图像。 ( 二 )磁化循环暗算外部磁场请点击此处查看该图的放大版本。
图4:第C 84纳米机械调查嵌入式Si衬底无序的Si表面的力-距离分析,7×7的表面,或C 84嵌入Si表面内,并且Si表面单自组装层,如使用UHV-AFM测定。重新打印与基准11的许可。 请点击此处查看该图的放大版本。
图5: 流程图建立仿真模型的戏剧性说明在MD模拟的设置从单个C层84和Si(111)7×7表面嵌入到硅(111)7×7型为C 84单层。详细步骤可以在协议的第5节中看到。 请点击此处查看该图的放大版本。
图6: 缩进力和接触应力分析压痕力(黑色)和联系的Ç84应力(蓝色)为压痕深度的函数。插入显示相应的快照,其中不同的颜色表示的所有原子的相应的冯米塞斯应变(εVM)。要清除显示应变局部化,只能用εVM原子> 0.08快照所示。 请点击此处查看该图的放大版本。
图7:。 硬度和模量减小分析 ( 一 )硬度和(b)模量减小的变化作为压痕深度对Si表面的C 84单层功能,请点击这里查看该图的放大版本。
图8: 加载刚度分析 。载入刚度与用AFM实验对于C 84 / Si的比较由MD模拟确定的距离的函数。从参考16修改。 请点击此处查看该图的放大版本。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
在这项研究中,我们证明Ç84的Si衬底上的自组装单层的通过新颖的退火工艺的制造( 图1)。这个方法也可以用于制备其它种类的纳米颗粒的嵌埋半导体基板。的C 84 -嵌入式Si衬底,其特征在用UHV-STM( 图2),场发射光谱仪,光致发光光谱,MFM和SQUID( 图3)的原子比例。
对应的C- 84纳米力学性能( 即 ,应力)的附着强度-嵌入式Si衬底可以用原子力显微镜( 图4)来测量。我们的结果表明,该Ç84的硬度-嵌入式Si衬底比得上碳化硅和硅的表面,使之适用作为用于切削工具,以及在MEMS器件的薄膜研磨材料。
jove_content“>在模拟部分中,冯米塞斯应变(εVM)分析能够探测原子结构,这是在实验中观察到非常困难的局部变形。但是,这是不可能表征相变。这里,我们提出了一些有用的指数,如配位数和HA索引23以检查相变,在压痕模型的设置,我们要指出的是,在俯视方向上的衬底的尺寸必须至少大三倍比探头用于消除尺寸效应和边界条件的限制,这将影响到动态和原子力流动的直径。另外,由于MD模拟的时间限制,来研究缩进的过程中,探头应以非常快的速度相比,在实验施加试样。我们注意到,这样的加载速度过高出来的长时间在OMIC扩散和迁移行为,但它仍然是合适的观察和根据机械负载24描述的塑性变形行为和材料性能,因为这些结果可以作为大致准静态的性质25被识别。一种替代理论,命名为并行副本动力学(PRD)26,已经发展到基本上加速仿真时间,但它需要大量的计算资源。
从MD模拟研究获得的数据与所述AFM压痕试验( 图8)的协议;另外,将C 84的硬度和模量减小-嵌入式Si衬底比得上其它的Si衬底。这些数据表明,Ç84 -嵌入式Si衬底可以具有在光电显著影响和稀磁性半导体(DMS)纳米器件。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
作者什么都没有透露。
Materials
Name | Company | Catalog Number | Comments |
Silicon wafer | Si(111). Type/Dopant: P/Boron; Resistivity: 0.05-0.1 Ohm·cm | ||
Carbon, C84 | Legend Star | C84 powder, 98% | |
Hydrochloric acid | Sigma-Aldrich | 84422 | RCA, 37% |
Ammonium | Choneye Pure Chemical | RCA, 25% | |
Hydrogen peroxide | Choneye Pure Chemical | RCA, 35% | |
Nitrogen | Ni Ni Air | high-pressure bottle, 95% | |
Tungsten | Nilaco | 461327 | wire, diameter 0.3 mm, tip |
Sodium hydroxide | UCW | 85765 | etching Tungsten wire for tip |
Acetone | Marcon Fine Chemicals | 99920 | suitable for liquid chromatography and UV-spectrophotometry |
Methanol | Marcon Fine Chemicals | 64837 | suitable for liquid chromatography and UV-spectrophotometry |
UHV-SPM | JEOL Ltd | JSPM-4500A | Ultrahigh Vacuum Scanning Tunneling Microscope and Ultrahigh Vacuum Atomic Force Microscope |
Power supply | Keithley | 237 | High-Voltage Source-Measure Unit |
SQUID | Quantum desigh | MPMS-7 | Magnetic field strength: ±7.0 Tesla, Temperature range: 2–400 K, Magnetic-dipole range: 5 × 10-7 – 300 emu |
ALPS | National Center for High-performance Computing, Taiwan | Advanced Large-scale Parallel Supercluster, 177Tflops; 25,600 CPU cores; 73,728 GB RAM; 1,074 TB storage |
References
- Kroto, H. W., Heath, J. R., O'Brien, S. C., Curl, R. F., Smalley, R. E.
C60: Buckminsterfullerene. Nature. 318, 162-163 (1985). - Zhu, Z. P., Gu, Y. D. Structure of carbon caps and formation of fullerenes. Carbon. 34, 173-178 (1996).
- Margadonna, S., et al. Crystal Structure of the Higher Fullerene C84. Chem. Mater. 10, 1742-1744 (1998).
- Diederich, F., et al. The Higher Fullerenes: Isolation and Characterization of C76, C84, C90, C94, and C70O, an Oxide of D5h-C70. Science. 252, 548-551 (1991).
- Lv, Z., Deng, Z., Xu, D., Li, X., Jia, Y. Efficient organic light-emitting diodes with C60 buffer layer. Displays. 30, 23-26 (2009).
- Tokunaga, K. On the difference in electronic properties between fullerene C60 and C60X2. Chem. Phys. Lett. 476, 253-257 (2009).
- Basa, D. K., Smith, F. W. Annealing and crystallization processes in a hydrogenated amorphous Si---C alloy film. Thin Solid Films. 192, 121-133 (1990).
- Neudeck, P. G., Powell, J. A. Performance limiting micropipe defects in silicon carbide wafers. IEEE Electron Device Lett. 15, 63-65 (1994).
- Huang, C. P., Su, C. C., Ho, M. S. Intramolecular Structures of C60 and C84 Molecules on Si(111)-7x7 Surfaces by Scanning Tunneling Microscopy. Appl. Surf. Sci. 254, 7712-7717 (2008).
- Huang, C. P., Su, C. C., Su, W. S., Hsu, C. F., Ho, M. S. Nano-measurements of electronic and mechanical properties of fullerene embedded Si(111) surfaces. Appl. Phys. Lett. 97, 061908 (2010).
- Huang, C. P., Hsu, C. F., Ho, M. S. Investigation of Fullerene Embedded Silicon Surfaces with Scanning Probe Microscopy. J. Nanosci. Nanotechnol. 10, 7145-7148 (2010).
- Huang, C. P., Su, W. S., Su, C. C., Ho, M. S. Characteristics of Si(111) surface with embedded C84 molecules. RSC Adv. 3 (111), 9234-9239 (2013).
- Method of forming self-assembled and uniform fullerene array on surface of substrate. US Patent. Ho, M. S., Huang, C. P. , US9109278 (2015).
- Plimpton, S. J. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comp. Phys. 117, 1-19 (1995).
- Su, W. S., et al. Using a functional C84 monolayer to improve the mechanical properties and alter substrate deformation. RSC Adv. 5, 47498-47505 (2015).
- Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO- the Open Visualization Tool. Modelling Simul. Mater. Sci. Eng. 18, 015012 (2010).
- Watson, P. R., Hove, M. A. V., Hermann, K. Open Surface Structure Database. , Available from: http://www.fhi-berlin.mpg.de/KHsoftware/oSSD (2014).
- Tersoff, J. New empirical model for the structural properties of silicon. Phys. Rev. Lett. 56, 632 (1986).
- Stuart, S. J., Tutein, A. B., Harrison, J. A. A reactive potential for hydrocarbons with intermolecular interactions. Journal of Chemical Physics. 112, 6472 (2000).
- Rapaport, D. C. The Art of Molecular Dynamics Simulations. , Cambridge University Press. (1997).
- Fowler, P. W. Cn Fullerenes. , Available from: http://www.nanotube.msu.edu (2013).
- Chandra, N., Namilae, S., Shet, C. Local elastic properties of carbon nanotubes in the presence of Stone-Wales defects. Phys. Rev. B. 69, 094101 (2004).
- Honeycutt, J. D., Andersen, H. C. Molecular Dynamics Study of Melting and Freezing of Small Lennard-Jones Clusters. J. Phys. Chem. 91, 4950 (1987).
- Schulz, M. J., Kelkar, A. D., Sundaresan, M. J. Nanoengineering of Structural, Functional and Smart Materials. , Taylor & Francis. 736 (2005).
- Ju, S. P., Wang, C. T., Chien, C. H., Huang, J. C., Jian, S. R. The Nanoindentation Responses of Nickel Surfaces with Different Crystal Orientations. Molecular Simulation. 33, 905-917 (2007).
- Mishin, Y., Suzuki, A., Uberuaga, B. P., Voter, A. F. Stick-slip behavior of grain boundaries studied by accelerated molecular dynamics. Phys. Rev. B. 75, (2007).