| 查看: 26229 | 回复: 160 | ||
| 【奖励】 本帖被评价126次，作者identation增加金币 111.05 个 | ||
[资源]
【原创】第一原理计算的一些心得
|
||
|
把自己关于第一原理的一些理解写出： （1）第一性原理其实是包括基于密度泛函的从头算和基于Hartree-Fock自洽计算的从头算，前者以电子密度作为基本变量（霍亨伯格-科洪定理），通过求解Kohn-Sham方程，迭代自洽得到体系的基态电子密度，然后求体系的基态性质；后者则通过自洽求解Hartree-Fock方程，获得体系的波函数，求基态性质； 评述：K-S方程的计算水平达到了H-F水平，同时还考虑了电子间的交换关联作用。 （2）关于DFT中密度泛函的Functional，其实是交换关联泛函 包括LDA，GGA，杂化泛函等等 一般LDA为局域密度近似，在空间某点用均匀电子气密度作为交换关联泛函的唯一变量，多数为参数化的CA-PZ方案； GGA为广义梯度近似，不仅将电子密度作为交换关联泛函的变量，也考虑了密度的梯度为变量，包括PBE,PW,RPBE等方案，BLYP泛函也属于GGA； 此外还有一些杂化泛函，B3LYP等。 （3）关于赝势 在处理计算体系中原子的电子态时，有两种方法，一种是考虑所有电子，叫做全电子法，比如WIEN2K中的FLAPW方法(线性缀加平面波)；此外还有一种方法是只考虑价电子，而把芯电子和原子核构成离子实放在一起考虑，即赝势法，一般赝势法是选取一个截断半径，截断半径以内，波函数变化较平滑，和真实的不同，截断半径以外则和真实情况相同，而且赝势法得到的能量本征值和全电子法应该相同。 赝势包括模守恒和超软，模守恒较硬，一般需要较大的截断能，超软势则可以用较小的截断能即可。另外，模守恒势的散射特性和全电子相同，因此一般红外，拉曼等光谱的计算需要用模守恒势。 赝势的测试标准应是赝势与全电子法计算结果的匹配度，而不是赝势与实验结果的匹配度，因为和实验结果的匹配可能是偶然的。 （4）关于收敛测试 （a）Ecut，也就是截断能，一般情况下，总能相对于不同Ecut做计算，当Ecut增大时总能变化不明显了即可；然而，在需要考虑体系应力时，还需对应力进行收敛测试，而且应力相对于Ecut的收敛要比总能更为苛刻，也就是某个截断能下总能已经收敛了，但应力未必收敛。 （b）K-point，即K网格，一般金属需要较大的K网格，采用超晶胞时可以选用相对较小的K网格，但实际上还是要经过测试。 （5）关于磁性 一般何时考虑自旋呢？举例子，例如BaTiO3中，Ba、Ti和O分别为+2，+4和-2价，离子全部为各个轨道满壳层的结构，就不必考虑自旋了；对于BaMnO3中，由于Mn+3价时d轨道还有电子，但未满，因此需考虑Mn的自旋，至于Ba和O则不必考虑。其实设定自旋就是给定一个原子磁矩的初始值，只在刚开始计算时作为初始值使用，具体的可参照磁性物理。 （6）关于几何优化 包括很多种了，比如晶格常数和原子位置同时优化，只优化原子位置，只优化晶格常数，还有晶格常数和原子位置分开优化等等。 在PRL一篇文章中见到过只优化原子位置，晶格常数用实验值的例子（PRL 100, 186402 (2008)）；也见到过晶格常数先优化，之后固定晶格常数优化原子位置的情况；更多的情况则是Full geometry optimization。 一般情况下，也有不优化几何结构直接计算电子结构的，但是对于缺陷形成能的计算则往往要优化。 （7）关于软件 软件大致分为基于平面波的软件，如CASTEP、PWSCF和ABINIT等等，计算量大概和体系原子数目的三次方相关；还有基于原子轨道线性组合的软件(LCAO)，比如openmx，siesta，dmol等，计算量和体系原子数目相关，一般可模拟较多原子数目的体系。 参考文献： Designing meaningful density functional theory calculations in materials science—a primer |
回复此楼» 收录本帖的淘帖专辑推荐
» 猜你喜欢
299求调剂
已经有5人回复
-
材料专业求调剂
已经有6人回复
-
311求调剂
已经有5人回复
-
317求调剂
已经有13人回复
-
296求调剂
已经有10人回复
-
274求调剂
已经有10人回复
-
22408 344分 求调剂 一志愿 华电计算机技术
已经有4人回复
-
材料专硕英一数二306
已经有7人回复
-
南京大学化学376求调剂
已经有6人回复
-
南昌大学材料专硕311分求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp5.2计算半导体Si光学性质的总结
已经有106人回复
- HSE（GW类似）计算能带总结
已经有51人回复
- 分享vasp的一些可用的免费资源（入门级）
已经有79人回复
- VASP计算时的几个细节设置问题，希望有经验的人过来指点一下！
已经有23人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【原创】PWSCF中的nbnd和degauss
已经有11人回复
- 【原创】DFT计算体会
已经有20人回复
- 【求助成功】怎么样提高NEB计算速度?
已经有19人回复
- 【活动】各种量子力学计算软件使用经验谈（回复就有奖励，讲述详细的奖励5-10金币）
已经有30人回复
- 【其他】第一性计算数据处理经验有奖征集
已经有34人回复
6楼2010-07-04 16:28:23
4楼2010-07-04 15:08:59
5楼2010-07-04 16:01:35
8楼2010-07-04 17:42:48
9楼2010-07-04 18:40:32
10楼2010-07-04 19:34:51
11楼2010-07-04 20:17:56
12楼2010-07-04 22:10:56
13楼2010-07-05 09:37:18
17楼2010-07-05 14:43:58
18楼2010-07-05 14:58:52
19楼2010-07-05 16:26:09
20楼2010-07-05 19:36:27
21楼2010-07-05 21:07:58
22楼2010-07-05 23:02:47
23楼2010-07-05 23:41:48
26楼2010-07-06 13:39:34
27楼2010-07-06 14:03:57
28楼2010-07-06 21:59:07
29楼2010-07-07 06:37:30
31楼2010-07-07 15:10:48
32楼2010-07-08 08:47:52
34楼2010-07-08 16:36:46
35楼2010-07-08 19:07:39
36楼2010-07-08 19:19:16
38楼2010-07-12 08:32:24
39楼2010-07-12 10:34:53
40楼2010-07-12 10:43:40
43楼2010-07-14 21:40:14
44楼2010-07-15 16:59:35
45楼2010-07-16 12:59:02
46楼2010-07-17 00:05:51
47楼2010-07-17 07:26:32
48楼2010-07-17 10:43:25
49楼2010-07-17 21:35:52
50楼2010-07-18 12:31:09
简单回复
hedaors2楼
2010-07-04 14:26
回复
2010-07-04 15:00
回复
carlaty7楼
2010-07-04 17:06
回复
dnaja14楼
2010-07-05 10:22
回复
hn1987051915楼
2010-07-05 10:27
回复
md94119416楼
2010-07-05 13:28
回复
return82724楼
2010-07-06 07:59
回复
yfwanhy25楼
2010-07-06 08:47
回复
dalongmao30楼
2010-07-07 14:48
回复
fhz33楼
2010-07-08 10:55
回复
32875624937楼
2010-07-08 20:12
回复
qwji41楼
2010-07-12 16:55
回复
顶一下。
fansw112942楼
2010-07-13 08:18
回复
