 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3469 | 回復(fù): 10 | ||||
xhzha鐵桿木蟲 (正式寫手)
|
[求助]
怎么消除虛頻，細(xì)節(jié)
|
|
大家好，用siesta算聲子振動譜，出現(xiàn)虛頻，如附件圖所示，怎么消除虛頻呢，我看到很多人討論過，說提高收斂精度，還有就是用虛頻對應(yīng)的簡正坐標(biāo)加到體系的坐標(biāo)上去，重新優(yōu)化結(jié)構(gòu)，再計算聲子。我精度設(shè)到10-8，已經(jīng)很高，但結(jié)果還是不好，所以我想試試用虛頻對應(yīng)簡正坐標(biāo)加到體系坐標(biāo)上去重新優(yōu)化試試，但是這個具體操作過程我不是很了解，一個虛頻的振動模式不是體系所有的原子的振動結(jié)果么？虛頻的簡正坐標(biāo)是根據(jù)對應(yīng)gamma點的虛頻的頻率值，解附件中的常微分方程得到的么，得到的簡正坐標(biāo)又加到體系的哪個原子上呢，是x,y,z的哪一個分量上呢，很多東西不是很明白，還望大家不膩賜教，謝謝! |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 ATK和電子輸運 | 碳納米管聲子 |
» 猜你喜歡
070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有5人回復(fù)
-
274求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有3人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
【同濟軟件】軟件（085405）考研求調(diào)劑
已經(jīng)有3人回復(fù)
-
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 如何查看計算精度？
已經(jīng)有16人回復(fù)
- 虛頻 是怎么回事？
已經(jīng)有7人回復(fù)
- 聲子譜圖（速度自相關(guān)函數(shù)的傅里葉變換之類的問題）
已經(jīng)有9人回復(fù)
- 消除七個虛頻的方法
已經(jīng)有20人回復(fù)
- 碳化硅細(xì)粉注漿成型
已經(jīng)有15人回復(fù)
- 過渡態(tài)中怎樣消除不必要的虛頻
已經(jīng)有3人回復(fù)
- 哭了，各種虛頻！
已經(jīng)有23人回復(fù)
- MS中如何調(diào)整鍵長？
已經(jīng)有6人回復(fù)
- 如何改變能帶結(jié)構(gòu)圖的精度
已經(jīng)有16人回復(fù)
- 急求細(xì)節(jié)！學(xué)院的副院長用assistant dean 對不對嗎？急著封裝推薦信。。。
已經(jīng)有8人回復(fù)
- 【求助】為什么有虛頻的結(jié)構(gòu)能量比沒有虛頻的結(jié)構(gòu)能量低？
已經(jīng)有6人回復(fù)
- 【求助】關(guān)于虛頻
已經(jīng)有8人回復(fù)
- 【求助】請教CASTEP過渡態(tài)里面如何消除虛頻？
已經(jīng)有16人回復(fù)
- 【求助】siesta計算聲子
已經(jīng)有16人回復(fù)
- 【求助】再怎么消虛頻�。�
已經(jīng)有10人回復(fù)
- 【求助】虛頻消除問題--虛頻振動坐標(biāo)能找到，怎么找不到相應(yīng)的原坐標(biāo)
已經(jīng)有6人回復(fù)
- 【求助】消除虛頻問題
已經(jīng)有8人回復(fù)
yinwei_li
木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2079.1
- 散金: 1532
- 紅花: 10
- 帖子: 896
- 在線: 194.2小時
- 蟲號: 443243
- 注冊: 2007-10-27
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
【答案】應(yīng)助回帖
★ ★
xhzha(金幣+5): 2011-09-20 09:53:08
zzy870720z(金幣+2): 謝謝指點 2011-09-20 12:39:38
xhzha(金幣+5): 2011-09-20 09:53:08
zzy870720z(金幣+2): 謝謝指點 2011-09-20 12:39:38
|
有虛頻說明你的結(jié)構(gòu)動力學(xué)不穩(wěn)定，原則上按照虛頻的模式的原子振動方式移動原子（冷凍聲子方法）可以找到動力學(xué)穩(wěn)定的結(jié)構(gòu)，但是你這個結(jié)構(gòu)虛的太離譜，基本類似能帶了，無法用冷凍聲子計算。 虛頻就虛頻，為什么要消除呢？只有那些由于計算精度不夠引起的虛頻才能想辦法消除，要是真是由于結(jié)構(gòu)不穩(wěn)定引起的虛頻，如果你人為的消除虛頻是不是有改數(shù)據(jù)的嫌疑呢？ |
2樓2011-09-20 08:49:11
xhzha
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 9454.6
- 散金: 30
- 紅花: 7
- 帖子: 409
- 在線: 577.7小時
- 蟲號: 1101711
- 注冊: 2010-09-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
|
謝謝您的回復(fù)，DFT優(yōu)化是在初始模型的基礎(chǔ)上小范圍移動原子找到局域內(nèi)穩(wěn)定結(jié)構(gòu)，可能并不是體系的最穩(wěn)定結(jié)構(gòu)，所以我就想按照虛頻的原子振動方式移動原子，改變初始構(gòu)型，優(yōu)化找到動力學(xué)最穩(wěn)定的結(jié)構(gòu)，所以用這種方法消除虛頻應(yīng)該還是合理的吧。但是也就是怎樣去移動原子這點我比較模糊，虛頻是個別原子造成的么？siesta虛頻的振動方向可以可視化么，希望您能夠講述移動原子的細(xì)節(jié)。非常感謝您，我算的是納米管的周期性結(jié)構(gòu)。 |
3樓2011-09-20 09:52:38
yinwei_li
木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2079.1
- 散金: 1532
- 紅花: 10
- 帖子: 896
- 在線: 194.2小時
- 蟲號: 443243
- 注冊: 2007-10-27
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
【答案】應(yīng)助回帖
★ ★
zzy870720z(金幣+2): 鼓勵討論 2011-09-20 12:39:51
xhzha(金幣+5): 2011-09-20 16:45:32
zzy870720z(金幣+2): 鼓勵討論 2011-09-20 12:39:51
xhzha(金幣+5): 2011-09-20 16:45:32
|
通過原子移動來找穩(wěn)定結(jié)構(gòu)的方法是有一定限制的，如果你的聲子譜只是在布里淵區(qū)某個特殊點有虛頻，找到此點對應(yīng)的本征矢（本征矢可表征原子如何振動），沿著本征矢逐步移動原子并計算每移動一步體系的能量的變化，你會發(fā)現(xiàn)能量的變化會出現(xiàn)一個低谷，那么這個能量低谷的位置就對應(yīng)新的穩(wěn)定結(jié)構(gòu)。計算這個穩(wěn)定結(jié)構(gòu)的聲子譜應(yīng)該是沒有虛頻的。但是你給出的聲子譜在整個布里淵區(qū)都有虛頻，用這種方法可能就不適應(yīng)了，因為你沒辦法同時按照所有布里淵區(qū)點的本征矢移動原子。 或許大家有別的方法，這個我就不太清楚了！ |
4樓2011-09-20 11:41:18
yinwei_li
木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2079.1
- 散金: 1532
- 紅花: 10
- 帖子: 896
- 在線: 194.2小時
- 蟲號: 443243
- 注冊: 2007-10-27
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
5樓2011-09-20 11:43:53
xhzha
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 9454.6
- 散金: 30
- 紅花: 7
- 帖子: 409
- 在線: 577.7小時
- 蟲號: 1101711
- 注冊: 2010-09-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
6樓2011-09-20 16:43:13
yinwei_li
木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2079.1
- 散金: 1532
- 紅花: 10
- 帖子: 896
- 在線: 194.2小時
- 蟲號: 443243
- 注冊: 2007-10-27
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
7樓2011-09-20 19:55:34
xhzha
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 9454.6
- 散金: 30
- 紅花: 7
- 帖子: 409
- 在線: 577.7小時
- 蟲號: 1101711
- 注冊: 2010-09-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
8樓2011-09-21 08:07:27
yinwei_li
木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2079.1
- 散金: 1532
- 紅花: 10
- 帖子: 896
- 在線: 194.2小時
- 蟲號: 443243
- 注冊: 2007-10-27
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
9樓2011-09-21 09:47:49
xhzha
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 9454.6
- 散金: 30
- 紅花: 7
- 帖子: 409
- 在線: 577.7小時
- 蟲號: 1101711
- 注冊: 2010-09-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
10樓2011-09-21 16:03:19
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xhzha 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 12/600 |
|
|
[考研] 344求調(diào)劑 +5 | knight344 2026-03-16 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|