 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3273 | 回復(fù): 12 | |||||||
[交流]
關(guān)于原子自洽，電子自洽，非自洽三步運(yùn)算的思考。 已有10人參與
|
|||||||
|
先向各位大牛道歉， 因?yàn)樯洗挝野l(fā)了個(gè)帖子，被某牛人認(rèn)為是侮辱了物理。 所以，這次要先道歉，如果我侮辱了物理，請(qǐng)各位嚴(yán)厲的指出。 我一定會(huì)進(jìn)行嚴(yán)格的自我批評(píng)。 1、原子自洽，也就是結(jié)構(gòu)優(yōu)化，其實(shí)vasp也寫了一個(gè)CHGCAR文件。但是，我們?yōu)槭裁捶且�，電子自洽得出的CHGCAR呢？ 小蟲(chóng)認(rèn)為，結(jié)構(gòu)優(yōu)化的時(shí)候，每個(gè)原子步，都是基于當(dāng)前的POSCAR，進(jìn)行一次電子自洽，自洽完后，計(jì)算各原子所受到的力，判斷原子要不要移動(dòng)移動(dòng)。如此反復(fù)多次，最后一次，發(fā)現(xiàn)達(dá)到推出的要求了（步數(shù)達(dá)到了，或者精度達(dá)到了）。這時(shí)，我們可以讓VASP寫一個(gè)CHGCAR，而且，這個(gè)CHGCAR，就是依據(jù)最后的原子坐標(biāo)‘電子自洽’出來(lái)的電荷密度。但是我們偏偏不用它，而要另外進(jìn)行一次電子自洽。各位看官，這是否多余了呢？不是。因?yàn)椋航Y(jié)構(gòu)優(yōu)化后，胞發(fā)生了變化，布里淵區(qū)發(fā)生了變化，所以截?cái)鄤?dòng)能對(duì)應(yīng)的K，在K空間中區(qū)域不再是一個(gè)球，可能變成了橢球，也就是說(shuō)，計(jì)算所用之基函數(shù)，是基于初始坐標(biāo)的，對(duì)結(jié)構(gòu)優(yōu)化后的坐標(biāo)，這個(gè)基函數(shù)有點(diǎn)不合時(shí)宜。所以，我們要進(jìn)行一次電子自洽，得到我們要的電荷密度。 結(jié)論：不用結(jié)構(gòu)優(yōu)化的CHGCAR，不是因?yàn)檫@個(gè)CHGCAR沒(méi)有自洽，而是因?yàn)榛瘮?shù)與原子坐標(biāo)不匹配得不怎么好。 2、所以，計(jì)算態(tài)密度一定要以‘電子自洽’得到的CHGCAR為輸入文件。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性原理 | VASP | 量子理論見(jiàn)解 | 科研 |
» 本帖已獲得的紅花（最新10朵）
liangqi05» 猜你喜歡
材料與化工考研調(diào)劑
已經(jīng)有5人回復(fù)
-
299求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料292調(diào)劑
已經(jīng)有5人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有4人回復(fù)
-
收08調(diào)劑生
已經(jīng)有3人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有9人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
086003食品工程求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856材料專碩353求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- VASP靜態(tài)自洽計(jì)算，能量升高7eV，然后局域磁矩消失，請(qǐng)問(wèn)是什么問(wèn)題？
已經(jīng)有8人回復(fù)
- 關(guān)于贗勢(shì)與自洽計(jì)算
已經(jīng)有15人回復(fù)
- vasp自洽計(jì)算K點(diǎn)不匹配問(wèn)題，急��！
已經(jīng)有15人回復(fù)
- 靜態(tài)自洽計(jì)算遇到的問(wèn)題
已經(jīng)有5人回復(fù)
- 自洽場(chǎng)的fortran程序
已經(jīng)有11人回復(fù)
- 靜態(tài)自洽計(jì)算能量時(shí)EDIFFG的問(wèn)題
已經(jīng)有14人回復(fù)
- 【求助】我在用ATK做自洽計(jì)算的時(shí)候，出現(xiàn)如下的錯(cuò)誤：
已經(jīng)有3人回復(fù)
- 【求助】VASP中，“結(jié)構(gòu)優(yōu)化”和“電子自洽計(jì)算”的關(guān)系？
已經(jīng)有8人回復(fù)
- 【求助】結(jié)構(gòu)優(yōu)化之后的靜態(tài)自洽計(jì)算
已經(jīng)有12人回復(fù)
- 【討論】自洽計(jì)算與非自洽計(jì)算的含義、所做的工作
已經(jīng)有10人回復(fù)
- 【求助】SCF(自洽場(chǎng)理論)是怎么回事
已經(jīng)有3人回復(fù)
- 【原創(chuàng)】并行自洽計(jì)算后得到和節(jié)點(diǎn)數(shù)一樣多的波函數(shù)文件wfc
已經(jīng)有4人回復(fù)
- 【求助】vasp的電子自洽和原子核優(yōu)化過(guò)程是否可以分開(kāi)？
已經(jīng)有4人回復(fù)
- 【求助】靜態(tài)自洽計(jì)算讀WAVECAR有用嗎？
已經(jīng)有9人回復(fù)
- 【求助】請(qǐng)教PWSCF晶體建模和自洽計(jì)算，不收斂
已經(jīng)有21人回復(fù)
- 【求助】態(tài)密度計(jì)算究竟是自洽計(jì)算，還是非自洽計(jì)算？
已經(jīng)有24人回復(fù)
- 【求助】ALGaN/GaN 薛定諤泊松方程自洽求解（插AlN層）
已經(jīng)有5人回復(fù)
gavinliu7390
木蟲(chóng) (著名寫手)
葉落鷹飛
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時(shí)
- 蟲(chóng)號(hào): 675997
- 注冊(cè): 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
★ ★ ★ ★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
fzx2008: 金幣+3, 謝謝指教！ 2012-11-01 10:40:52
dyc_2008: 回帖置頂 2012-11-02 13:00:06
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
fzx2008: 金幣+3, 謝謝指教！ 2012-11-01 10:40:52
dyc_2008: 回帖置頂 2012-11-02 13:00:06
|
1, 這個(gè)CHGCAR是相當(dāng)有用的。 其實(shí)在做scf的時(shí)候，初始電荷密度是有一定的選擇的。如果初始電荷密度不太合理的話，scf需要很多次才能收斂。就是因?yàn)殚_(kāi)始的時(shí)候我們不知道什么樣的電荷密度是收斂的，所以一般我們都選用原子的電荷密度（基于原子軌道）。的確，在relax原子的時(shí)候，每一次由于原子位置或晶格的變化，會(huì)導(dǎo)致電荷密度或波函數(shù)有一個(gè)小的偏離，但這種偏離不大。所以我們?cè)趦?yōu)化的時(shí)候還是用了上一次的CHGCAR。 驗(yàn)證這一點(diǎn)實(shí)際很簡(jiǎn)單，你可以做兩個(gè)計(jì)算，一個(gè)是直接優(yōu)化。 另一個(gè)是做一步優(yōu)化，然后把CONTCAR拷貝成POSCAR，接著優(yōu)化。如此反復(fù)。第二種辦法的時(shí)間絕對(duì)比第一種時(shí)間長(zhǎng)。 2, 盡管所有的教程都是推薦把CHGCAR作為輸入。 但本人認(rèn)為， 如果你不計(jì)較是否多花時(shí)間的話，可以不這么做。直接用非常稠密的k點(diǎn)做dos就可以。至于為什么，它的思想跟第一類似。先用少量的k點(diǎn)做scf，然后用非常dense的k點(diǎn)，這樣時(shí)間會(huì)少一些。 但我認(rèn)為，以現(xiàn)在機(jī)器的計(jì)算能力，完全沒(méi)必要。 |
4樓2012-11-01 08:45:46
★ ★ ★ ★ ★
fzx2008: 金幣+5, 有見(jiàn)解！ 2012-10-31 18:39:56
fzx2008: 金幣+5, 有見(jiàn)解！ 2012-10-31 18:39:56
|
自己頂自己，亂談自洽非自洽 結(jié)構(gòu)自洽就是我們照合影時(shí)，優(yōu)化每個(gè)人所站的位置 靜態(tài)自洽，就是你站的位置定了。大家再來(lái)優(yōu)化一下各自的pose 好把整張照片弄的河蟹一點(diǎn)。在微觀世界，可以認(rèn)為 ： 電荷密度 就是這張我們需要的合影 至于態(tài)密度，就是來(lái)分析下，共有多少個(gè)pose，沒(méi)個(gè)pose有多少人在擺。 V手勢(shì)3人 吐舌頭 5人 。 。 。 高難度的一字馬， 0 人 更到難度的倒立，0人 難度越高，在微觀世界可認(rèn)為是能級(jí)越高。雖然有這個(gè)態(tài)的存在，但是這個(gè)態(tài)上卻沒(méi)有電子。 |
2樓2012-10-31 18:28:19
3樓2012-10-31 19:08:48
stefsui
木蟲(chóng) (小有名氣)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 2252.3
- 散金: 114
- 紅花: 1
- 帖子: 155
- 在線: 472.7小時(shí)
- 蟲(chóng)號(hào): 872985
- 注冊(cè): 2009-10-15
- 專業(yè): 理論和計(jì)算化學(xué)
5樓2012-11-01 23:49:00
真情無(wú)價(jià)108
木蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1594.9
- 散金: 90
- 紅花: 2
- 帖子: 296
- 在線: 179.2小時(shí)
- 蟲(chóng)號(hào): 1516328
- 注冊(cè): 2011-11-29
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
6樓2012-11-02 11:21:56
7樓2012-11-03 09:39:38
貝芷苡
銅蟲(chóng) (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 262.1
- 散金: 1
- 帖子: 146
- 在線: 49.2小時(shí)
- 蟲(chóng)號(hào): 2016934
- 注冊(cè): 2012-09-20
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
8樓2013-01-03 10:44:27
paopaotu326
金蟲(chóng) (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1023
- 散金: 5
- 紅花: 5
- 帖子: 675
- 在線: 178.6小時(shí)
- 蟲(chóng)號(hào): 951482
- 注冊(cè): 2010-02-02
- 專業(yè): 金屬材料
9樓2013-01-11 11:07:12
撥云見(jiàn)日no1
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 262.3
- 帖子: 38
- 在線: 69.8小時(shí)
- 蟲(chóng)號(hào): 2463781
- 注冊(cè): 2013-05-14
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
10樓2013-05-20 10:10:41
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 dyc_2008 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料292調(diào)劑 +5 | 橘頌思美人 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑 +5 | 研研，接電話 2026-03-24 | 6/300 |
|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
|
[考研]
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22 | 5/250 |
|
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|