平面波基组下GW和Bethe-Salpeter方程的三次标度算法实现及应用

来源 :中国科学技术大学 | 被引量 : 0次 | 上传用户:nihaohaoya
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
基于密度泛函理论的第一性计算方法已经被大量应用于预测和研究材料的各种物理和化学性质,已经成为独立于实验和理论的第三种研究方法。密度泛函理论的成功在于创造性利用密度作为能量的泛函,将3N维的Schr(?)dinger方程转为了三维问题,大大降低了计算量。但是密度泛函理论是基于变分原理的基态理论,这使得其在应用到激发态过程中存在很多问题,比如无法准确地预测分子和固体体系的光谱和激发态能级等。同时由于到目前为止没有精确的交换关联泛函形式,这使得密度泛函理论在预测半导体材料的带隙非常依赖于交换关联泛函的选取,特别是LDA和GGA都会严重低估能隙。而基于波函数的Hartree-Fock方法由于其也是基于变分原理来得到基态能量,且只考虑了单电子近似,使得其无法准确地描述激发态性质和大大高估了半导体体系的带隙。基于格林函数的多体微扰理论在一定程度上解决上述困难,特别是基于单体格林函数的的GW近似可以很好地描述半导体体系的准粒子性质,得到和实验一致的带隙,基于两体格林函数的Bethe-Salpeter方程(BSE)可以描述激子的运动。本论文主要研究周期性体系,因此从平面波基组下的密度泛函理论出发,主要研究GW和BSE的三次标度算法实现及其在二维材料中的简单应用,大致分为如下六个章节。首先第一章主要回顾了量子多体系统和密度泛函理论的发展历史。简单介绍了量子多体系统和求解量子多体系统的困难,进一步介绍了为了简化计算而将核和电子分开处理的Born-Oppenheimer(BO)近似。在BO近似下发展的Hartree-Fock近似,在一定程度可以减少计算量,但是其方程仍然为3N维的。通过利用密度作为能量的泛函可以将多体系统Schr(?)dinger方程从3N维问题转化为3维问题,大大减少了计算量。但是密度泛函理论的最大问题是寻找精确的交换关联泛函,目前已经有一些常用的交换关联泛函。在本章的最后介绍一些常用的密度泛函软件。第二章主要是介绍基于格林函数的多体微扰理论(MBPT)。在本章首先介绍了激发态电子结构和激子分类,然后回顾了 GW和BSE的发展历史。通过利用Feynman传播子引入格林函数和准粒子的概念,并给出了单体和两体格林函数的具体形式。在GW介绍中,重点介绍了格林函数的运动方程、Dyson方程和Hedin方程。基于格林函数的运动方程,通过引入自能项可以得到准粒子方程,至此就可以自洽求解。由于自洽求解过程非常耗时而且效果也不是特别好,因此发展了 G0W0方法。紧接着介绍了三种目前常见的G0W0近似的方法。接着还介绍了基于两体格林函数的Bethe-Salpeter方程(BSE)及其与含时密度泛函(TDDFT)的区别。最后介绍了一些常用的GW和BSE软件。第三章内容主要是利用低秩分解算法和Cauchy积分在平面波基组下实现三次标度的G0W0计算。G0W0近似在实际应用中是非常耗时且占用内存,通常比相同体系下DFT计算大一到两个数量级。通常情况下G0W0计算复杂度为O(Ne4~5),这主要是因为G0W0计算中包含多个两电子四中心积分,因此如何降低两电子四中心问题是本章研究的重点。首先介绍了插值可分离密度拟合(ISDF)算法,而实现ISDF需要两步,首先是利用QRCP或者K-means获得插值点,然后通过最小二乘法得到插值矢量。通过利用ISDF方法对Khatri-Rao积进行低秩分解,再利用交换积分顺序的方法将自能项的计算降低为三次标度。但是由于介电函数部分含有占据态和空态的耦合项,这里可以利用Cauchy积分来解耦,实现了三次标度的静态的COHSEX的计算。结果部分我们分析了 ISDF低秩分解算法的精度和计算效率,可以在不损失精度的情况下取得接近10倍的加速效果。同时还发现做了低秩分解以后介电函数求逆变成了主要的耗时项,通过利用Sherman-Morrison-Woodburg(SMW)公式将介电函数求逆转化为多个小矩阵的乘积的形式,进一步降低介电矩阵求逆的内存和时间消耗。本章最后还介绍了高精度的G0W0近似GPP在低秩分解算法下的精度。第四章主要是介绍了利用低秩分解算法和迭代对角化算法在平面波基组下三次标度的BSE计算的实现方法。本节首先介绍了构造BSE的流程,分析了BSE计算中瓶颈部分,可以发现BSE计算由两个最主要的过程构成,分别为构造Bethe-Salpeter Hamiltonian(BSH)和对角化BSH得到本征值和本征矢。首先利用ISDF低秩分解算法降低了构造BSH中两电子四中心积分计算量大的问题。由于构成BSH中包含特殊项-屏蔽的交换项,这使得其无法直接利用迭代对角算法,这里我们创造性的提出通过对子空间的波矢做重排,再利用隐式迭代对角化实现了三次标度的BSE计算。结果部分我们分析了利用ISDF和迭代算法以后的计算精度和计算效率,发现利用ISDF低秩分解算法和隐式迭代相结合的方法,可以取得接近25倍的加速效果。第五章给了一个GW+BSE在二维铁磁半导体中的简单应用。通过分析不同石墨烯量子点的磁性性质,选择了本身就是铁磁的三角形之字边界的石墨烯量子点作为基础单元,通过碳原子间的四元环链接构造了周期性结构的二维六方形周期性格点。通过第一性原理计算确认其为铁磁性的半导体,而且还是特殊的双极磁性半导体。然后利用GW+BSE计算了其准粒子能带结构和激发态性质,分析了光谱和激子波函数。本章最后利用声子谱和AIMD,分析了其稳定性,并通过经验公式和基于Heisenberg模型的蒙特卡洛计算得到磁转变的居里温度,确定其为室温铁磁半导体。第六章对毕业论文的主要内容做了总结,并展望了激发态电子结构未来的发展趋势。在论文附录部分重点介绍二维材料的力学性质。正文部分主要是激发态电子结构的研究,而材料的基态电子结构同样重要。通过第一性原理计算,单层β-Te具有高的载流子迁移率、合适的带隙和稳定性。利用应力可以改变β-Te实现电子迁移率从沿着扶手椅方向最大转换到之字型方向最大,从而实现β-Te电子电导方向的反转。为了证实了其可以能够承受的应力,进一步研究β-Te的力学性质和稳定性,这为β-Te在柔性电子学中的应用提供了理论支撑。前人研究发现单层黑磷或者类黑磷材料具有面外负泊松比,我们通过计算发现类黑磷材料同样具有面内负泊松比,这证实了类黑磷材料是双拉胀材料。
其他文献
热塑性材料具有重量轻、导热系数低、易加工、价格低廉等优点,广泛应用于建筑外墙保温系统、管道、家具、玩具、包装、医疗健康及高科技等领域中。但易燃、热值高、受热熔融软化的特点使其具有较大的火灾危险性,近年来与之相关的火灾事故屡见报道,造成了重大的人员伤亡和财产损失。本论文以PMMA(聚甲基丙烯酸甲酯,Poly(methyl methacrylate))材料为主要研究对象,应用小尺寸实验与理论分析相结合
近年来,学龄前儿童数量不断增加,研究学龄前儿童群体的疏散运动特性,对提高学龄前儿童疏散效率与安全水平具有重要意义。瓶颈作为行人通行设施中的典型结构,因其对行人流的约束已获得了广泛关注。然而,当前的瓶颈疏散研究忽略了学龄前儿童与成年人的疏散特性差异,研究学龄前儿童的瓶颈疏散特性能够为儿童活动场所的安全及应急预案设计提供支撑。本文通过开展可控实验,研究了学龄前儿童瓶颈疏散的疏散运动特性,量化了瓶颈宽度
<正>腹痛是常见病多发病,疼痛的来源可以是腹腔脏器,也可以是腹部壁层神经引起,目前内外科医师更注重于前者,而极易将最浅表的腹壁疼痛误为腹腔脏器病变所致。现将作者近三年来54例脊柱相关性腹痛的诊治体会报告如下。
会议
上转换纳米粒子(upconversionnanoparticles,UCNPs)能够将低能量近红外(near infrared,NIR)激发光转换成高能量可见或紫外光。近红外激发的UCNPs具有近红外激发、光稳定性优异、信噪比高和多色发射等优点,在生物医学、环境检测、食品安全等领域具有广泛应用。然而,传统上转换发光探针具有发光共振能量转移(luminescence resonance energy
目前受控热核聚变能被认为是能够有效解决人类未来能源需求的清洁新能源,燃烧等离子体、托卡马克聚变堆工程技术、聚变堆各类材料是影响托卡马克聚变堆走向产业化应用的三个主要方面。在磁约束热核聚变装置中,面向等离子体部件的表面材料会与等离子体发生作用,形成杂质进入等离子体,杂质辐射对高参数等离子体稳定运行产生重要影响。前期大多数中大型托卡马克实验装置均使用碳瓦作为第一壁和偏滤器靶板表面材料,在未来托卡马克聚
随着海上运输业的蓬勃发展,船舶火灾安全问题日渐凸显。根据事故统计船舶火灾多发于机舱,而七成以上的船舶机舱火灾事故是发生在海上航行过程中。海上航行中的船舶不可避免会受到海风的作用,事故分析表明环境风影响下的船舶火灾危害性更大。前人研究中通常将船舶舱室简化为顶部开口舱室,但鲜少考虑环境风这一外界因素的影响。因此,研究环境风对顶部开口舱室火灾的影响,能够为此类船舶火灾事故的预防和应急救援对策的制定提供理
黄热病毒家族的寨卡病毒(ZikaVirus,ZIKV)由于导致人类神经系统疾病如小头畸形和Guillain-Barré综合征等而受到广泛关注。寨卡病毒是一种单链正义RNA病毒,其RNA基因组编码的非结构蛋白3(NS3)包含一个蛋白酶结构域和一个解旋酶结构域,N端丝氨酸蛋白酶结构域需要与NS2B蛋白因子相互作用才能发挥活性;C端解旋酶结构域在核糖核苷三磷酸(Nucleoside triphospha
单视角深度估计任务是计算机视觉领域中极具挑战性的问题之一,其目标是给定单张彩色图像,预测图像上每个像素点的深度值。由于三维空间到图像平面的透视投影,该任务存在天然的歧义性,是一个非适定问题。近年来,随着深度学习技术的发展,单视角深度估计的算法性能得到大幅提升,但仍然面临三个困难:(1)通过深度卷积神经网络预测的深度图缺乏对场景精细结构的描述;(2)深度标注数据获取难度大;(3)单视角深度估计网络的
近年来,太赫兹(THz)电磁波在医疗生物检测、材料分析、无线通信、大气监测等广泛领域展现出越来越重要的应用前景。THz辐射源是支撑THz技术应用和发展的关键器件。然而,由于THz波处在电子学和光子学的中间地带,基于传统电子学方法和光子学方法的THz源在输出功率、频率、便携性等方面难以满足实用需求。自由电子驱动的电磁辐射源在THz频段可以实现大功率输出,是一类重要的THz源。但是这类装置通常体积庞大
服役于人体生理环境中的植入金属材料往往由于腐蚀、磨损等作用而产生金属离子释放。某些适量释放的金属离子虽然在促进组织再生等方面发挥积极作用,但不可否认,过量的离子积累将对人体产生毒副作用,进而直接影响材料的生物相容性。通常,植入金属离子释放行为与材料的耐蚀能力息息相关,而作为评价植入金属材料耐蚀性的通用、实用手段,腐蚀电化学测试技术已广泛应用。因此,探究模拟人体环境中金属离子释放的动力学过程,阐述离