24小时热门版块排行榜

返回列表

【奖励】本帖被评价126次，作者identation增加金币 111.05 个

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

identation

金虫 (著名写手)

[资源] 【原创】第一原理计算的一些心得

把自己关于第一原理的一些理解写出：
（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

回复此楼

» 收录本帖的淘帖专辑推荐

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

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人回复

1楼2010-07-04 14:21:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★★★★★ 五星级,优秀推荐

非常好

但对于这句话的后半部分存疑（“模守恒势的散射特性和全电子相同，因此一般红外，拉曼等光谱的计算需要用模守恒势。）”
模守恒势的散射特性和全电子相同。此处的散射性质好像是指的电子散射，而不是光子散射。而且，USPP已经广泛用于红外、拉曼的频率计算，还有红外强度的计算。

[ Last edited by xirainbow on 2010-7-4 at 19:36 ]