 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2776 | 回復: 11 | ||
孫遠慧木蟲 (小有名氣)
|
[求助]
VASP優(yōu)化幾何時，ISIF=3時，如何只改變晶格，而不改變位矢方向. 已有2人參與
|
|
|
各位大神，我在ms建了一個立方晶系，之后用ISIF=2優(yōu)化幾何，發(fā)現(xiàn)OUTCAR里面的FORCE on cell下面的YY的值并不是負值，而是正的11點多（圖一）。之后試了ISIF=3優(yōu)化幾何發(fā)現(xiàn)YY降到負的了。但是伴隨產(chǎn)生的另一個新問題就是CONTCAR的坐標發(fā)生變化，如圖（圖二、三）所示。有沒有人知道設置什么參數(shù)，可以只改變a,b,c的長度，而不改變各個位矢的夾角，就是原來是零的位置優(yōu)化完之后還是零，但是可以改變晶胞的三個邊長。希望各位大神指教��！ 1.jpg 2.jpg 3.jpg |
回復此樓» 猜你喜歡
招收調(diào)劑碩士
已經(jīng)有12人回復
-
296求調(diào)劑
已經(jīng)有8人回復
-
招收博士1-2人
已經(jīng)有3人回復
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有11人回復
-
279分求調(diào)劑 一志愿211
已經(jīng)有9人回復
-
274求調(diào)劑
已經(jīng)有8人回復
-
290求調(diào)劑
已經(jīng)有4人回復
-
一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑
已經(jīng)有4人回復
-
295復試調(diào)劑
已經(jīng)有6人回復
-
工科材料085601 279求調(diào)劑
已經(jīng)有8人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求教各位大俠，優(yōu)化石墨烯的晶格常數(shù)時，IBRION與ISIF是如何設置的呢？
已經(jīng)有4人回復
- vasp 優(yōu)化大概的模型 isif=?合適啊
已經(jīng)有18人回復
- VASP結(jié)構(gòu)優(yōu)化ISIF=3，為什么晶格常數(shù)不變
已經(jīng)有12人回復
- 在vasp計算里，有沒有好的只優(yōu)化原子位置和晶胞體積的方法？
已經(jīng)有8人回復
- 求助成功：VASP進行體積弛豫時用ISIF＝7，如何看弛豫后的結(jié)果
已經(jīng)有7人回復
- 【求助】vasp優(yōu)化時選擇了isif=3 怎么限制晶格的一個維度呢？
已經(jīng)有16人回復
孫遠慧
木蟲 (小有名氣)
- 應助: 20 (小學生)
- 金幣: 1450.3
- 散金: 282
- 紅花: 5
- 帖子: 229
- 在線: 224.3小時
- 蟲號: 1375586
- 注冊: 2011-08-22
- 性別: GG
- 專業(yè): 理論和計算化學
9樓2014-01-10 16:30:02
kjsun
木蟲 (正式寫手)
- 應助: 34 (小學生)
- 金幣: 5557.5
- 散金: 100
- 紅花: 5
- 帖子: 434
- 在線: 469.8小時
- 蟲號: 285070
- 注冊: 2006-10-14
- 性別: GG
- 專業(yè): 催化化學
★ ★
sunyang1988: 金幣+2, 謝謝交流 2014-01-17 13:52:49
sunyang1988: 金幣+2, 謝謝交流 2014-01-17 13:52:49
|
vasp 書上的例子稍微變通一下就能滿足你的要求。 對比vasp的例子，仔細看下面的變化 rm WAVECAR for i in 3.7 3.8 3.9 4.0 4.1 do cat >POSCAR <<! fcc: 1 $i 0.0 0.0 0.0 3.7 0.0 0.0 0.0 3.7 1 cartesian 0 0 0 ! echo "a= $i" ; vamp E=`tail -1 OSZICAR` ; echo $i $E >>SUMMARY.fcc done cat SUMMARY.fcc |
11樓2014-01-17 09:16:31
2樓2014-01-09 11:24:36
孫遠慧
木蟲 (小有名氣)
- 應助: 20 (小學生)
- 金幣: 1450.3
- 散金: 282
- 紅花: 5
- 帖子: 229
- 在線: 224.3小時
- 蟲號: 1375586
- 注冊: 2011-08-22
- 性別: GG
- 專業(yè): 理論和計算化學
|
現(xiàn)在呢？ 第一張圖就是FORCE on cell 里面YY豎行和 in total橫行對應的值是11點多，第二張到第三張的是 —————————————————————————— 20.000000000 0.000000000 0.000000000 0.000000000 45.000000000 0.000000000 0.000000000 0.000000000 7.571900000 —————————————————————————— 變成圖三 —————————————————————————— 19.1459010205 -0.007219492 0.0378856652 -0.0141628126 45.608278082 0.0187786055 0.01463468126 0.031684302 7.5643364137 —————————————————————————— 我想讓圖二中20 45 7.5719這三個數(shù)據(jù)變化，而原來是零的不變化，該添加什么參數(shù)��？ 2014-01-09_11-59-19.jpg |
3樓2014-01-09 12:09:31
孫遠慧
木蟲 (小有名氣)
- 應助: 20 (小學生)
- 金幣: 1450.3
- 散金: 282
- 紅花: 5
- 帖子: 229
- 在線: 224.3小時
- 蟲號: 1375586
- 注冊: 2011-08-22
- 性別: GG
- 專業(yè): 理論和計算化學
4樓2014-01-09 12:11:48
kjsun
木蟲 (正式寫手)
- 應助: 34 (小學生)
- 金幣: 5557.5
- 散金: 100
- 紅花: 5
- 帖子: 434
- 在線: 469.8小時
- 蟲號: 285070
- 注冊: 2006-10-14
- 性別: GG
- 專業(yè): 催化化學
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
孫遠慧: 金幣+5, ★★★很有幫助 2014-01-10 09:28:20
sunyang1988: 金幣+2, 謝謝交流 2014-01-10 11:17:44
感謝參與，應助指數(shù) +1
孫遠慧: 金幣+5, ★★★很有幫助 2014-01-10 09:28:20
sunyang1988: 金幣+2, 謝謝交流 2014-01-10 11:17:44
|
如果你看過vasp的源程序，就知道如果不改源程序的話，是不可能只改變邊長，不改變角度的。原因是ISIF=3相對于ISIF=0，是將晶格參數(shù)的9個變量作為和原子坐標一樣的參數(shù)并列，然后優(yōu)化。也就是他們和原子坐標等同，這樣他就有額外的9個獨立的自由度。你要的是限制某些自由度（角度），vasp里面是沒有類似的限制的。 解決辦法， 1 （推薦）. 如果想只改變邊長不改變角度，最好的辦法應該是寫個小腳本，腳本的寫法看vasp手冊最后。 2. （不推薦），改變源程序，找到你的用的算法那部分就是IBRION=1or 2or3 對應的腳本部分，加入你的限制條件，重新編譯。 |
5樓2014-01-09 12:26:26
孫遠慧
木蟲 (小有名氣)
- 應助: 20 (小學生)
- 金幣: 1450.3
- 散金: 282
- 紅花: 5
- 帖子: 229
- 在線: 224.3小時
- 蟲號: 1375586
- 注冊: 2011-08-22
- 性別: GG
- 專業(yè): 理論和計算化學
|
是110頁那個嗎？腳本能控制在計算過程中只改變20 45 7.5719這三個值的大小，而不改變其余六個的值嗎？不需要成比例變化，我只要保證還是原來的立方晶系就行。能不能麻煩你給我發(fā)一個滿足要求的腳本，或者告訴我在書上多少頁？syh5852767@163.com，謝了。。 |
6樓2014-01-10 10:40:25
kjsun
木蟲 (正式寫手)
- 應助: 34 (小學生)
- 金幣: 5557.5
- 散金: 100
- 紅花: 5
- 帖子: 434
- 在線: 469.8小時
- 蟲號: 285070
- 注冊: 2006-10-14
- 性別: GG
- 專業(yè): 催化化學
★
sunyang1988: 金幣+1, 謝謝交流 2014-01-17 13:52:36
sunyang1988: 金幣+1, 謝謝交流 2014-01-17 13:52:36
|
現(xiàn)在有網(wǎng)絡，我好久不看書了，不知道多少頁。 腳本參照這個例子 http://cms.mpi.univie.ac.at/vasp/guide/node174.html |
7樓2014-01-10 12:09:15
8樓2014-01-10 13:04:40
孫遠慧
木蟲 (小有名氣)
- 應助: 20 (小學生)
- 金幣: 1450.3
- 散金: 282
- 紅花: 5
- 帖子: 229
- 在線: 224.3小時
- 蟲號: 1375586
- 注冊: 2011-08-22
- 性別: GG
- 專業(yè): 理論和計算化學
10樓2014-01-16 22:05:48
相關版塊跳轉(zhuǎn)
我要訂閱樓主 孫遠慧 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +4 | 申子申申 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +8 | 手機用戶 2026-03-17 | 8/400 |
|
|
|
[考研] 329求調(diào)劑 +6 | 想上學吖吖 2026-03-19 | 6/300 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +17 | enenenhui 2026-03-13 | 18/900 |
|
|
|
[論文投稿]
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學材料求調(diào)劑，總分315（英一） +3 | sbdksD 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊 2026-03-18 | 5/250 |
|
|
|
[考研] 本科鄭州大學物理學院，一志愿華科070200學碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|