【獎勵】本帖被評價14次，作者chemvasp增加金幣 10.4 個

chemvasp

新蟲 (初入文壇)

[資源] VASP并行可執(zhí)行軟件包，可對晶胞參數(shù)進(jìn)行部分優(yōu)化

下載地址：http://ishare.iask.sina.com.cn/f/17331535.html

在優(yōu)化表面，二維或一維體系時，必須固定某個或某幾個晶胞方向不優(yōu)化。
VASP本身沒有這個功能，但是提供了constr_cell_relax.F文件來設(shè)置。
根據(jù)文件中的提示，添加幾行代碼后，就可以任意固定晶胞方向了。

使用方法：
在目錄下建立文件OPTCELL，第一行輸入三個數(shù)字（比如110），
數(shù)字1表示優(yōu)化，數(shù)字0表示不優(yōu)化；三個數(shù)字依次對應(yīng)于
xyz方向。比如110表示z軸不優(yōu)化；001表示x和y軸都不優(yōu)化。
注1. 無OPTCELL文件則完全是原版軟件的功能；
注2. 只適用與正交晶胞，所有角都是直角。

附：修改后的constr_cell_relax.F文件
   SUBROUTINE CONSTR_CELL_RELAX(FCELL)
   USE PREC
   REAL(Q) FCELL(3,3), SAVE(3)
   LOGICAL FILFLG
   INTEGER ICELL(3)

   INQUIRE(FILE='OPTCELL',EXIST=FILFLG)
   IF (FILFLG) THEN
      OPEN(67,FILE='OPTCELL',FORM='FORMATTED',STATUS='OLD')
      READ(67,"(3I1)"

(ICELL(I),I=1,3)
      CLOSE(67)
      DO I=1,3
      SAVE(I)=FCELL(I,I)
      ENDDO
      FCELL=0.0d0
      DO I=1,3
      IF (ICELL(I)==1) FCELL(I,I)=SAVE(I)
      ENDDO
   ENDIF

!    just one simple example
!    relaxation in x directions only
!    SAVE=FCELL(1,1)
!    FCELL=0 ! F90 style: set the whole array to zero
!    FCELL(1,1)=SAVE

   RETURN
   END SUBROUTINE

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

第一性原理計算-tg	vasp計算資料	第一性原理計算經(jīng)驗	材料計算模擬軟件
VASP資料	量化軟件	有效質(zhì)量	網(wǎng)貼收集
vasp	soft study	vasp

» 猜你喜歡

學(xué)碩274求調(diào)劑已經(jīng)有11人回復(fù)
一志愿鄭大材料工程290求調(diào)劑已經(jīng)有12人回復(fù)
環(huán)境工程 085701，267求調(diào)劑已經(jīng)有8人回復(fù)
083000學(xué)碩274求調(diào)劑已經(jīng)有10人回復(fù)
一志愿085600中科院寧波所276分求調(diào)劑已經(jīng)有4人回復(fù)
材料專業(yè)申博已經(jīng)有3人回復(fù)
307求調(diào)劑已經(jīng)有20人回復(fù)
285求調(diào)劑已經(jīng)有8人回復(fù)
材料與化工（0856）304求B區(qū)調(diào)劑已經(jīng)有8人回復(fù)
370求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

關(guān)于vasp收斂和優(yōu)化已經(jīng)有24人回復(fù)
ansys中接觸單元的參數(shù)如何定義已經(jīng)有10人回復(fù)
vasp.5.2.12 說明書已經(jīng)有4人回復(fù)
正交設(shè)計試驗后，是否以及怎么樣進(jìn)行后續(xù)的進(jìn)一步優(yōu)化？已經(jīng)有6人回復(fù)
在vasp計算里，有沒有好的只優(yōu)化原子位置和晶胞體積的方法？已經(jīng)有8人回復(fù)
vasp 那些參數(shù)設(shè)置可以影響到磁性計算的磁矩已經(jīng)有12人回復(fù)
dmol在對晶胞參數(shù)進(jìn)行優(yōu)化時是否應(yīng)該固定晶胞角度（fix angles）已經(jīng)有12人回復(fù)
利用CASTEP對建立的晶胞進(jìn)行幾何優(yōu)化的問題（急�。。。� 已經(jīng)有13人回復(fù)
【求助】VASP怎么做晶格參數(shù)的結(jié)構(gòu)優(yōu)化已經(jīng)有17人回復(fù)
【求助】vasp優(yōu)化晶格參數(shù)a 已經(jīng)有13人回復(fù)
【求助】侯老師的vasp入門的晶格優(yōu)化部分沒看懂，問一下已經(jīng)有25人回復(fù)
【求助】能否發(fā)一份用vasp計算超晶胞某個表面的輸入文件給我已經(jīng)有7人回復(fù)
【求助】Castep中結(jié)構(gòu)優(yōu)化后晶胞參數(shù)差別較大 vasp高手也幫我看一下已經(jīng)有22人回復(fù)
【求助】vasp編譯時，報錯：cannot find -lmkl_em64t 已經(jīng)有11人回復(fù)

1樓 2011-07-22 15:11:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

huangll99

木蟲 (職業(yè)作家)

應(yīng)助: 43 (小學(xué)生)
金幣: 775.9
帖子: 4051
在線: 1447.4小時
蟲號: 1521828

★★★★★ 五星級,優(yōu)秀推薦

請問樓主，如果三個數(shù)都選1，那么晶胞的邊的夾角會不會變化？

贊一下(1人)

回復(fù)此樓

7樓2012-06-28 09:13:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

cavediger

金蟲 (正式寫手)

1ST強帖: 9
應(yīng)助: 6 (幼兒園)
金幣: 709
帖子: 545
在線: 230.2小時
蟲號: 1255502

★★★★★ 五星級,優(yōu)秀推薦

挺好，有創(chuàng)新，支持

回復(fù)此樓

2樓2011-07-22 17:08:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

雁兒霏霏

木蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2645.4
帖子: 883
在線: 160.1小時
蟲號: 934537

★★★★★ 五星級,優(yōu)秀推薦

強烈支持

回復(fù)此樓

3樓2011-07-23 08:08:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wutianhuan

鐵桿木蟲 (著名寫手)

應(yīng)助: 3 (幼兒園)
金幣: 6697.4
帖子: 1723
在線: 216.2小時
蟲號: 284387

★★★★★ 五星級,優(yōu)秀推薦

謝謝，相信會很有用的，呵呵

贊一下

回復(fù)此樓

4樓2011-07-23 12:21:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dycding

鐵桿木蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 5724.4
帖子: 11
在線: 592.6小時
蟲號: 115379

★★★ 三星級,支持鼓勵

請問我在CentOS5.5上按說明方法安裝時出現(xiàn)了無法執(zhí)行二進(jìn)制文件vasp52mpi.x（cannot execute binary file）這是什么原因？有什么方法可以解決嗎？謝謝！

贊一下

回復(fù)此樓

5樓2011-07-30 15:47:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dragon5061

至尊木蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 11731.9
帖子: 842
在線: 531.2小時
蟲號: 631715

引用回帖:

572389樓: Originally posted by chemvasp at 2011-07-22 15:11:45:
下載地址：http://ishare.iask.sina.com.cn/f/17331535.html

在優(yōu)化表面，二維或一維體系時，必須固定某個或某幾個晶胞方向不優(yōu)化。
VASP本身沒有這個功能，但是提供了constr_cell_relax.F文件來設(shè)置。
根據(jù)文

如果不是非正交晶胞的話，該如何實現(xiàn)呢？

贊一下

回復(fù)此樓

6樓2012-05-22 11:17:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ljzhou86

金蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 1000.9
帖子: 597
在線: 640.7小時
蟲號: 996220

引用回帖:

572389樓: Originally posted by chemvasp at 2011-07-22 15:11:45
下載地址：http://ishare.iask.sina.com.cn/f/17331535.html

在優(yōu)化表面，二維或一維體系時，必須固定某個或某幾個晶胞方向不優(yōu)化。
VASP本身沒有這個功能，但是提供了constr_cell_relax.F文件來設(shè)置。
根據(jù)文 ...

你好，我優(yōu)化一個2D構(gòu)型，(011)，x軸不優(yōu)化，放寬y,z軸優(yōu)化，但最終x,y,z都在優(yōu)化，似乎這于這種情況不適合，你檢查下代碼是哪有問題？當(dāng)然如果是(110)這樣的話，程序就沒有問題。

贊一下

回復(fù)此樓

9樓2012-08-01 00:07:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ljzhou86

金蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 1000.9
帖子: 597
在線: 640.7小時
蟲號: 996220

★★★★★ 五星級,優(yōu)秀推薦

贊一下

回復(fù)此樓

10樓2012-08-01 00:34:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hbl626

鐵蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 19.8
帖子: 330
在線: 54.2小時
蟲號: 546036

★★★★★ 五星級,優(yōu)秀推薦

太牛了，謝謝

回復(fù)此樓

11樓2012-08-18 16:58:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liancsh

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 124.6
帖子: 101
在線: 98.9小時
蟲號: 927362

★★★★★ 五星級,優(yōu)秀推薦

非常感謝！

回復(fù)此樓

12樓2012-11-07 08:27:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

duxiu727

金蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 648.6
帖子: 88
在線: 91.1小時
蟲號: 1677460

★★★★★ 五星級,優(yōu)秀推薦

引用回帖:

6樓: Originally posted by dragon5061 at 2012-05-22 11:17:08
如果不是非正交晶胞的話，該如何實現(xiàn)呢？...

同問，如果非正交呢

贊一下

回復(fù)此樓

13樓2013-10-21 22:12:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zyqh

鐵桿木蟲 (著名寫手)

應(yīng)助: 33 (小學(xué)生)
金幣: 6708
帖子: 2377
在線: 593.4小時
蟲號: 2762045

★ 一星級,一般

頂頂更健康

回復(fù)此樓

14樓2015-01-01 22:55:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangjl580

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 3.8
帖子: 42
在線: 15.5小時
蟲號: 674885

引用回帖:

10樓: Originally posted by ljzhou86 at 2012-08-01 00:34:09
你好，我優(yōu)化一個2D構(gòu)型，(011)，x軸不優(yōu)化，放寬y,z軸優(yōu)化，但最終x,y,z都在優(yōu)化，似乎這于這種情況不適合，你檢查下代碼是哪有問題？當(dāng)然如果是(110)這樣的話，程序就沒有問題。

我也發(fā)現(xiàn)了 011肯定是有問題的代碼有問題暫時還沒學(xué)程序，看不懂哪里有問題

贊一下

回復(fù)此樓

15樓2015-01-31 19:44:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

藍(lán)凌小未

銅蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1062.9
帖子: 460
在線: 106.4小時
蟲號: 2585093

★★★★★ 五星級,優(yōu)秀推薦

引用回帖:

你好~我也遇到了這個問題，請問你是如何解決的，謝謝大神蟲友~

贊一下

回復(fù)此樓

16樓2019-05-14 16:01:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

MichaelLai

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 259.6
帖子: 134
在線: 49小時
蟲號: 5358514

★★★★★ 五星級,優(yōu)秀推薦

請問老師，這個攻略里面重新編譯VASP是什么意思哇...具體如何操作呢

贊一下

回復(fù)此樓

17樓2024-07-22 16:30:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyikeco

禁蟲

本帖內(nèi)容被屏蔽

18樓2024-07-23 00:35:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復(fù)

sg184089268樓

2012-07-25 16:41 回復(fù)

五星好評

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求調(diào)劑 +3	圖鑒212 2026-03-30	3/150	2026-03-31 00:31 by jp9609
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +10	dongfang59 2026-03-27	10/500	2026-03-30 23:42 by 果果媽咪
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +3	披星河 2026-03-30	3/150	2026-03-30 23:38 by 果果媽咪
[有機交流] 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機理 10+3	kaobao456 2026-03-29	4/200	2026-03-30 23:16 by nBu鋰
[考研] 277跪求調(diào)劑 +8	1915668 2026-03-27	12/600	2026-03-30 21:01 by dophin1985
[考研] 332求調(diào)劑 +14	032500 2026-03-25	14/700	2026-03-30 20:03 by 無際的草原
[考研] 105500藥學(xué)求調(diào)劑，一志愿山東大學(xué)藥學(xué)，348分 +3	gr哈哈哈 2026-03-28	3/150	2026-03-30 18:56 by 源_2020
[考研] 310求調(diào)劑 +10	爭取九點睡 2026-03-30	10/500	2026-03-30 16:45 by ztnimte
[考研] 化學(xué)308分調(diào)劑 +7	你好明天你好 2026-03-30	8/400	2026-03-30 13:02 by 你好明天你好
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑 +5	篤行致遠(yuǎn). 2026-03-28	5/250	2026-03-30 06:56 by ilovexiaobin
[考研] 求調(diào)劑 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐兒
[考研] 329求調(diào)劑 +10	鈕恩雪 2026-03-25	10/500	2026-03-29 13:32 by peike
[考研] 352分-085602-一志愿985 +5	海納百川Ly 2026-03-29	5/250	2026-03-29 09:57 by Sjndkwm
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 343求調(diào)劑 +5	愛羈絆 2026-03-28	5/250	2026-03-28 20:53 by 唐沐兒
[考研] 一志愿上海理工能源動力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4	哇呼哼呼哼 2026-03-26	4/200	2026-03-27 17:02 by zhyzzh
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3	zzffylgg 2026-03-24	3/150	2026-03-27 10:42 by shangxh
[考研] 各位老師您好：本人初試372分 +5	jj涌77 2026-03-25	6/300	2026-03-25 14:15 by mapenggao