 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3280 | 回復(fù): 12 | |||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||||
[交流]
關(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呢？ 小蟲認(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)‘電子自洽’出來的電荷密度。但是我們偏偏不用它，而要另外進(jìn)行一次電子自洽。各位看官，這是否多余了呢？不是。因?yàn)椋航Y(jié)構(gòu)優(yōu)化后，胞發(fā)生了變化，布里淵區(qū)發(fā)生了變化，所以截?cái)鄤?dòng)能對(duì)應(yīng)的K，在K空間中區(qū)域不再是一個(gè)球，可能變成了橢球，也就是說，計(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沒有自洽，而是因?yàn)榛瘮?shù)與原子坐標(biāo)不匹配得不怎么好。 2、所以，計(jì)算態(tài)密度一定要以‘電子自洽’得到的CHGCAR為輸入文件。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性原理 | VASP | 量子理論見解 | 科研 |
» 本帖已獲得的紅花（最新10朵）
liangqi05» 猜你喜歡
材料考研調(diào)劑
已經(jīng)有7人回復(fù)
-
270調(diào)劑
已經(jīng)有7人回復(fù)
-
275學(xué)碩081000服從調(diào)劑到其他專業(yè)，保不住本專業(yè)了
已經(jīng)有6人回復(fù)
-
0710生物學(xué)，325求調(diào)劑
已經(jīng)有3人回復(fù)
-
296材料專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
385分 生物學(xué)（071000）求調(diào)劑
已經(jīng)有5人回復(fù)
-
377求調(diào)劑
已經(jīng)有3人回復(fù)
-
【求調(diào)劑】新能源材料本科，一志愿211，初試321
已經(jīng)有6人回復(fù)
-
316求調(diào)劑
已經(jīng)有11人回復(fù)
-
321求調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- VASP靜態(tài)自洽計(jì)算，能量升高7eV，然后局域磁矩消失，請(qǐng)問是什么問題？
已經(jīng)有8人回復(fù)
- 關(guān)于贗勢與自洽計(jì)算
已經(jīng)有15人回復(fù)
- vasp自洽計(jì)算K點(diǎn)不匹配問題，急�。�
已經(jīng)有15人回復(fù)
- 靜態(tài)自洽計(jì)算遇到的問題
已經(jīng)有5人回復(fù)
- 自洽場的fortran程序
已經(jīng)有11人回復(fù)
- 靜態(tài)自洽計(jì)算能量時(shí)EDIFFG的問題
已經(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(自洽場理論)是怎么回事
已經(jīng)有3人回復(fù)
- 【原創(chuàng)】并行自洽計(jì)算后得到和節(jié)點(diǎn)數(shù)一樣多的波函數(shù)文件wfc
已經(jīng)有4人回復(fù)
- 【求助】vasp的電子自洽和原子核優(yōu)化過程是否可以分開？
已經(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
木蟲 (著名寫手)
葉落鷹飛
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時(shí)
- 蟲號(hào): 675997
- 注冊(cè): 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
★ ★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
fzx2008: 金幣+3, 謝謝指教！ 2012-11-01 10:40:52
dyc_2008: 回帖置頂 2012-11-02 13:00:06
小木蟲: 金幣+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)殚_始的時(shí)候我們不知道什么樣的電荷密度是收斂的，所以一般我們都選用原子的電荷密度（基于原子軌道）。的確，在relax原子的時(shí)候，每一次由于原子位置或晶格的變化，會(huì)導(dǎo)致電荷密度或波函數(shù)有一個(gè)小的偏離，但這種偏離不大。所以我們?cè)趦?yōu)化的時(shí)候還是用了上一次的CHGCAR。 驗(yàn)證這一點(diǎn)實(shí)際很簡單，你可以做兩個(gè)計(jì)算，一個(gè)是直接優(yōu)化。 另一個(gè)是做一步優(yōu)化，然后把CONTCAR拷貝成POSCAR，接著優(yōu)化。如此反復(fù)。第二種辦法的時(shí)間絕對(duì)比第一種時(shí)間長。 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ì)算能力，完全沒必要。 |
4樓2012-11-01 08:45:46
★ ★ ★ ★ ★
fzx2008: 金幣+5, 有見解！ 2012-10-31 18:39:56
fzx2008: 金幣+5, 有見解！ 2012-10-31 18:39:56
|
自己頂自己，亂談自洽非自洽 結(jié)構(gòu)自洽就是我們照合影時(shí)，優(yōu)化每個(gè)人所站的位置 靜態(tài)自洽，就是你站的位置定了。大家再來優(yōu)化一下各自的pose 好把整張照片弄的河蟹一點(diǎn)。在微觀世界，可以認(rèn)為 ： 電荷密度 就是這張我們需要的合影 至于態(tài)密度，就是來分析下，共有多少個(gè)pose，沒個(gè)pose有多少人在擺。 V手勢3人 吐舌頭 5人 。 。 。 高難度的一字馬， 0 人 更到難度的倒立，0人 難度越高，在微觀世界可認(rèn)為是能級(jí)越高。雖然有這個(gè)態(tài)的存在，但是這個(gè)態(tài)上卻沒有電子。 |
2樓2012-10-31 18:28:19
stefsui
木蟲 (小有名氣)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 2252.3
- 散金: 114
- 紅花: 1
- 帖子: 155
- 在線: 472.7小時(shí)
- 蟲號(hào): 872985
- 注冊(cè): 2009-10-15
- 專業(yè): 理論和計(jì)算化學(xué)
5樓2012-11-01 23:49:00
真情無價(jià)108
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1594.9
- 散金: 90
- 紅花: 2
- 帖子: 296
- 在線: 179.2小時(shí)
- 蟲號(hào): 1516328
- 注冊(cè): 2011-11-29
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
6樓2012-11-02 11:21:56
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【求調(diào)劑】新能源材料本科，一志愿211，初試321 +6 | 求調(diào)劑學(xué)校， 2026-04-02 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 282求調(diào)劑 +7 | 呼吸都是減肥 2026-04-01 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +10 | 陽陽哇塞 2026-04-01 | 10/500 |
|
|
|
[考研] 302求調(diào)劑一志愿北航070300，本科鄭大化學(xué) +8 | 圣日耳曼條 2026-04-01 | 11/550 |
|
|
|
[考研] 270求調(diào)劑 +7 | 小杰pp 2026-03-31 | 8/400 |
|
|
|
[考研] 286求調(diào)劑 +16 | PolarBear11 2026-03-26 | 16/800 |
|
|
|
[考研] 0703一志愿南師大334求調(diào)劑 +4 | seven7yu 2026-03-30 | 4/200 |
|
|
|
[考研] 330分求調(diào)劑 +11 | qzenlc 2026-03-29 | 11/550 |
|
|
|
[考研] 309分085801求調(diào)劑 +7 | 學(xué)員Gtwj7W 2026-03-31 | 7/350 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-29 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +8 | 嗯嗯嗯嗯2 2026-03-27 | 8/400 |
|
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
|
[考研] 求化學(xué)調(diào)劑 +12 | wulanna 2026-03-28 | 12/600 |
|
|
|
[考研] 354求調(diào)劑 +3 | lxb598 2026-03-31 | 4/200 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
|
[考研] 085602化工求調(diào)劑（331分） +8 | 111@127 2026-03-30 | 8/400 |
|
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|