 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3269 | 回復: 12 | |||||||
[交流]
關于原子自洽，電子自洽，非自洽三步運算的思考。 已有10人參與
|
|||||||
|
先向各位大牛道歉， 因為上次我發(fā)了個帖子，被某牛人認為是侮辱了物理。 所以，這次要先道歉，如果我侮辱了物理，請各位嚴厲的指出。 我一定會進行嚴格的自我批評。 1、原子自洽，也就是結構優(yōu)化，其實vasp也寫了一個CHGCAR文件。但是，我們?yōu)槭裁捶且茫娮幼郧⒌贸龅腃HGCAR呢？ 小蟲認為，結構優(yōu)化的時候，每個原子步，都是基于當前的POSCAR，進行一次電子自洽，自洽完后，計算各原子所受到的力，判斷原子要不要移動移動。如此反復多次，最后一次，發(fā)現達到推出的要求了（步數達到了，或者精度達到了）。這時，我們可以讓VASP寫一個CHGCAR，而且，這個CHGCAR，就是依據最后的原子坐標‘電子自洽’出來的電荷密度。但是我們偏偏不用它，而要另外進行一次電子自洽。各位看官，這是否多余了呢？不是。因為：結構優(yōu)化后，胞發(fā)生了變化，布里淵區(qū)發(fā)生了變化，所以截斷動能對應的K，在K空間中區(qū)域不再是一個球，可能變成了橢球，也就是說，計算所用之基函數，是基于初始坐標的，對結構優(yōu)化后的坐標，這個基函數有點不合時宜。所以，我們要進行一次電子自洽，得到我們要的電荷密度。 結論：不用結構優(yōu)化的CHGCAR，不是因為這個CHGCAR沒有自洽，而是因為基函數與原子坐標不匹配得不怎么好。 2、所以，計算態(tài)密度一定要以‘電子自洽’得到的CHGCAR為輸入文件。 |
回復此樓» 收錄本帖的淘帖專輯推薦
 第一性原理 | VASP | 量子理論見解 | 科研 |
» 本帖已獲得的紅花（最新10朵）
liangqi05» 猜你喜歡
279分求調劑 一志愿211
已經有12人回復
-
297求調劑
已經有9人回復
-
085700資源與環(huán)境308求調劑
已經有13人回復
-
材料 336 求調劑
已經有4人回復
-
華東師范大學-071000生物學-293分-求調劑
已經有3人回復
-
280求調劑
已經有8人回復
-
324分 085600材料化工求調劑
已經有4人回復
-
一志愿華中科技大學，080502，354分求調劑
已經有5人回復
-
一志愿武理材料305分求調劑
已經有7人回復
-
271材料工程求調劑
已經有8人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- VASP靜態(tài)自洽計算，能量升高7eV，然后局域磁矩消失，請問是什么問題？
已經有8人回復
- 關于贗勢與自洽計算
已經有15人回復
- vasp自洽計算K點不匹配問題，急��！
已經有15人回復
- 靜態(tài)自洽計算遇到的問題
已經有5人回復
- 自洽場的fortran程序
已經有11人回復
- 靜態(tài)自洽計算能量時EDIFFG的問題
已經有14人回復
- 【求助】我在用ATK做自洽計算的時候，出現如下的錯誤：
已經有3人回復
- 【求助】VASP中，“結構優(yōu)化”和“電子自洽計算”的關系？
已經有8人回復
- 【求助】結構優(yōu)化之后的靜態(tài)自洽計算
已經有12人回復
- 【討論】自洽計算與非自洽計算的含義、所做的工作
已經有10人回復
- 【求助】SCF(自洽場理論)是怎么回事
已經有3人回復
- 【原創(chuàng)】并行自洽計算后得到和節(jié)點數一樣多的波函數文件wfc
已經有4人回復
- 【求助】vasp的電子自洽和原子核優(yōu)化過程是否可以分開？
已經有4人回復
- 【求助】靜態(tài)自洽計算讀WAVECAR有用嗎？
已經有9人回復
- 【求助】請教PWSCF晶體建模和自洽計算，不收斂
已經有21人回復
- 【求助】態(tài)密度計算究竟是自洽計算，還是非自洽計算？
已經有24人回復
- 【求助】ALGaN/GaN 薛定諤泊松方程自洽求解（插AlN層）
已經有5人回復
★ ★ ★ ★ ★
fzx2008: 金幣+5, 有見解！ 2012-10-31 18:39:56
fzx2008: 金幣+5, 有見解！ 2012-10-31 18:39:56
|
自己頂自己，亂談自洽非自洽 結構自洽就是我們照合影時，優(yōu)化每個人所站的位置 靜態(tài)自洽，就是你站的位置定了。大家再來優(yōu)化一下各自的pose 好把整張照片弄的河蟹一點。在微觀世界，可以認為 ： 電荷密度 就是這張我們需要的合影 至于態(tài)密度，就是來分析下，共有多少個pose，沒個pose有多少人在擺。 V手勢3人 吐舌頭 5人 。 。 。 高難度的一字馬， 0 人 更到難度的倒立，0人 難度越高，在微觀世界可認為是能級越高。雖然有這個態(tài)的存在，但是這個態(tài)上卻沒有電子。 |
2樓2012-10-31 18:28:19
3樓2012-10-31 19:08:48
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強帖: 3
- 應助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時
- 蟲號: 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
fzx2008: 金幣+3, 謝謝指教！ 2012-11-01 10:40:52
dyc_2008: 回帖置頂 2012-11-02 13:00:06
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
fzx2008: 金幣+3, 謝謝指教！ 2012-11-01 10:40:52
dyc_2008: 回帖置頂 2012-11-02 13:00:06
|
1, 這個CHGCAR是相當有用的。 其實在做scf的時候，初始電荷密度是有一定的選擇的。如果初始電荷密度不太合理的話，scf需要很多次才能收斂。就是因為開始的時候我們不知道什么樣的電荷密度是收斂的，所以一般我們都選用原子的電荷密度（基于原子軌道）。的確，在relax原子的時候，每一次由于原子位置或晶格的變化，會導致電荷密度或波函數有一個小的偏離，但這種偏離不大。所以我們在優(yōu)化的時候還是用了上一次的CHGCAR。 驗證這一點實際很簡單，你可以做兩個計算，一個是直接優(yōu)化。 另一個是做一步優(yōu)化，然后把CONTCAR拷貝成POSCAR，接著優(yōu)化。如此反復。第二種辦法的時間絕對比第一種時間長。 2, 盡管所有的教程都是推薦把CHGCAR作為輸入。 但本人認為， 如果你不計較是否多花時間的話，可以不這么做。直接用非常稠密的k點做dos就可以。至于為什么，它的思想跟第一類似。先用少量的k點做scf，然后用非常dense的k點，這樣時間會少一些。 但我認為，以現在機器的計算能力，完全沒必要。 |
4樓2012-11-01 08:45:46
stefsui
木蟲 (小有名氣)
- 應助: 11 (小學生)
- 金幣: 2252.3
- 散金: 114
- 紅花: 1
- 帖子: 155
- 在線: 472.7小時
- 蟲號: 872985
- 注冊: 2009-10-15
- 專業(yè): 理論和計算化學
5樓2012-11-01 23:49:00
真情無價108
木蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 1594.9
- 散金: 90
- 紅花: 2
- 帖子: 296
- 在線: 179.2小時
- 蟲號: 1516328
- 注冊: 2011-11-29
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
6樓2012-11-02 11:21:56
7樓2012-11-03 09:39:38
貝芷苡
銅蟲 (小有名氣)
- 應助: 3 (幼兒園)
- 金幣: 262.1
- 散金: 1
- 帖子: 146
- 在線: 49.2小時
- 蟲號: 2016934
- 注冊: 2012-09-20
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
8樓2013-01-03 10:44:27
paopaotu326
金蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 1023
- 散金: 5
- 紅花: 5
- 帖子: 675
- 在線: 178.6小時
- 蟲號: 951482
- 注冊: 2010-02-02
- 專業(yè): 金屬材料
9樓2013-01-11 11:07:12
撥云見日no1
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 262.3
- 帖子: 38
- 在線: 69.8小時
- 蟲號: 2463781
- 注冊: 2013-05-14
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
10樓2013-05-20 10:10:41
相關版塊跳轉
我要訂閱樓主 dyc_2008 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 材料專業(yè)求調劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南昌大學，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 藥學383 求調劑 +3 | 藥學chy 2026-03-15 | 5/250 |
|
|
|
[考研] 316求調劑 +5 | 梁茜雯 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿武理材料工程348求調劑 +3 | ￣^￣゜汗 2026-03-19 | 4/200 |
|
|
|
[考研] 0817 化學工程 299分求調劑 有科研經歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 一志愿華中農業(yè)071010，總分320求調劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 材料學求調劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
|
[考研] 環(huán)境工程調劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
|
[考研] 298-一志愿中國農業(yè)大學-求調劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
|
[考研] 085600材料與化工求調劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 344求調劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 303求調劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學調劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 0703化學336分求調劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 290求調劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 080500，材料學碩302分求調劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|