 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3277 | 回復(fù): 12 | |||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||||||
[交流]
關(guān)于原子自洽，電子自洽，非自洽三步運算的思考。 已有10人參與
|
|||||||
|
先向各位大牛道歉， 因為上次我發(fā)了個帖子，被某牛人認(rèn)為是侮辱了物理。 所以，這次要先道歉，如果我侮辱了物理，請各位嚴(yán)厲的指出。 我一定會進行嚴(yán)格的自我批評。 1、原子自洽，也就是結(jié)構(gòu)優(yōu)化，其實vasp也寫了一個CHGCAR文件。但是，我們?yōu)槭裁捶且�，電子自洽得出的CHGCAR呢？ 小蟲認(rèn)為，結(jié)構(gòu)優(yōu)化的時候，每個原子步，都是基于當(dāng)前的POSCAR，進行一次電子自洽，自洽完后，計算各原子所受到的力，判斷原子要不要移動移動。如此反復(fù)多次，最后一次，發(fā)現(xiàn)達到推出的要求了（步數(shù)達到了，或者精度達到了）。這時，我們可以讓VASP寫一個CHGCAR，而且，這個CHGCAR，就是依據(jù)最后的原子坐標(biāo)‘電子自洽’出來的電荷密度。但是我們偏偏不用它，而要另外進行一次電子自洽。各位看官，這是否多余了呢？不是。因為：結(jié)構(gòu)優(yōu)化后，胞發(fā)生了變化，布里淵區(qū)發(fā)生了變化，所以截斷動能對應(yīng)的K，在K空間中區(qū)域不再是一個球，可能變成了橢球，也就是說，計算所用之基函數(shù)，是基于初始坐標(biāo)的，對結(jié)構(gòu)優(yōu)化后的坐標(biāo)，這個基函數(shù)有點不合時宜。所以，我們要進行一次電子自洽，得到我們要的電荷密度。 結(jié)論：不用結(jié)構(gòu)優(yōu)化的CHGCAR，不是因為這個CHGCAR沒有自洽，而是因為基函數(shù)與原子坐標(biāo)不匹配得不怎么好。 2、所以，計算態(tài)密度一定要以‘電子自洽’得到的CHGCAR為輸入文件。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性原理 | VASP | 量子理論見解 | 科研 |
» 本帖已獲得的紅花（最新10朵）
liangqi05» 猜你喜歡
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有23人回復(fù)
-
08工科，295，接受跨專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有12人回復(fù)
-
334分 一志愿武理 材料求調(diào)劑
已經(jīng)有16人回復(fù)
-
0703化學(xué)321分求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料專碩 085600求調(diào)劑
已經(jīng)有7人回復(fù)
-
297求調(diào)劑
已經(jīng)有15人回復(fù)
-
11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè)
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有9人回復(fù)
-
化工專碩求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- VASP靜態(tài)自洽計算，能量升高7eV，然后局域磁矩消失，請問是什么問題？
已經(jīng)有8人回復(fù)
- 關(guān)于贗勢與自洽計算
已經(jīng)有15人回復(fù)
- vasp自洽計算K點不匹配問題，急！！
已經(jīng)有15人回復(fù)
- 靜態(tài)自洽計算遇到的問題
已經(jīng)有5人回復(fù)
- 自洽場的fortran程序
已經(jīng)有11人回復(fù)
- 靜態(tài)自洽計算能量時EDIFFG的問題
已經(jīng)有14人回復(fù)
- 【求助】我在用ATK做自洽計算的時候，出現(xiàn)如下的錯誤：
已經(jīng)有3人回復(fù)
- 【求助】VASP中，“結(jié)構(gòu)優(yōu)化”和“電子自洽計算”的關(guān)系？
已經(jīng)有8人回復(fù)
- 【求助】結(jié)構(gòu)優(yōu)化之后的靜態(tài)自洽計算
已經(jīng)有12人回復(fù)
- 【討論】自洽計算與非自洽計算的含義、所做的工作
已經(jīng)有10人回復(fù)
- 【求助】SCF(自洽場理論)是怎么回事
已經(jīng)有3人回復(fù)
- 【原創(chuàng)】并行自洽計算后得到和節(jié)點數(shù)一樣多的波函數(shù)文件wfc
已經(jīng)有4人回復(fù)
- 【求助】vasp的電子自洽和原子核優(yōu)化過程是否可以分開？
已經(jīng)有4人回復(fù)
- 【求助】靜態(tài)自洽計算讀WAVECAR有用嗎？
已經(jīng)有9人回復(fù)
- 【求助】請教PWSCF晶體建模和自洽計算，不收斂
已經(jīng)有21人回復(fù)
- 【求助】態(tài)密度計算究竟是自洽計算，還是非自洽計算？
已經(jīng)有24人回復(fù)
- 【求助】ALGaN/GaN 薛定諤泊松方程自洽求解（插AlN層）
已經(jīng)有5人回復(fù)
paopaotu326
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1023
- 散金: 5
- 紅花: 5
- 帖子: 675
- 在線: 178.6小時
- 蟲號: 951482
- 注冊: 2010-02-02
- 專業(yè): 金屬材料
9樓2013-01-11 11:07:12
★ ★ ★ ★ ★
fzx2008: 金幣+5, 有見解！ 2012-10-31 18:39:56
fzx2008: 金幣+5, 有見解！ 2012-10-31 18:39:56
|
自己頂自己，亂談自洽非自洽 結(jié)構(gòu)自洽就是我們照合影時，優(yōu)化每個人所站的位置 靜態(tài)自洽，就是你站的位置定了。大家再來優(yōu)化一下各自的pose 好把整張照片弄的河蟹一點。在微觀世界，可以認(rèn)為 ： 電荷密度 就是這張我們需要的合影 至于態(tài)密度，就是來分析下，共有多少個pose，沒個pose有多少人在擺。 V手勢3人 吐舌頭 5人 。 。 。 高難度的一字馬， 0 人 更到難度的倒立，0人 難度越高，在微觀世界可認(rèn)為是能級越高。雖然有這個態(tài)的存在，但是這個態(tài)上卻沒有電子。 |
2樓2012-10-31 18:28:19
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時
- 蟲號: 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
★ ★ ★ ★
小木蟲: 金幣+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是相當(dāng)有用的。 其實在做scf的時候，初始電荷密度是有一定的選擇的。如果初始電荷密度不太合理的話，scf需要很多次才能收斂。就是因為開始的時候我們不知道什么樣的電荷密度是收斂的，所以一般我們都選用原子的電荷密度（基于原子軌道）。的確，在relax原子的時候，每一次由于原子位置或晶格的變化，會導(dǎo)致電荷密度或波函數(shù)有一個小的偏離，但這種偏離不大。所以我們在優(yōu)化的時候還是用了上一次的CHGCAR。 驗證這一點實際很簡單，你可以做兩個計算，一個是直接優(yōu)化。 另一個是做一步優(yōu)化，然后把CONTCAR拷貝成POSCAR，接著優(yōu)化。如此反復(fù)。第二種辦法的時間絕對比第一種時間長。 2, 盡管所有的教程都是推薦把CHGCAR作為輸入。 但本人認(rèn)為， 如果你不計較是否多花時間的話，可以不這么做。直接用非常稠密的k點做dos就可以。至于為什么，它的思想跟第一類似。先用少量的k點做scf，然后用非常dense的k點，這樣時間會少一些。 但我認(rèn)為，以現(xiàn)在機器的計算能力，完全沒必要。 |
4樓2012-11-01 08:45:46
stefsui
木蟲 (小有名氣)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 2252.3
- 散金: 114
- 紅花: 1
- 帖子: 155
- 在線: 472.7小時
- 蟲號: 872985
- 注冊: 2009-10-15
- 專業(yè): 理論和計算化學(xué)
5樓2012-11-01 23:49:00
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)碩333求調(diào)劑 +12 | 北道巷 2026-03-24 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 295材料工程專碩求調(diào)劑 +8 | 1428151015 2026-03-27 | 8/400 |
|
|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-30 | 5/250 |
|
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +5 | 藍月亮亮 2026-03-30 | 5/250 |
|
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-23 | 8/400 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +8 | hanamiko 2026-03-27 | 8/400 |
|
|
[考研] 材料工程專碩求調(diào)劑 +6 | hyl3153942 2026-03-29 | 6/300 |
|
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我！��！ +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|