| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1360 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
[交流]
vasp晶格優(yōu)化時(shí)的對(duì)稱(chēng)性問(wèn)題 已有1人參與
|
|||
|
在VASP進(jìn)行晶格優(yōu)化的時(shí)候考慮到體系的對(duì)稱(chēng)性,這一點(diǎn)已經(jīng)過(guò)驗(yàn)證。在Aantase體系中,使用LDA+U方法,若將原子的位置取到晶格矢量的恰好1/2,1/4,3/4的時(shí)候,為了滿(mǎn)足晶格的對(duì)稱(chēng)性要求,會(huì)使得晶格體積膨脹很大(ISIF=3),誤差嚴(yán)重。這個(gè)問(wèn)題詳見(jiàn)http://www.gaoyang168.com/bbs/viewthread.php?tid=3152699&fpage=1。 后來(lái)我在每一個(gè)原子坐標(biāo)之后加入一個(gè)變化(正負(fù)0.0000001),則優(yōu)化的時(shí)候?qū)ΨQ(chēng)性不會(huì)被考慮(至少不會(huì)考慮得那么多),但是計(jì)算的速度慢了很多,因?yàn)榭紤]對(duì)稱(chēng)性的時(shí)候只有少數(shù)幾個(gè)坐標(biāo)會(huì)發(fā)生移動(dòng)。 但這是也有一個(gè)小小的疑問(wèn)。優(yōu)化的時(shí)候晶格矢量就不再是正交的了,也就是在CONTCAR中的晶格矢量的那三行中,原來(lái)為零的非對(duì)角元素有了一定的增量。這樣計(jì)算出來(lái)的結(jié)果是可信的么?大家在優(yōu)化的時(shí)候都考慮了對(duì)稱(chēng)性么?在POSCAR中輸入的原子的坐標(biāo)一般是如何取得的呢,是直接根據(jù)幾何位置將坐標(biāo)算出來(lái)么?對(duì)于1/n這樣的高對(duì)稱(chēng)點(diǎn)又應(yīng)該如何處理呢(因?yàn)榫拖裆厦嫠f(shuō)的,在某些時(shí)候,考慮了對(duì)稱(chēng)性會(huì)對(duì)晶格常數(shù)有很大影響)? |

金蟲(chóng) (正式寫(xiě)手)
|
我覺(jué)得樓主你多慮了。 你優(yōu)化的晶格常數(shù)和文獻(xiàn)差不多的話(huà)體積也應(yīng)該差不多的。 你優(yōu)化的時(shí)候晶格矢量就不再是正交的了,這點(diǎn)我覺(jué)得不可信。我只有一次把POSCAR原子順序顛倒了,得到的的基矢才變得不正交了(優(yōu)化Rutial SnO2)。 我優(yōu)化的時(shí)候都沒(méi)有考慮對(duì)稱(chēng)性,在POSCAR中輸入的原子的坐標(biāo)是用Material Studio (很有用)導(dǎo)出來(lái)的。 我覺(jué)得你讓VASP自己優(yōu)化就好了,沒(méi)必要加那么多限制(你在每一個(gè)原子坐標(biāo)之后加入一個(gè)變化(正負(fù)0.0000001))。如果一個(gè)原子是在高對(duì)稱(chēng)位置上,那么弛豫的時(shí)候它的位置一般都不會(huì)變化的。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑 +4 | 亞克西good 2026-03-26 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料與化工328分調(diào)劑 +7 | 。,。,。,。i 2026-03-23 | 7/350 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 277跪求調(diào)劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|
[考研] 343求調(diào)劑 +6 | 愛(ài)羈絆 2026-03-29 | 6/300 |
|
|
[考研] 348求調(diào)劑 +5 | 小懶蟲(chóng)不懶了 2026-03-28 | 5/250 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 壓國(guó)家一區(qū)線(xiàn),求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見(jiàn)α 2026-03-25 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 材料/農(nóng)業(yè)專(zhuān)業(yè),07/08開(kāi)頭均可,過(guò)線(xiàn)就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|