VASPsol隱式溶劑模型實操指南從真空DFT到溶劑化計算的完整避坑路線【免費下載鏈接】VASPsolSolvation model for the plane wave DFT code VASP.項目地址: https://gitcode.com/gh_mirrors/va/VASPsol如果有一件事能讓一位材料計算研究者的周末徹底泡湯那就是真空里算得完美無缺的結(jié)構(gòu)被導(dǎo)師一句那在水溶液里呢問得啞口無言。VASPsol 正是為這類場景準備的——它是面向平面波 DFT 代碼 VASP 的隱式溶劑模型實現(xiàn)用連續(xù)介質(zhì)代替成千上萬個顯式水分子把靜電、空化、色散三種溶劑效應(yīng)并入自洽場迭代。本文不堆概念直接按決策 → 安裝 → 調(diào)參 → 踩坑 → 復(fù)現(xiàn)的順序帶你走完從真空計算遷移到溶劑化計算的完整路徑。一個讓我重新認識溶劑效應(yīng)的下午故事要從我接手一個電化學(xué)還原 CO? 的課題說起。催化劑表面的吸附能在真空模型下算得很漂亮數(shù)據(jù)也擬合得很好直到組會上有人問你的反應(yīng)發(fā)生在水溶液里真空的結(jié)果有多少參考價值我當時的直覺是補一個顯式溶劑模型往盒子塞幾十個水分子再跑一遍分子動力學(xué)取樣。結(jié)果很現(xiàn)實——機時翻了幾倍收斂曲線像過山車一個體系折騰兩周還沒跑完。后來才意識到我需要的不是更真實的溶劑而是一個能抓住溶劑核心物理效應(yīng)的低成本替代方案。這就是隱式溶劑模型的定位不逐顆模擬溶劑分子而是把溶劑看作一個具有介電常數(shù)、表面張力的連續(xù)介質(zhì)溶質(zhì)周圍自動生成一個空腔空腔之外就是均勻的溶劑。VASPsol 把這件事直接做進了 VASP 的自洽循環(huán)里。顯式與隱式一場關(guān)于性價比的理性選擇先別急著動手花兩分鐘想清楚該走哪條路。兩種方案沒有絕對優(yōu)劣只有匹配不匹配。對比維度顯式溶劑模型VASPsol 隱式溶劑模型建模方式數(shù)百個溶劑分子需平衡構(gòu)型與動力學(xué)取樣連續(xù)介質(zhì) 介電空腔無額外自由度計算成本數(shù)倍至數(shù)十倍機時取樣要求高比真空計算只多約 30%–50% 開銷適用體系溶質(zhì)–溶劑短程相互作用主導(dǎo)的場景周期體系、表面、界面、反應(yīng)路徑能拿到什么溶劑構(gòu)型、氫鍵網(wǎng)絡(luò)、動力學(xué)信息溶劑化能、束縛電荷、靜電響應(yīng)上手難度需要力場、平衡、統(tǒng)計處理鏈條長改幾個 INCAR 參數(shù)即可決策指南可以濃縮成三句話要研究氫鍵網(wǎng)絡(luò)和溶劑重組動力學(xué)請選顯式模型要算吸附能、反應(yīng)能壘、溶劑化能這些能量學(xué)問題隱式模型是更劃算的起點如果你的體系是金屬或半導(dǎo)體表面這類周期性體系隱式模型幾乎就是唯一現(xiàn)實的選擇——這正是 VASPsol 設(shè)計時的重點場景。先看懂它在算什么再動手配置VASPsol 的物理內(nèi)核并不神秘它把溶劑效應(yīng)拆成三股力量靜電相互作用溶劑的高介電常數(shù)會屏蔽溶質(zhì)的電荷分布誘導(dǎo)出束縛電荷這部分由相對介電常數(shù)EB_K控制空化能在溶劑中挖出一個容納溶質(zhì)的空腔需要做功由表面張力參數(shù)TAU控制色散相互作用溶質(zhì)與溶劑分子之間的范德華吸引也一并計入。關(guān)鍵的一點是這三者不是算完一遍就結(jié)束的。束縛電荷會反過來改變?nèi)苜|(zhì)的電子結(jié)構(gòu)電子結(jié)構(gòu)又決定束縛電荷的分布——所以 VASPsol 在每步電子迭代里都要重新求解這套耦合方程直到變化小于收斂標準EDIFFSOL默認取EDIFF/100。理解了這條自洽鏈后面遇到收斂問題時的排查方向就清楚了。把 VASPsol 裝進 VASP一張分版本對照表安裝方式完全取決于你的 VASP 版本對照下表操作即可。VASP 版本要做的操作注意事項5.2.12 / 5.3.3 / 5.3.5應(yīng)用 patches/ 下對應(yīng)補丁拷貝solvation.F調(diào)整 Makefile確保solvation.o排在pot.o之前然后重新編譯≥ 5.4.1只需把 src/solvation.F 拷貝到 VASP 的src/目錄最省事的一條路線≥ 5.4.4在上述基礎(chǔ)上給編譯器加-Dsol_compat預(yù)處理選項新版 VASP 重構(gòu)了余誤差函數(shù)定義不加會編譯報錯≥ 6.1.0應(yīng)用VASPsol6.patch并確保.objects中solvation.o在pot.o之前補丁同時修正了電解質(zhì)模型需要的FERMI_SHIFT輸出對最常見的中高版本安裝本質(zhì)上就是一條命令git clone https://gitcode.com/gh_mirrors/va/VASPsol cp src/solvation.F /path/to/vasp.5.4.x/src/然后按官方流程重新編譯。編譯通過后先跑一次測試計算確認 stdout 里出現(xiàn) All solvation model paramters are set 這類輸出再正式開工。30 分鐘跑通第一個溶劑化計算一個容易忽略的官方建議是永遠從真空波函數(shù)出發(fā)。完整流程分四步。先做一次真空計算結(jié)構(gòu)優(yōu)化或靜態(tài)均可但 INCAR 里務(wù)必寫明LWAVE .TRUE.保存好 WAVECAR在 INCAR 中加入溶劑化參數(shù)最簡配置只有兩行LSOL .TRUE. # 啟用隱式溶劑默認即為水 EB_K 80.0 # 溶劑的相對介電常數(shù)水的參考值 78.4~80設(shè)置ISTART 1直接從真空 WAVECAR 續(xù)算運行結(jié)束后從輸出中讀取溶劑化能vasp_std vasp.out grep SOL: OUTCAR輸出形如SOL: 1 0.12345E01 0.23456E00 0.14691E01 45一行里的三個能量依次是靜電貢獻、空化能貢獻和總?cè)軇┗軉挝?eV。想換溶劑就改EB_K想關(guān)掉空化能就設(shè)TAU 0——第一個溶劑化計算就這么跑通了。參數(shù)旋鈕怎么擰一張決策表VASPsol 的 INCAR 參數(shù)不多但每個旋鈕都有明確的適用場景。我把它們整理成一張什么時候動它的決策表。參數(shù)控制什么什么場景下需要調(diào)整參考取值LSOL溶劑化總開關(guān)所有溶劑化計算.TRUE.EB_K溶劑介電常數(shù)換溶劑時必改水 78.4有機溶劑 10–40TAU表面張力只想考察靜電效應(yīng)時忽略空化能設(shè) 0LAMBDA_D_KDebye 長度?模擬電解質(zhì)溶液時0.1 M 約 3.0 ?EDIFFSOL溶劑化迭代收斂標準大體系提速時放寬初算 1E-5精算 1E-6LRHOB是否輸出束縛電荷密度想可視化電荷分布時.TRUE.SIGMA_K/NC_K空腔寬度與開啟密度一般不動默認擬合自水默認即可 一條貫穿始終的經(jīng)驗空化能的數(shù)值強烈依賴實空間網(wǎng)格能否分辨空腔表面。官方文檔反復(fù)強調(diào)PREC Accurate且ENCUT要留足余量這不是建議是要求。真空計算里夠用的截斷能搬到溶劑化計算往往要再抬高一截。踩坑記錄四個讓我折騰到深夜的問題與其抽象地講 FAQ不如分享我真實排過的四個坑按出現(xiàn)頻率排序?？右荒芰刻鴣硖ゾ褪遣皇諗?。癥狀是電子迭代步數(shù)拉滿仍不收斂或能量在小數(shù)點后第三位來回振蕩。我排查的順序是先看ENCUT是否足夠、PREC是否設(shè)為 Accurate再從真空 WAVECAR 起算ISTART 1而非冷啟動最后才考慮放寬EDIFFSOL。前兩步解決了 80% 的情況?？佣Y(jié)果跟別人對不上差一大截。十有八九是沒從真空 WAVECAR 續(xù)算。隱式溶劑的自洽迭代很依賴一個合理的初始電荷密度真空收斂波函數(shù)就是最好的起點。務(wù)必確認真空計算時LWAVE .TRUE.溶劑化計算時ISTART 1?？尤_了 LRHOB 之后大體系慢到懷疑人生。LRHOB .TRUE.會把束縛電荷密度寫入 RHOB 文件但實現(xiàn)方式?jīng)Q定了每一步 SCF 迭代都會寫一次文件。體系一大I/O 開銷直接吞掉性能。正確姿勢是先不帶它算完并收斂再開一次靜態(tài)單點計算專門輸出 RHOB用于可視化分析?？铀碾娊赓|(zhì)模型算出來電位參考點不對。這是最隱蔽的一個。VASP 的靜電參考電位設(shè)定是整個模擬盒平均為零而不是電解質(zhì)體相區(qū)域為零。開啟LAMBDA_D_K后需要把 stdout 里打印的FERMI_SHIFT加到費米能級上如果你還平移了參考電位別忘了對體系能量補一個Q·V修正其中Q是凈電荷V是平移量。不補這兩步帶電體系的絕對能量沒有可比性。兩個值得先復(fù)現(xiàn)的場景場景一分子的溶劑化能。examples/H2O/ 提供了水分子在真空與溶劑兩種條件下的完整輸入。復(fù)現(xiàn)路徑先跑 Vacuum 目錄下的真空計算再用 Solvation 目錄的 INCAR 從 WAVECAR 續(xù)算最后做減法ΔE_solv E(溶劑中) ? E(真空中)這個差值就是分子溶劑化能也是理解溶劑到底穩(wěn)不穩(wěn)定這個分子的第一塊敲門磚。場景二電解質(zhì)溶液中的表面過程。examples/CO/ 和 examples/PbS_100/ 展示了表面吸附體系的配置方式。若涉及電化學(xué)環(huán)境在 INCAR 里加入LAMBDA_D_K即可切換到線性化 Poisson-Boltzmann 電解質(zhì)模型LSOL .TRUE. EB_K 80.0 LAMBDA_D_K 3.0 # 約對應(yīng) 0.1 M 電解質(zhì)溶液注意 Debye 長度隨濃度縮放0.01 M 約 9.6 ?0.1 M 約 3.0 ?1.0 M 約 0.96 ?。濃度越高屏蔽越強需要的真空層和盒長也隨之變化。從復(fù)制代碼到獨當一面一條可復(fù)制的學(xué)習(xí)路線圖最后給想把 VASPsol 用扎實的讀者一條五階段路線每個階段都有明確的驗收標準跑通示例把 examples/ 里的 H2O 案例從頭跑一遍能正確讀出SOL:行的三個能量——驗收標準是結(jié)果與參考輸出量級一致復(fù)現(xiàn)文獻數(shù)值用論文里的分子或表面體系復(fù)現(xiàn)溶劑化能——這一步能幫你校準自己的ENCUT、EDIFFSOL選擇遷移到自有體系把自己的結(jié)構(gòu)從真空遷移到溶劑化流程理解EB_K與TAU的取舍進入電解質(zhì)領(lǐng)域開啟LAMBDA_D_K學(xué)會用FERMI_SHIFT和Q·V修正處理帶電體系進階分析用LRHOB輸出束縛電荷把靜電響應(yīng)可視化跟實驗結(jié)果對照。每一步之間不必追求一步到位。隱式溶劑模型的魅力在于它把溶劑效應(yīng)這個龐大的問題壓縮成了幾個可以逐個吃透的參數(shù)你完全可以從最小的分子開始一步步走到帶電表面。引用方式如果 VASPsol 在你的研究中發(fā)揮了作用請引用下面兩篇論文第二篇對應(yīng)電解質(zhì)模型部分article{VASPsol2014-Dielectric, title {Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways}, author {K. Mathew and R. Sundararaman and K. Letchworth-Weaver and T. A. Arias and R. G. Hennig}, journal {J. Chem. Phys.}, volume 140, pages {084106}, year 2014 } article{VASPsol2019-Electrolyte, title {Implicit self-consistent electrolyte model in plane-wave density-functional theory}, author {K. Mathew and V. S. C. Kolluru and S. Mula and S. N. Steinmann and R. G. Hennig}, journal {J. Chem. Phys.}, volume 151, pages {234101}, year 2019 }更詳細的使用說明見官方文檔 docs/USAGE.md?，F(xiàn)在從克隆倉庫、編譯 VASP 開始把第一個水分子泡進溶劑里吧——從真空到溶液可能只差三行 INCAR?！久赓M下載鏈接】VASPsolSolvation model for the plane wave DFT code VASP.項目地址: https://gitcode.com/gh_mirrors/va/VASPsol創(chuàng)作聲明：本文部分內(nèi)容由AI輔助生成（AIGC），僅供參考 VASPsol隱式溶劑模型實操指南：從真空DFT到溶劑化計算的完整避坑路線
VASPsol隱式溶劑模型實操指南從真空DFT到溶劑化計算的完整避坑路線【免費下載鏈接】VASPsolSolvation model for the plane wave DFT code VASP.項目地址: https://gitcode.com/gh_mirrors/va/VASPsol如果有一件事能讓一位材料計算研究者的周末徹底泡湯那就是真空里算得完美無缺的結(jié)構(gòu)被導(dǎo)師一句那在水溶液里呢問得啞口無言。VASPsol 正是為這類場景準備的——它是面向平面波 DFT 代碼 VASP 的隱式溶劑模型實現(xiàn)用連續(xù)介質(zhì)代替成千上萬個顯式水分子把靜電、空化、色散三種溶劑效應(yīng)并入自洽場迭代。本文不堆概念直接按決策 → 安裝 → 調(diào)參 → 踩坑 → 復(fù)現(xiàn)的順序帶你走完從真空計算遷移到溶劑化計算的完整路徑。一個讓我重新認識溶劑效應(yīng)的下午故事要從我接手一個電化學(xué)還原 CO? 的課題說起。催化劑表面的吸附能在真空模型下算得很漂亮數(shù)據(jù)也擬合得很好直到組會上有人問你的反應(yīng)發(fā)生在水溶液里真空的結(jié)果有多少參考價值我當時的直覺是補一個顯式溶劑模型往盒子塞幾十個水分子再跑一遍分子動力學(xué)取樣。結(jié)果很現(xiàn)實——機時翻了幾倍收斂曲線像過山車一個體系折騰兩周還沒跑完。后來才意識到我需要的不是更真實的溶劑而是一個能抓住溶劑核心物理效應(yīng)的低成本替代方案。這就是隱式溶劑模型的定位不逐顆模擬溶劑分子而是把溶劑看作一個具有介電常數(shù)、表面張力的連續(xù)介質(zhì)溶質(zhì)周圍自動生成一個空腔空腔之外就是均勻的溶劑。VASPsol 把這件事直接做進了 VASP 的自洽循環(huán)里。顯式與隱式一場關(guān)于性價比的理性選擇先別急著動手花兩分鐘想清楚該走哪條路。兩種方案沒有絕對優(yōu)劣只有匹配不匹配。對比維度顯式溶劑模型VASPsol 隱式溶劑模型建模方式數(shù)百個溶劑分子需平衡構(gòu)型與動力學(xué)取樣連續(xù)介質(zhì) 介電空腔無額外自由度計算成本數(shù)倍至數(shù)十倍機時取樣要求高比真空計算只多約 30%–50% 開銷適用體系溶質(zhì)–溶劑短程相互作用主導(dǎo)的場景周期體系、表面、界面、反應(yīng)路徑能拿到什么溶劑構(gòu)型、氫鍵網(wǎng)絡(luò)、動力學(xué)信息溶劑化能、束縛電荷、靜電響應(yīng)上手難度需要力場、平衡、統(tǒng)計處理鏈條長改幾個 INCAR 參數(shù)即可決策指南可以濃縮成三句話要研究氫鍵網(wǎng)絡(luò)和溶劑重組動力學(xué)請選顯式模型要算吸附能、反應(yīng)能壘、溶劑化能這些能量學(xué)問題隱式模型是更劃算的起點如果你的體系是金屬或半導(dǎo)體表面這類周期性體系隱式模型幾乎就是唯一現(xiàn)實的選擇——這正是 VASPsol 設(shè)計時的重點場景。先看懂它在算什么再動手配置VASPsol 的物理內(nèi)核并不神秘它把溶劑效應(yīng)拆成三股力量靜電相互作用溶劑的高介電常數(shù)會屏蔽溶質(zhì)的電荷分布誘導(dǎo)出束縛電荷這部分由相對介電常數(shù)EB_K控制空化能在溶劑中挖出一個容納溶質(zhì)的空腔需要做功由表面張力參數(shù)TAU控制色散相互作用溶質(zhì)與溶劑分子之間的范德華吸引也一并計入。關(guān)鍵的一點是這三者不是算完一遍就結(jié)束的。束縛電荷會反過來改變?nèi)苜|(zhì)的電子結(jié)構(gòu)電子結(jié)構(gòu)又決定束縛電荷的分布——所以 VASPsol 在每步電子迭代里都要重新求解這套耦合方程直到變化小于收斂標準EDIFFSOL默認取EDIFF/100。理解了這條自洽鏈后面遇到收斂問題時的排查方向就清楚了。把 VASPsol 裝進 VASP一張分版本對照表安裝方式完全取決于你的 VASP 版本對照下表操作即可。VASP 版本要做的操作注意事項5.2.12 / 5.3.3 / 5.3.5應(yīng)用 patches/ 下對應(yīng)補丁拷貝solvation.F調(diào)整 Makefile確保solvation.o排在pot.o之前然后重新編譯≥ 5.4.1只需把 src/solvation.F 拷貝到 VASP 的src/目錄最省事的一條路線≥ 5.4.4在上述基礎(chǔ)上給編譯器加-Dsol_compat預(yù)處理選項新版 VASP 重構(gòu)了余誤差函數(shù)定義不加會編譯報錯≥ 6.1.0應(yīng)用VASPsol6.patch并確保.objects中solvation.o在pot.o之前補丁同時修正了電解質(zhì)模型需要的FERMI_SHIFT輸出對最常見的中高版本安裝本質(zhì)上就是一條命令git clone https://gitcode.com/gh_mirrors/va/VASPsol cp src/solvation.F /path/to/vasp.5.4.x/src/然后按官方流程重新編譯。編譯通過后先跑一次測試計算確認 stdout 里出現(xiàn) All solvation model paramters are set 這類輸出再正式開工。30 分鐘跑通第一個溶劑化計算一個容易忽略的官方建議是永遠從真空波函數(shù)出發(fā)。完整流程分四步。先做一次真空計算結(jié)構(gòu)優(yōu)化或靜態(tài)均可但 INCAR 里務(wù)必寫明LWAVE .TRUE.保存好 WAVECAR在 INCAR 中加入溶劑化參數(shù)最簡配置只有兩行LSOL .TRUE. # 啟用隱式溶劑默認即為水 EB_K 80.0 # 溶劑的相對介電常數(shù)水的參考值 78.4~80設(shè)置ISTART 1直接從真空 WAVECAR 續(xù)算運行結(jié)束后從輸出中讀取溶劑化能vasp_std vasp.out grep SOL: OUTCAR輸出形如SOL: 1 0.12345E01 0.23456E00 0.14691E01 45一行里的三個能量依次是靜電貢獻、空化能貢獻和總?cè)軇┗軉挝?eV。想換溶劑就改EB_K想關(guān)掉空化能就設(shè)TAU 0——第一個溶劑化計算就這么跑通了。參數(shù)旋鈕怎么擰一張決策表VASPsol 的 INCAR 參數(shù)不多但每個旋鈕都有明確的適用場景。我把它們整理成一張什么時候動它的決策表。參數(shù)控制什么什么場景下需要調(diào)整參考取值LSOL溶劑化總開關(guān)所有溶劑化計算.TRUE.EB_K溶劑介電常數(shù)換溶劑時必改水 78.4有機溶劑 10–40TAU表面張力只想考察靜電效應(yīng)時忽略空化能設(shè) 0LAMBDA_D_KDebye 長度?模擬電解質(zhì)溶液時0.1 M 約 3.0 ?EDIFFSOL溶劑化迭代收斂標準大體系提速時放寬初算 1E-5精算 1E-6LRHOB是否輸出束縛電荷密度想可視化電荷分布時.TRUE.SIGMA_K/NC_K空腔寬度與開啟密度一般不動默認擬合自水默認即可 一條貫穿始終的經(jīng)驗空化能的數(shù)值強烈依賴實空間網(wǎng)格能否分辨空腔表面。官方文檔反復(fù)強調(diào)PREC Accurate且ENCUT要留足余量這不是建議是要求。真空計算里夠用的截斷能搬到溶劑化計算往往要再抬高一截。踩坑記錄四個讓我折騰到深夜的問題與其抽象地講 FAQ不如分享我真實排過的四個坑按出現(xiàn)頻率排序?？右荒芰刻鴣硖ゾ褪遣皇諗?。癥狀是電子迭代步數(shù)拉滿仍不收斂或能量在小數(shù)點后第三位來回振蕩。我排查的順序是先看ENCUT是否足夠、PREC是否設(shè)為 Accurate再從真空 WAVECAR 起算ISTART 1而非冷啟動最后才考慮放寬EDIFFSOL。前兩步解決了 80% 的情況?？佣Y(jié)果跟別人對不上差一大截。十有八九是沒從真空 WAVECAR 續(xù)算。隱式溶劑的自洽迭代很依賴一個合理的初始電荷密度真空收斂波函數(shù)就是最好的起點。務(wù)必確認真空計算時LWAVE .TRUE.溶劑化計算時ISTART 1?？尤_了 LRHOB 之后大體系慢到懷疑人生。LRHOB .TRUE.會把束縛電荷密度寫入 RHOB 文件但實現(xiàn)方式?jīng)Q定了每一步 SCF 迭代都會寫一次文件。體系一大I/O 開銷直接吞掉性能。正確姿勢是先不帶它算完并收斂再開一次靜態(tài)單點計算專門輸出 RHOB用于可視化分析?？铀碾娊赓|(zhì)模型算出來電位參考點不對。這是最隱蔽的一個。VASP 的靜電參考電位設(shè)定是整個模擬盒平均為零而不是電解質(zhì)體相區(qū)域為零。開啟LAMBDA_D_K后需要把 stdout 里打印的FERMI_SHIFT加到費米能級上如果你還平移了參考電位別忘了對體系能量補一個Q·V修正其中Q是凈電荷V是平移量。不補這兩步帶電體系的絕對能量沒有可比性。兩個值得先復(fù)現(xiàn)的場景場景一分子的溶劑化能。examples/H2O/ 提供了水分子在真空與溶劑兩種條件下的完整輸入。復(fù)現(xiàn)路徑先跑 Vacuum 目錄下的真空計算再用 Solvation 目錄的 INCAR 從 WAVECAR 續(xù)算最后做減法ΔE_solv E(溶劑中) ? E(真空中)這個差值就是分子溶劑化能也是理解溶劑到底穩(wěn)不穩(wěn)定這個分子的第一塊敲門磚。場景二電解質(zhì)溶液中的表面過程。examples/CO/ 和 examples/PbS_100/ 展示了表面吸附體系的配置方式。若涉及電化學(xué)環(huán)境在 INCAR 里加入LAMBDA_D_K即可切換到線性化 Poisson-Boltzmann 電解質(zhì)模型LSOL .TRUE. EB_K 80.0 LAMBDA_D_K 3.0 # 約對應(yīng) 0.1 M 電解質(zhì)溶液注意 Debye 長度隨濃度縮放0.01 M 約 9.6 ?0.1 M 約 3.0 ?1.0 M 約 0.96 ?。濃度越高屏蔽越強需要的真空層和盒長也隨之變化。從復(fù)制代碼到獨當一面一條可復(fù)制的學(xué)習(xí)路線圖最后給想把 VASPsol 用扎實的讀者一條五階段路線每個階段都有明確的驗收標準跑通示例把 examples/ 里的 H2O 案例從頭跑一遍能正確讀出SOL:行的三個能量——驗收標準是結(jié)果與參考輸出量級一致復(fù)現(xiàn)文獻數(shù)值用論文里的分子或表面體系復(fù)現(xiàn)溶劑化能——這一步能幫你校準自己的ENCUT、EDIFFSOL選擇遷移到自有體系把自己的結(jié)構(gòu)從真空遷移到溶劑化流程理解EB_K與TAU的取舍進入電解質(zhì)領(lǐng)域開啟LAMBDA_D_K學(xué)會用FERMI_SHIFT和Q·V修正處理帶電體系進階分析用LRHOB輸出束縛電荷把靜電響應(yīng)可視化跟實驗結(jié)果對照。每一步之間不必追求一步到位。隱式溶劑模型的魅力在于它把溶劑效應(yīng)這個龐大的問題壓縮成了幾個可以逐個吃透的參數(shù)你完全可以從最小的分子開始一步步走到帶電表面。引用方式如果 VASPsol 在你的研究中發(fā)揮了作用請引用下面兩篇論文第二篇對應(yīng)電解質(zhì)模型部分article{VASPsol2014-Dielectric, title {Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways}, author {K. Mathew and R. Sundararaman and K. Letchworth-Weaver and T. A. Arias and R. G. Hennig}, journal {J. Chem. Phys.}, volume 140, pages {084106}, year 2014 } article{VASPsol2019-Electrolyte, title {Implicit self-consistent electrolyte model in plane-wave density-functional theory}, author {K. Mathew and V. S. C. Kolluru and S. Mula and S. N. Steinmann and R. G. Hennig}, journal {J. Chem. Phys.}, volume 151, pages {234101}, year 2019 }更詳細的使用說明見官方文檔 docs/USAGE.md?，F(xiàn)在從克隆倉庫、編譯 VASP 開始把第一個水分子泡進溶劑里吧——從真空到溶液可能只差三行 INCAR?！久赓M下載鏈接】VASPsolSolvation model for the plane wave DFT code VASP.項目地址: https://gitcode.com/gh_mirrors/va/VASPsol創(chuàng)作聲明：本文部分內(nèi)容由AI輔助生成（AIGC），僅供參考