【免费下载】第三节第一性原理计算简介

合集下载

第一性原理计算方法讲义

第一性原理计算方法讲义

第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现,使得人们对于物质微观结构的认识日益深入。

原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。

目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。

1964年,Hohenberg和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT)。

它建立在非均匀作为基本变量。

1965年,Kohn和Sham 电子气理论基础之上,以粒子数密度()r提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。

第一性原理简介

第一性原理简介

第⼀性原理简介1什么是第⼀性原理?根据原⼦核和电⼦互相作⽤的原理及其基本运动规律,运⽤量⼦⼒学原理,从具体要求出发,经过⼀些近似处理后直接求解薛定谔⽅程的算法,称为第⼀性原理。

⼴义的第⼀原理包括两⼤类,以Hartree-Fock⾃洽场计算为基础的从头算和密度泛函理论(DFT)计算。

从定义可以看出第⼀性原理涉及到量⼦⼒学、薛定谔⽅程、Hartree-Fock⾃洽场、密度泛函理论等许多对我来说很陌⽣的物理化学定义。

因此我通过向师兄请教和上⽹查资料⼀点点的了解并学习这些知识。

2第⼀性原理的作⽤以密度泛函理论(DFT)为基础以及在此基础上发展起来的简单⽽具有⼀定精度的局域密度近似(LDA)和⼴义梯度近似(GGA)的第⼀性原理电⼦结构计算⽅法,与传统的解析⽅法⼀样,不但能够给出描述体系微观电⼦特性的物理量如波函数、态密度、费⽶⾯、电⼦间互作⽤势等,以及在此基础上所得到的体现体系宏观物理特性的参量如结合能、电离能、⽐热、电导、光电⼦谱、穆斯堡尔谱等等,⽽且它还可以帮助⼈们预⾔许多新的物理现象和物理规律。

密度泛函计算的⼀些结果能够与实验直接进⾏⽐较,⼀些应⽤程序的发展乃⾄商业软件的发布,导致了基于密度泛函理论的第⼀原理计算⽅法的⼴泛应⽤。

密度泛函理论(DFT)为第⼀性原理中的⼀类,在物理系、化学、材料科学以及其他⼯程领域中,密度泛函理论(DFT)及其计算已经快速发展成为材料建模模拟的⼀种“标准⼯具”。

密度泛函理论可以计算预测固体的晶体结构、晶格参数、能带结构、态密度(DOS)、光学性能、磁性能以及原⼦集合的总能等等。

3第⼀性原理怎么⽤?⽬前我所学到的利⽤第⼀性原理的软件为Material Studio、V ASP软件。

其中Materials Studio(简称MS)是专门为材料科学领域研究者开发的⼀款可运⾏在PC上的模拟软件。

使化学及材料科学的研究者们能更⽅便地建⽴三维结构模型,并对各种晶体、⽆定型以及⾼分⼦材料的性质及相关过程进⾏深⼊的研究。

第一性原理计算方法讲义

第一性原理计算方法讲义

第一性原理计算方法讲义标准化管理部编码-[99968T-6889628-J68568-1689N]第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现,使得人们对于物质微观结构的认识日益深入。

原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。

目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。

1964年,Hohenberg和Kohn提出了严格的密度泛函理论(DensityFunctional Theory, DFT)。

第一性原理计算流程与原理公式

第一性原理计算流程与原理公式

第一性原理计算流程与原理公式下载温馨提示:该文档是我店铺精心编制而成,希望大家下载以后,能够帮助大家解决实际的问题。

文档下载后可定制随意修改,请根据实际需要进行相应的调整和使用,谢谢!并且,本店铺为大家提供各种各样类型的实用资料,如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等,如想了解不同资料格式和写法,敬请关注!Download tips: This document is carefully compiled by theeditor. I hope that after you download them,they can help yousolve practical problems. The document can be customized andmodified after downloading,please adjust and use it according toactual needs, thank you!In addition, our shop provides you with various types ofpractical materials,such as educational essays, diaryappreciation,sentence excerpts,ancient poems,classic articles,topic composition,work summary,word parsing,copy excerpts,other materials and so on,want to know different data formats andwriting methods,please pay attention!一、流程1. 确定研究体系:明确要研究的物质或材料的结构和组成。

第一性原理计算简述

第一性原理计算简述

第一性原理计算简述第一性原理,英文First Principle,是一个计算物理或计算化学专业名词,广义的第一性原理计算指的是一切基于量子力学原理的计算。

我们知道物质由分子组成,分子由原子组成,原子由原子核和电子组成。

量子力学计算就是根据原子核和电子的相互作用原理去计算分子结构和分子能量(或离子),然后就能计算物质的各种性质。

从头算(ab initio)是狭义的第一性原理计算,它是指不使用经验参数,只用电子质量,光速,质子中子质量等少数实验数据去做量子计算。

但是这个计算很慢,所以就加入一些经验参数,可以大大加快计算速度,当然也会不可避免的牺牲计算结果精度。

那为什么使用“第一性原理”这个字眼呢?据说这是来源于“第一推动力”这个宗教词汇。

第一推动力是牛顿创立的,因为牛顿第一定律说明了物质在不受外力的作用下保持静止或匀速直线运动。

如果宇宙诞生之初万事万物应该是静止的,后来却都在运动,是怎么动起来的呢?牛顿相信这是由于上帝推了一把,并且牛顿晚年致力于神学研究。

现代科学认为宇宙起源于大爆炸,那么大爆炸也是有原因的吧。

所有这些说不清的东西,都归结为宇宙“第一推动力”问题。

科学不相信上帝,我们不清楚“第一推动力”问题只是因为我们科学知识不完善。

第一推动一定由某种原理决定。

这个可以成为“第一原理”。

爱因斯坦晚年致力与“大统一场理论”研究,也是希望找到统概一切物理定律的“第一原理”,可惜,这是当时科学水平所不能及的。

现在也远没有答案。

但是为什么称量子力学计算为第一性原理计算?大概是因为这种计算能够从根本上计算出来分子结构和物质的性质,这样的理论很接近于反映宇宙本质的原理,就称为第一原理了。

广义的第一原理包括两大类,以Hartree-Fork自洽场计算为基础的ab initio从头算,和密度泛函理论(DFT)计算。

也有人主张,ab initio专指从头算,而第一性原理和所谓量子化学计算特指密度泛函理论计算。

根据原子核和电子互相作用的原理及其基本运动规律,运用量子力学原理,从具体要求出发,经过一些近似处理后直接求解薛定谔方程的算法,习惯上称为第一原理。

第一性原理计算 ppt课件

第一性原理计算  ppt课件
j 1
ppt课件
9
• 除非i等于j,否则Lagrange乘法因子为0 .
Fi i ij j
•每一个电子都被假设在包括原子核和其他电子的固定区域中运动 •每一个电子方程得到的解都会影响系统中其他电子的解 自洽场假设
解的思路 •得到一个H-F本征方程的试探解,用来计算库仑和交还能。 •解H-F方,给出第二套解。 •一次类推。 •SCF方法逐渐得到对应越来越低能量的单电子的解。直到一点,在这一 点所有的电子的方程不再改变
ij Sij 0
ij
ppt课件
6
• Lagrange乘法因子被写为-ij
E
ij Sij 0
ij
• 多电子系统和单电子系统的最大区别在于存在电子之间的
作用,包括库仑和交换积分
• 找到多电子系统最好的波函数,保留轨道图像
• 找到一个解,当一个电子的轨道变化时由于电子耦合作用 会影响其他自旋轨道上的电子的运动,在上述情况下,这
ppt课件
12
• 完整的轨道,我们要用角度部分乘以径向函数部分
1s (r) 3 / exp(r)
2s (r) 5 / 3 rexp(r)
2 pz (r) 5 / exp(r)cos
ppt课件
13
• Slater给出了一系列的轨道指数的确定原则
Z
n*
Z为原子数,为屏蔽常数,n*为有效的主量子数
n* :
n=1,2,3时取与n相同的值;n=4,5,6s时,取3.7,4.0,4.2
屏蔽常数按下列规则确定 :
把轨道划分如下组: (1s);(2s,2p);(3s,3p);(3d);(4s,4p);(4f);(5s,5p);(5d)

第一性原理计算原理和方法

第一性原理计算原理和方法

第二章 计算方法及其基本原理介绍化学反应的本质是旧键的断裂和新建的形成,参与成键原子的电子壳层重新组合是导致生成稳定多原子化学键的明显特征。

因此阐述化学键的理论应当描写电子壳层的相互作用与重排,借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学,为解决这一问题提供了一般的方法,然而,对于一些实际的体系,不引入一些近似,就不可能求解其Schrodinger 方程。

这些近似使一般量子力学方程简化为现代电子计算机可以求解的方程。

这些近似和关于分子波函数的方程形成计算量子化学的数学基础。

2.1 SCF-MO 方法的基本原理分子轨道的自洽场计算方法(SCF-MO)是各种计算方法的理论基础和核心部分,因此在介绍本文计算工作所用方法之前,有必要对其关键的部分作一简要阐述。

2.1.1 Schrodinger 方程及一些基本近似 为了后面介绍各种具体在自洽场分子轨道(SCF MO)方法方便,这里将主要阐明用于本文量子化学计算的一些重要的基本近似,给出SCF MO 方法的一些基本方程,并对这些方程作简略说明,因为在大量的文献和教材中对这些方程已有系统的推导和阐述[1-5]。

确定任何一个分子的可能稳定状态的电子结构和性质,在非相对论近似下,须求解定态Schrodinger 方程''12121212122ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =⎥⎥⎦⎤⎢⎢⎣⎡-++∇-∇-∑∑∑∑∑∑≠≠ (2.1) 其中分子波函数依赖于电子和原子核的坐标,Hamilton 算符包含了电子p 的动能和电子p 与q 的静电排斥算符,R AB =R 图2-1分子体系的坐标∑∑≠+∇-=p q p pqp e r H 12121ˆ2 (2.2) 以及原子核的动能∑∇-=A A AN M H 2121ˆ (2.3) 和电子与核的相互作用及核排斥能∑∑≠+-=p A B A AB B A pAA eN R Z Z r Z H ,21ˆ (2.4) 式中Z A 和M A 是原子核A 的电荷和质量,r pq =|r p -r q |,r pA =|r p -R A |和R AB =|R A -R B |分别是电子p 和q 、核A 和电子p 及核A 和B 间的距离(均以原子单位表示之)。

第一性原理计算

第一性原理计算

5.1 原子的H-F计算以及Slater规则
• 如果假设电子的分布实球形对称的 • 近似分析函数
Rn(lr)Ylm (,)
•氢原子得到的径向函数不能直接用于多电子原子。 •因为内壳层电子对原子核电荷具有屏蔽作用。 •如果考虑屏蔽效应把轨道指数作适当地修正,仍可以采用氢原子的波函 数形式。
• Slater
第一性原理计算
• Hiicore:电子在裸露原子核作用场中运动的能量。
• 对于在轨道 i运动的电子:如果一个轨道具有两个电子, 则对N/2个电子而言单个电子能量为2 Hiicore.

N /2
2
H
c ii
ore
i 1
• 电子与电子的作用项
• i和j中存在4个电子。
• 在一个轨道上的2个电子以库仑作用的方式与另一个轨道 上的两个电子发生作用。记为4Jij。
•需要猜想密度矩阵P。最简单的方法是采用空矩阵 •H-F计算的结果是一系列K原子轨道,K是计算的基函数的数量 •N个电子被填充到这些轨道上,从最低能量的轨道开始,一个轨道2个电子进行填充。
H-F方程给出了一系列的轨道轨道能量,i,
N/2
i Hiciore (2JijKij)
j1
整个基态的电子能量
N/2
d1(1)2 [Jj(1)Kj(1)]1)(v(1)
j1
库仑算符Jj(1)由自旋轨道j的交互作用引起
Jj(1) d2j(2)r112j(2)
利用轨道函数的基函数线性组合的方法
k
Jj(1)d2 cj
1
1k (2)r12 1cj
(2)
交换项可以写为
K j(1 )i(1 ) [d2 k 1 c j (2 )r 1 12j(2 )k ] 1 c j (2 )

第一性原理计算方法讲义

第一性原理计算方法讲义

第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现,使得人们对于物质微观结构的认识日益深入。

原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。

目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。

1964年,Hohenberg和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT)。

它建立在非均匀作为基本变量。

1965年,Kohn和Sham 电子气理论基础之上,以粒子数密度()r提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。

第一性原理计算方法讲义

第一性原理计算方法讲义

第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现,使得人们对于物质微观结构的认识日益深入。

原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。

目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。

1964年,Hohenberg和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT)。

它建立在非均匀作为基本变量。

1965年,Kohn和Sham 电子气理论基础之上,以粒子数密度()r提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。

第一性原理计算方法讲义

第一性原理计算方法讲义

第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现,使得人们对于物质微观结构的认识日益深入。

原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。

目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。

1964年,Hohenberg和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT)。

它建立在非均匀作为基本变量。

1965年,Kohn和Sham 电子气理论基础之上,以粒子数密度()r提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。

第一性原理计算简述

第一性原理计算简述

第一性原理‎计算简述第一性原理‎,英文Fir‎s t Princ‎iple,是一个计算‎物理或计算‎化学专业名‎词,广义的第一‎性原理计算‎指的是一切‎基于量子力‎学原理的计‎算。

我们知道物‎质由分子组‎成,分子由原子‎组成,原子由原子‎核和电子组‎成。

量子力学计‎算就是根据‎原子核和电‎子的相互作‎用原理去计‎算分子结构‎和分子能量‎(或离子),然后就能计‎算物质的各‎种性质。

从头算(ab initi‎o)是狭义的第‎一性原理计‎算,它是指不使‎用经验参数‎,只用电子质‎量,光速,质子中子质‎量等少数实‎验数据去做‎量子计算。

但是这个计‎算很慢,所以就加入‎一些经验参‎数,可以大大加‎快计算速度‎,当然也会不‎可避免的牺‎牲计算结果‎精度。

那为什么使‎用“第一性原理‎”这个字眼呢‎?据说这是来‎源于“第一推动力‎”这个宗教词‎汇。

第一推动力‎是牛顿创立‎的,因为牛顿第‎一定律说明‎了物质在不‎受外力的作‎用下保持静‎止或匀速直‎线运动。

如果宇宙诞‎生之初万事‎万物应该是‎静止的,后来却都在‎运动,是怎么动起‎来的呢?牛顿相信这‎是由于上帝‎推了一把,并且牛顿晚‎年致力于神‎学研究。

现代科学认‎为宇宙起源‎于大爆炸,那么大爆炸‎也是有原因‎的吧。

所有这些说‎不清的东西‎,都归结为宇‎宙“第一推动力‎”问题。

科学不相信‎上帝,我们不清楚‎“第一推动力‎”问题只是因‎为我们科学‎知识不完善‎。

第一推动一‎定由某种原‎理决定。

这个可以成‎为“第一原理”。

爱因斯坦晚‎年致力与“大统一场理‎论”研究,也是希望找‎到统概一切‎物理定律的‎“第一原理”,可惜,这是当时科‎学水平所不‎能及的。

现在也远没‎有答案。

但是为什么‎称量子力学‎计算为第一‎性原理计算‎?大概是因为‎这种计算能‎够从根本上‎计算出来分‎子结构和物‎质的性质,这样的理论‎很接近于反‎映宇宙本质‎的原理,就称为第一‎原理了。

广义的第一‎原理包括两‎大类,以Hart‎r ee-Fork自‎洽场计算为‎基础的ab‎initi‎o从头算,和密度泛函‎理论(DFT)计算。

第一性原理计算方法讲义

第一性原理计算方法讲义

第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现,使得人们对于物质微观结构的认识日益深入。

原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。

目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。

1964年,Hohenberg和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT)。

它建立在非均匀作为基本变量。

1965年,Kohn和Sham 电子气理论基础之上,以粒子数密度()r提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。

第一性原理简介

第一性原理简介

1什么是第一性原理?根据原子核和电子互相作用的原理及其基本运动规律,运用量子力学原理,从具体要求出发,经过一些近似处理后直接求解薛定谔方程的算法,称为第一性原理。

广义的第一原理包括两大类,以Hartree-Fock自洽场计算为基础的从头算和密度泛函理论(DFT)计算。

从定义可以看出第一性原理涉及到量子力学、薛定谔方程、Hartree-Fock自洽场、密度泛函理论等许多对我来说很陌生的物理化学定义。

因此我通过向师兄请教和上网查资料一点点的了解并学习这些知识。

2第一性原理的作用以密度泛函理论(DFT)为基础以及在此基础上发展起来的简单而具有一定精度的局域密度近似(LDA)和广义梯度近似(GGA)的第一性原理电子结构计算方法,与传统的解析方法一样,不但能够给出描述体系微观电子特性的物理量如波函数、态密度、费米面、电子间互作用势等,以及在此基础上所得到的体现体系宏观物理特性的参量如结合能、电离能、比热、电导、光电子谱、穆斯堡尔谱等等,而且它还可以帮助人们预言许多新的物理现象和物理规律。

密度泛函计算的一些结果能够与实验直接进行比较,一些应用程序的发展乃至商业软件的发布,导致了基于密度泛函理论的第一原理计算方法的广泛应用。

密度泛函理论(DFT)为第一性原理中的一类,在物理系、化学、材料科学以及其他工程领域中,密度泛函理论(DFT)及其计算已经快速发展成为材料建模模拟的一种“标准工具”。

密度泛函理论可以计算预测固体的晶体结构、晶格参数、能带结构、态密度(DOS)、光学性能、磁性能以及原子集合的总能等等。

3第一性原理怎么用?目前我所学到的利用第一性原理的软件为Material Studio、V ASP软件。

其中Materials Studio(简称MS)是专门为材料科学领域研究者开发的一款可运行在PC上的模拟软件。

使化学及材料科学的研究者们能更方便地建立三维结构模型,并对各种晶体、无定型以及高分子材料的性质及相关过程进行深入的研究。

第一性原理计算

第一性原理计算

第一性原理根据原子核和电子互相作用的原理及其基本运动规律,运用量子力学原理,从具体要求出发,经过一些近似处理后直接求解薛定谔方程的算法,习惯上称为第一原理第一性原理通常是跟计算联系在一起的,是指在进行计算的时候除了告诉程序你所使用的原子和他们的位置外,没有其他的实验的,经验的或者半经验的参量,且具有很好的移植性。

作为评价事物的依据,第一性原理和经验参数是两个极端。

第一性原理是某些硬性规定或推演得出的结论,而经验参数则是通过大量实例得出的规律性的数据,这些数据可以来自第一性原理(称为理论统计数据),也可以来自实验(称为实验统计数据)。

但是就某个特定的问题,第一性原理和经验参数没有明显的界限,必须特别界定。

如果某些原理或数据来源于第一性原理,但推演过程中加入了一些假设(这些假设当然是很有说服力的),那么这些原理或数据就称为“半经验的”。

第一性原理,英文First Principle,是一个计算物理或计算化学专业名词,广义的第一性原理计算指的是一切基于量子力学原理的计算。

我们知道物质由分子组成,分子由原子组成,原子由原子核和电子组成。

量子力学计算就是根据原子核和电子的相互作用原理去计算分子结构和分子能量(或离子),然后就能计算物质的各种性质。

从头算(ab initio)是狭义的第一性原理计算,它是指不使用经验参数,只用电子质量,光速,质子中子质量等少数实验数据去做量子计算。

但是这个计算很慢,所以就加入一些经验参数,可以大大加快计算速度,当然也会不可避免的牺牲计算结果精度。

那为什么使用“第一性原理”这个字眼呢?据说这是来源于“第一推动力”这个宗教词汇。

第一推动力是牛顿创立的,因为牛顿第一定律说明了物质在不受外力的作用下保持静止或匀速直线运动。

如果宇宙诞生之初万事万物应该是静止的,后来却都在运动,是怎么动起来的呢?牛顿相信这是由于上帝推了一把,并且牛顿晚年致力于神学研究。

现代科学认为宇宙起源于大爆炸,那么大爆炸也是有原因的吧。

第一性原理计算简述

第一性原理计算简述

第一性原理计算简述第一性原理计算的基本思想是将材料系统中的电子行为完全用量子力学方法描述,并且只基于一些常见的物理规律进行计算,而不依赖于实验数据或经验性参数。

这种方法被认为是计算物理学中最精确的方法之一,可以提供高度精确的材料性质和行为预测。

第一性原理计算的核心是薛定谔方程的求解。

薛定谔方程描述了一个系统的波函数随时间的演化,通过求解该方程可以得到系统的能量、波函数及其他的物理量。

然而,由于薛定谔方程的复杂性,直接求解它对于复杂的体系是不现实的。

因此,第一性原理计算采用了一些近似方法对薛定谔方程进行处理。

在第一性原理计算中,常用的近似方法包括密度泛函理论(DFT)和平面波基组方法。

密度泛函理论是一种计算材料中电子的方法,它通过将波函数的描述换成了电子密度的描述,从而大大简化了计算。

平面波基组方法是一种将波函数展开成平面波的形式,并与周期性边界条件相适应的方法,用于对材料中的电子进行离散化处理。

除了薛定谔方程的求解方法,第一性原理计算还需要一些模型和算法来处理实际系统中的一些问题。

例如,需要考虑电子之间的相互作用,常用的方法有赝势(pseudopotential)和Hartree-Fock方法。

赝势方法将复杂的电子-电子相互作用简化为一个有效的势能,从而加快了计算速度。

Hartree-Fock方法是一种处理多电子系统中电子之间相互作用的方法,它将多体态用单体态的乘积形式进行描述,并采用自洽迭代的方式求解能量。

第一性原理计算可以用于多种材料的性质和行为的预测和解释。

例如,可以通过计算系统的能带结构来预测材料的导电性质;可以通过计算材料的弹性常数来预测其力学性质;可以通过计算材料的反应势垒来预测化学反应的速率等。

此外,第一性原理计算还可以用于设计新的材料,例如预测新的材料的晶体结构和电子性质,从而为材料科学的研究和应用提供宝贵的理论指导。

然而,第一性原理计算也存在一些限制和挑战。

首先,计算方法的复杂性限制了其应用范围。

第一性原理简介

第一性原理简介

2第一性原理的作用以密度泛函理论(DFT)为基础以及在此基础上发展起来的简单而具有一定精度的局域密度近似(LDA)和广义梯度近似(GGA)的第一性原理电子结构计算方法,与传统的解析方法一样,不但能够给出描述体系微观电子特性的物理量如波函数、态密度、费米面、电子间互作用势等,以及在此基础上所得到的体现体系宏观物理特性的参量如结合能、电离能、比热、电导、光电子谱、穆斯堡尔谱等等,而且它还可以帮助人们预言许多新的物理现象和物理规律。

密度泛函计算的一些结果能够与实验直接进行比较,一些应用程序的发展乃至商业软件的发布,导致了基于密度泛函理论的第一原理计算方法的广泛应用。

密度泛函理论(DFT)为第一性原理中的一类,在物理系、化学、材料科学以及其他工程领域中,密度泛函理论(DFT)及其计算已经快速发展成为材料建模模拟的一种“标准工具”。

密度泛函理论可以计算预测固体的晶体结构、晶格参数、能带结构、态密度(DOS)、光学性能、磁性能以及原子集合的总能等等。

3第一性原理怎么用?目前我所学到的利用第一性原理的软件为Material Studio、V ASP软件。

其中Materials Studio(简称MS)是专门为材料科学领域研究者开发的一款可运行在PC上的模拟软件。

使化学及材料科学的研究者们能更方便地建立三维结构模型,并对各种晶体、无定型以及高分子材料的性质及相关过程进行深入的研究。

模拟的内容包括了催化剂、聚合物、固体及表面、晶体与衍射、化学反应等材料和化学研究领域的主要课题。

模块简介Materials Studio采用了大家非常熟悉的Microsoft标准用户界面,允许用户通过各种控制面板直接对计算参数和计算结果进行设置和分析。

目前,Materials Studio软件包括如下功能模块:Materials Visualizer:提供了搭建分子、晶体及高分子材料结构模型所需要的所有工具,可以操作、观察及分析结构模型,处理图表、表格或文本等形式的数据,并提供软件的基本环境和分析工具以及支持Materials Studio的其他产品。

第一性原理简介

第一性原理简介

1什么是第一性原理?根据原子核和电子互相作用的原理及其基本运动规律,运用量子力学原理,从具体要求出发,经过一些近似处理后直接求解薛定谔方程的算法,称为第一性原理。

广义的第一原理包括两大类,以Hartree-Fock自洽场计算为基础的从头算和密度泛函理论(DFT)计算。

从定义可以看出第一性原理涉及到量子力学、薛定谔方程、Hartree-Fock自洽场、密度泛函理论等许多对我来说很陌生的物理化学定义。

因此我通过向师兄请教和上网查资料一点点的了解并学习这些知识。

2第一性原理的作用以密度泛函理论(DFT)为基础以及在此基础上发展起来的简单而具有一定精度的局域密度近似(LDA)和广义梯度近似(GGA)的第一性原理电子结构计算方法,与传统的解析方法一样,不但能够给出描述体系微观电子特性的物理量如波函数、态密度、费米面、电子间互作用势等,以及在此基础上所得到的体现体系宏观物理特性的参量如结合能、电离能、比热、电导、光电子谱、穆斯堡尔谱等等,而且它还可以帮助人们预言许多新的物理现象和物理规律。

密度泛函计算的一些结果能够与实验直接进行比较,一些应用程序的发展乃至商业软件的发布,导致了基于密度泛函理论的第一原理计算方法的广泛应用。

密度泛函理论(DFT)为第一性原理中的一类,在物理系、化学、材料科学以及其他工程领域中,密度泛函理论(DFT)及其计算已经快速发展成为材料建模模拟的一种“标准工具”。

密度泛函理论可以计算预测固体的晶体结构、晶格参数、能带结构、态密度(DOS)、光学性能、磁性能以及原子集合的总能等等。

3第一性原理怎么用?目前我所学到的利用第一性原理的软件为Material Studio、V ASP软件。

其中Materials Studio(简称MS)是专门为材料科学领域研究者开发的一款可运行在PC上的模拟软件。

使化学及材料科学的研究者们能更方便地建立三维结构模型,并对各种晶体、无定型以及高分子材料的性质及相关过程进行深入的研究。

第一性原理计算方法讲义

第一性原理计算方法讲义

第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现,使得人们对于物质微观结构的认识日益深入。

原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。

目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。

1964年,Hohenberg和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT)。

它建立在非均匀作为基本变量。

1965年,Kohn和Sham 电子气理论基础之上,以粒子数密度()r提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

对应。

直接的结果是导致任意可观察量的算符d的基态期望值与基态的电子密度存在唯一泛函:第二定理:对于为哈密顿量片的情况,基态总能的泛函M的形式如下:[p]=〈叫f + v\力+〈叫q平〉(1.67)其中Hohenberg-Kohn密度泛函[p]对任何多电子体系都适用。

Ey…[p\在厂`所对应的基态密度下达到最小值(等于基态的总能)。

以上定理思想的三个关键词为可逆性(一一对应,p<r^V找丨)、普适性和变分方法(最小值)。

下面对三个关键词进行阐述:第一,基态密度与外源的势一一对应。

很显然,对于一个给定的多电子体系存在一个对应的外源势,通过哈密顿量(1.65)和薛定愕方程得到一个唯一的基态多粒子波函数。

从该波函数出发,容易求得相应的电子密度。

因此,任一外源势都可以通过明确的方法得到一个唯一的与之对应的基态密度。

并且,基态密度中包含信息同波函数一样多,包括所有可能知道的关于原子、分子和固体的信息。

所有可观测量都能用唯一的方式从密度中得到,也就是说,均为密度的泛函。

第二,尸“^厂]的普适性。

等式(1.68)可以容易地通过密度算符写出,正印证了在得到基态密度后,外源势对总能的贡献可以精确地计算得到。

Hohenberg-Kohn泛函的表达式不能确切给出,但是,由于不含任何原子以及其位置的信息,所以对任一多电子体系均适用。

这意味着,原则上,尸的表达式适用于所有原子,分子或固体体系。

第三,第二定理使得可以通过变分法求得体系的基态密度。

能使得五`[p]最小的密度正是对应外源势厂⑶(F)的基态密度。

1.3.2.2 Kohn-Sham 方程发表于1965年的Kohn-Sham方程[8]将密度泛函理论变为了实际可操作的工具,求解基态密度成为了实际可用方法。

首先,需要重写Hohenberg-Kohn泛函。

定义关联能(correlation energy)为存在于精确解中,但不存在于Hartree-Fock解中的总能的一部分。

总能泛函五,[p]和分别对应于精确解和Hartree-Fock 哈密顿量:E^=T + V其中r和r分别是精确解的动能和电子之间相互作用的势能泛函,7?是无相互作用电子气的动能泛函,)?代表Hartree贡献,Fx代表交换作用的贡献。

用(1.69)式减(1.70)式,可以得到关联作用贡献的表达式为:IVc=T-T。

(1,71)交换作用的贡献被定义为存在于Hartree-Fock解中,但不存在于Hartree解中的总能的一部分。

显然,Hartree泛函可以由下式给出Eh=T,+V?(1.72)从而得到V^ = V-V?综上,可以将Hohenberg-Kohn泛函重写为FfiK =T -\-V + T^-T^=ro+F+(r -『0)= TQ + V + V^ + V?-Vff(1.74)= TO+Vh+KH^-^H)其中是交换关联泛函,但是其具体形式不得而知。

若假设知道其具体形式,则可精确求解能量的泛函:Ek? [P] = TAPVVh[PV Kc [P] + 厂W(1-75)如此,可以用Hohenberg-Kohn第二定理求解基态密度。

(1.75)式也可以理解为在两个外源势作用下:一个来自原子核,另一个来自交换关联效应,无相互作用的经典电子气的能量泛函。

相应的哈密顿量被称为Kohn-Sham哈密顿量,可以如下写出片卢H+t+t`方2g2 ^(P')(1.76)=-—Vf+ — f dF ` + 厂,c + 厂2m^ `J 丨厂一厂丨|"`其中交换关联势由泛函的导数给出Y(1.77)5pKohn-Sham理论可以如下正式给出:任一N电子系统精确的基态密度是=⑷(1.78其中单粒子波函数於(?)是Kohn-Sham方程的TV□个最低能态的解^KS^i :左M(1_ 79)到此,求解基态密度不需要再运用Hohenberg-Kohn第二定理,只需要求解熟悉的类似薛定愕方程的单粒子方程。

通过替代常规的薛定愕方程,避免了需要求解由于电子间相互作用导致的含有親合的复杂系统的微分方程。

需要注意的是,(1.79)式中的单粒子波函并不是电子的真实波函数。

它们描述的是数学的准粒子,并没有实际的物理意义。

只有准粒子的密度与真实体系的电子密度相同。

单粒子能量A也不是单电子能量。

入:/Vi (少构造;^可以看到,Hartree算符以及交换关联作用算符厂?均依赖于密度/?(?),该密度由直接由要求解的波函数^!>,00所决定。

因此,陷入了一个自洽问题:等式的解(<zUF))决定原始的等式(//ks中的^^和P。

),然而,在解未知情况下方程不能被给出。

如图1.5所示,解决上述问题,可以先猜测一个密度P。

,由此构造哈密顿量Hfu。

可以求解得到//£,,的本征值和本征函数^>1,从可以得到/?,。

很可能P,与猜测的P。

不同。

然后继续从P,构造//fn,求得广2,……如此循环下去,最后会收敛到密度P,,由它构造的方程再次求得P,,那么,该密度则为所要求解的值。

1.3.2.3交换关联泛函为了求解多电子问题,首先采用了波恩-奥本海默近似,然后通过Kohn-Sham方法,将相互稱合的多体问题变成了求解虚拟系统中无相互作用单粒子在外源势中的运动,其粒子的密度与真实的多电子系统的密度一致,并且,可以通过自洽方法求得该密度。

然而,在这一过程中,并不知道交换关联泛函的具体形式,因此,必须对其形式进行猜测,给出进一步的近似。

其中,被广泛使用的交换关联泛函近似叫做局域密度近似(Local DensityApproximation, LDA)。

该近似认为交换关联泛函的形式如下:ET=\p{r)s,c[p{r))dr(1.80)其中,自由电子气函数的数值形式是可以求得的。

该近似认为特定密度对交换关联能的贡献可以被分成许多具有恒定密度的小空间。

每一个小空间对总的交换关联能的贡献等于同等体积下均匀电子气的贡献,且均匀电子气的密度总体上与原材料在该体积的密度一致。

局域密度近似被认为适用于电荷密度变化缓慢的系统,事实上,计算发现对于其它真实系统也足够精确。

在局域密度近似下,各个极小体积对交换关联能的贡献只取决于该体积内的局域密度,因此,自然地想到可以改进该近似,使得对交换关联能的贡献还取决于邻近体积。

也就是说,可以考虑密度的梯度。

如此得到的近似被称为广义梯度近似(Generalized Gradient Approximation, GGA)=I p{7)slr (p{f))F,, (p(r),|Vp(F)|)c/r其中Fa是无量纲的标量泛函,是均匀电子气的交换关联能。

通过对以各阶梯度为参数进行多项式展开,将交换关联能与密度的梯度产生联系。

因此,该近似更适用于处理电子密度不均匀的原子、分子体系。

目前,已经构造出许多形式的GGA泛函,包括B88 (Becke) [9],PW91 (Perdew和Wang)[10]和PBE (Perdew, Burke 和Enzerhof)["]。

1.3.3方程求解按照密度泛函理论,多电子问题被最终近似为单电子问题,通过求解如下的Kohn-Sham方程即可得到体系的基态密度二士+ + r士⑴(1-82)方程的求解,使用的数学方法是寻找一组基函数以及系数将上式中的展开波函数属于一个无穷维的函数空间,因此,P在理论上应该是无限的。

实际中,会将波函数用有限个数的基函数展开。

虽然,有限组的基函数永远不能确切描述但是,可以尽量找寻接近一组?^??的基函数。

一般地,存在三种选择基函数的方法,分别是:>平面波和格点方法(Plane waves and grids)。

这种方法采用平面波作为基函数展开,是解微分方程的通用方法,比如薛定谔方程和泊松方程。

平面波尤其适用于周期性的晶体,釆用该方法可以对所求解的问题提供很多直观的理解,以及一些利于实际计算的简单算法。

另一方面,格点方法极为适合实空间中的有限系统,在许多科学和工程领域均有应用。

现代电4^结构的计算算法用到平面波和包含快速傅立叶变化的格点算法。

>局域轨道方法(Localized orbitals)。

这种方法采用局域的原子轨道作为基底函数,在化学中被广泛应用。

相对于平面波方法,其基函数的数目会极大减少,因此,特别适用于较大体系的计算模拟,比如大分子,以及纳米器件设计。

>缀加函数方法(Augmentedfunctions)。

这种方法将空间分为两个区域,一个为原子位置处的球形区域内,另一个则是各个球形之间的区域。

球形区域之间的区域,波函数变化平缓可以用平面波或者变化平缓的函数表示,在球形区域内部波函数变化剧烈可以用径向函数乘以球谐函数表示,两种函数应在球形边界处连续。

原问题转化为寻找满足边界条件的函数。

假定基函数己经确定,则可以处理(1.82)式的本征值问题,将(1.83)式带入(1.82)式中,并且左乘〈^>^|(丨=1,...,户),得到通过对角化哈密顿量对应的矩阵,可以得到尸个本征值以及对应的P个用所选基函数表示的本征函数,从而完成对方程(1.82)的求解。

相关文档
最新文档