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

本帖產(chǎn)生 1 個 1ST強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

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ù)此樓

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

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

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有9人回復(fù)
293求調(diào)劑已經(jīng)有14人回復(fù)
0703化學(xué)336分求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有9人回復(fù)
環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請已經(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ì)試驗(yàn)后，是否以及怎么樣進(jìn)行后續(xù)的進(jìn)一步優(yōu)化？已經(jīng)有6人回復(fù)
在vasp計(jì)算里，有沒有好的只優(yōu)化原子位置和晶胞體積的方法？已經(jīng)有8人回復(fù)
vasp 那些參數(shù)設(shè)置可以影響到磁性計(jì)算的磁矩已經(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ì)算超晶胞某個表面的輸入文件給我已經(jīng)有7人回復(fù)
【求助】Castep中結(jié)構(gòu)優(yōu)化后晶胞參數(shù)差別較大 vasp高手也幫我看一下已經(jīng)有22人回復(fù)
【求助】vasp編譯時，報(bào)錯：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

金蟲 (正式寫手)

★★★★★ 五星級,優(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)秀推薦

強(qiáng)烈支持

回復(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ù)文

如果不是非正交晶胞的話，該如何實(shí)現(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
如果不是非正交晶胞的話，該如何實(shí)現(xiàn)呢？...

同問，如果非正交呢

贊一下

回復(fù)此樓

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

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

zyqh

鐵桿木蟲 (著名寫手)

應(yīng)助: 33 (小學(xué)生)
金幣: 6698
帖子: 2375
在線: 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的回帖