掺杂硅烯的第一性原理计算(DOC)

题目掺杂硅烯的第一性原理计算

学生姓名张正龙学号1210014048 所在学院物理与电信工程学院

专业班级物理学1202

指导教师潘峰

完成地点陕西理工学院

2016 年6 月5 日

掺杂硅烯的第一性原理计算

张正龙

（陕西理工学院 物理与电信工程学院 物理学专业 2012级（2）班 陕西 汉中 723001）

指导教师：潘峰

[摘 要] 硅烯是一种新型的二维材料，是一个非常有意义的研究体系，其结构与石墨烯的结构类似。由于石墨烯固有的“零带隙”，使其通过一些化学方法打开的带隙难以控制。而硅烯本身不仅具有独特的电子性质，而且能通过吸附、掺杂等原子层面的方法对其电子结构进行设计与改性。制备硅烯和研究硅烯的电子结构已经成为了目前研究的一个重点，而且打开并且硅烯的能量带隙也是研究其电子性质的重点。本文主要介绍了完美硅烯的结构和掺杂原子等原子层面的修饰手段对硅烯电子结构进行设计和修改，并分析对其电子结构的影响。本文主要利用基于第一性原理的密度泛函理论计算方法并借助MS 计算软件的计算和图像显示，研究了掺杂钾原子对打开硅烯能隙的影响以及研究如何打开能隙的新途径。

[关键词]第一性原理；密度泛函理论；硅烯；能隙

引言

硅烯，是一种新型的二维材料，是一种单元子层的二维硅薄膜[1]，其结构与石墨烯的结构类似，目前已经能在实验上成功制备。通过理论计算，我们知道了位于同一个平面的硅原子组成的二维材料不是很稳定，硅原子更加倾向于形成原子不共面的上下翘曲结构。所以与石墨烯不同的是，石墨烯的结构是一种平面的二维蜂窝状结构，而硅烯也是六元环结构，只是六元环中有三个硅原子向上翘曲。硅烯的这种弱翘曲结构让其产生了许多异于石墨烯的优良电子性质。硅烯具有更强的自旋-轨道耦合，可以在它的狄拉克点打开能量带隙，并且比石墨烯由于自旋打开的能量带隙大得多，也对电场有更强的反应，能与外来的原子、分子发生相互的作用。

硅烯具有极高的载流子迁移率，极适合制作成为高效率的场效应管，而打开硅烯的能量带隙是它面向电子器件应用的关键，是以硅烯能量带隙的打开成为了实现应用的关键问题。在本文中，我们通过Material Studio 软件用碱金属原子钾原子对硅烯中的硅原子进行替换，通过分析其能带图和态密度图，进而计算硅烯的能量带隙，研究是否打开硅烯中的能量带隙，并且通过和完美硅烯的能量带隙图的对比，研究打开能量带隙的新方法。

1 计算方法及计算软件 1.1 第一性原理计算方法

量子力学第一性原理计算方法 [2]，仅仅只需要五个物理基本常数来计算：即电子质量、电子电量、普朗克常数、光速、和玻尔兹曼常数，通过对薛定谔方程求解获得材料的电子结构，就预测微观体系的性质和状态而不需要依赖于任何实验经验常数，进而预测材料的组成成分、结构与性能之间的关系。

第一性原理算法对于两个电子以上的多电子体系，很难求出解析的解，需要通过计算机的计算能力来计算，但是半经验法就有优势就是复杂体系效果会更加好一些。实践证明：这两种方法是互补的，不排斥的，可以一同使用，通过第一性原理计算计算机模拟得到的参数带入半经验法去描述物质的性质，进而去发现和解决问题。 我们得到两个电子以上的薛定谔方程为， ˆ(,)(,)H H r R E r R ψψ= (1.1)

可以计算出的哈密顿量为，

ˆ()()(,)N e ee NN Ne

H T T T r V R V r R =++++ (1.2) 222j N R j j T M ∇∑（R ）=-

(1.3)

222i e r i i T m ∇∑（r ）=- (1.4) 2

ee ,1()2i i i i e V r r r ''

=-∑ (1.5) ,1()()2NN N j j j j V R V R R ''=-∑

(1.6) ,(,)()Ne e N j j i j V r R V

r R -=--∑ (1.7)

以上方程表示的是描述固体的非相对量子力学的基础，其中T N 表示原子核的动能，T ee 表示电子的动能，V ee 表示电子与电子之间的库伦相互作用，V NN 表示原子核与原子核之间的相互作用，V Ne 表示原子核与电子之间的相互作用。

1.2 密度泛函理论(DFT)

通过学习，我们知道量子力学一种最流行的表达方式就是波动力学薛定谔的形式，这是波函数及其运动方程的核心-薛定谔方程。对于给定的系统中，我们可以得到全部的信息都包含在该系统的波函数。对一个外势场ʋ(r)中的N 电子体系，量子力学的波动力学范式可以表示成

12()(,......)N r r r r υ⇒ψ⇒可观测量 （1.8） 即对给定的外势，将其带入薛定谔方程计算可以得到电子波函数，进一步通过波函数来计算力学量算符的期望值就可以得到我们所需要的全部可观测量的期望值。电子密度是这些可观测量中的一个，

333*2322()...(,...)(,...)N N N n r N d r d r d r r r r r r r =ψψ⎰⎰⎰ （1.9）

目前有很多方法可以用来解薛定谔方程，比如物理中基于费曼图和格林函数的微扰方法。

当用量子力学处理真实的物理化学体系时，传统的波动力学方法便显的有点不能满足条件，不能满足我们所需要处理的问题。在大多数情况下，我们仅仅只需要与实验相关的信息，如能量、密度等简单容易通过测量得到的信息。所以，我们希望可以通过使用一些简单的物理量来构造并且完善我们的理论。单粒子格林函数G (r,t ;r',t')就是这样的一个物理量。我们可以通过Dyson 方程求解单粒子格林函数，从单粒子格林函数中，我们得到单粒子密度矩阵为，

(,)lim (,;,)t t r r i G r r t t γ'→'''=- （1.10）

或者等价的

333*2322(,)...(,...)(,...)N N N r r N d r d r d r r r r r r r γ''=ψψ⎰⎰⎰ （1.11）

通过密度矩阵，我们可以求得所有单体算符的期望值。因为从格林函数到密度矩阵的计算过程中丢到了一些信息，我们不能简单直接的从密度矩阵求出哈密顿量的期望值，我们需要先获取实验上直接可观测的粒子密度,

n()(,)|r r r r r γ'→'=

（1.12）