DFT 中的自洽計(jì)算與非自洽計(jì)算
核心區(qū)別在于:是否通過(guò)迭代更新電子密度以滿足自洽條件,二者的計(jì)算目的、流程和參數(shù)設(shè)置截然不同。
一、 核心定義與本質(zhì)區(qū)別
1、自洽計(jì)算(SCF)
自洽計(jì)算的核心目標(biāo)是求解基態(tài)電子密度。
它基于 Kohn-Sham 方程,通過(guò)迭代循環(huán)實(shí)現(xiàn):先假設(shè)一個(gè)初始電子密度,計(jì)算對(duì)應(yīng)的勢(shì)函數(shù),求解單電子方程得到波函數(shù),再由波函數(shù)重新計(jì)算電子密度;重復(fù)此過(guò)程,直到兩次迭代的電子密度(或總能量)的差異小于設(shè)定閾值,達(dá)到自洽收斂。
收斂后得到的電子密度、總能量、電荷密度等,是體系的基態(tài)性質(zhì)。
2、非自洽計(jì)算(Non-SCF)
非自洽計(jì)算不更新電子密度,而是固定電子密度(通常取自 SCF 計(jì)算的收斂結(jié)果),在此基礎(chǔ)上求解 Kohn-Sham 方程,計(jì)算特定需求的物理量。
因?yàn)殡娮用芏炔辉俚拢杂?jì)算速度遠(yuǎn)快于 SCF,且不需要滿足自洽收斂條件。
二、關(guān)鍵差異對(duì)比

三、 典型應(yīng)用場(chǎng)景
1、自洽計(jì)算的應(yīng)用
優(yōu)化晶胞參數(shù)和原子位置(結(jié)構(gòu)弛豫,IBRION≠0時(shí)默認(rèn)包含 SCF);
計(jì)算體系的基態(tài)總能量、形成能、結(jié)合能;
獲得收斂的電荷密度,為后續(xù) Non-SCF 計(jì)算提供輸入。
2、非自洽計(jì)算的應(yīng)用
高精度能帶結(jié)構(gòu)計(jì)算:先通過(guò) SCF 得到基態(tài)密度,再選取高對(duì)稱(chēng) k 點(diǎn)路徑,用 Non-SCF 計(jì)算能帶(需設(shè)置KPOINTS為高對(duì)稱(chēng)路徑);
態(tài)密度(DOS)計(jì)算:使用密集的 k 點(diǎn)網(wǎng)格,固定電荷密度計(jì)算 DOS,提高精度;
電荷密度差分、局域態(tài)密度(LDOS):基于固定密度分析電子分布的變化;
光學(xué)性質(zhì)、介電函數(shù)等進(jìn)階性質(zhì)的計(jì)算。
四、 二者的關(guān)聯(lián)
在 VASP 的實(shí)際計(jì)算中,Non-SCF 計(jì)算幾乎總是依賴(lài)于 SCF 計(jì)算的結(jié)果:
1、先進(jìn)行 SCF 計(jì)算,得到收斂的電荷密度文件(CHGCAR)和波函數(shù)文件(WAVECAR);
2、保持體系結(jié)構(gòu)不變,修改INCAR參數(shù)(如設(shè)置ICHARG=11)、KPOINTS文件(如高對(duì)稱(chēng)路徑),執(zhí)行 Non-SCF 計(jì)算,調(diào)用CHGCAR中的電荷密度進(jìn)行后續(xù)求解。
提交
吸附能與結(jié)合能的理解和計(jì)算
Gaussian 計(jì)算分子偶極矩
如何提升COMSOL模型計(jì)算效率
量子化學(xué):HOMO/LUMO定義到應(yīng)用
量子化學(xué):什么是表面重構(gòu)?

投訴建議