 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1851 | 回復: 1 | ||||
oxox6085專家顧問 (正式寫手)
|
[求助]
求助�。�！L-J勢能如何設(shè)置？老是算不對！
|
|
請教大家一個關(guān)于L-J的問題： 我第一次用GULP，第一步用它先優(yōu)化一下CdSe晶體的結(jié)構(gòu)（Zincblende），用文獻中的potential（這個potential應該很可靠，很多人用），基本是修改了一下example1的輸入文件，但是算出來的結(jié)果總是和實驗對不上，晶體的晶格常數(shù)總是特別小，有時候甚至輸出文件直接告訴我說cell的變化超過范圍，然后就報錯了。。那個potential很簡單，就是：庫侖+LJ勢，庫侖作用gulp應該直接用eward summation，精度默認10^-8應該沒有問題。所以，我總感覺我的LJ勢輸入的有問題。以下是potential的文獻中數(shù)據(jù)： 文中的sigma的單位是唉，epsilon是K，gulp里sigma是唉，但是epsilon用的是eV，所以需要轉(zhuǎn)化一下epsilon的單位：1(1/kB) eV=11604.5 K，所以我的輸入文件如下： =======================分割線======================= opti prop conp title CdSe-ZB end cell 6.0500 6.05000 6.05000 90.0 90.0 90.0 frac Cd 0.000000 0.000000 0.000000 Se 0.250000 0.250000 0.250000 space 216 species Cd 1.18 Se -1.18 lennard zero Cd Cd 0.001448 1.98 0.0 12.0 Se Se 0.001285 5.24 0.0 12.0 Cd Se 0.001364 3.61 0.0 12.0 output xr CdSe-ZB output marvin CdSe-ZB.mvn =======================分割線======================= 文獻里面公式中有系數(shù)4，對應的是L-J勢中的zero，那個很小的數(shù)值（比如0.001448）就是經(jīng)過單位轉(zhuǎn)換的epsilon，Cd-Se之間的參數(shù)是通過combining rules手工輸入的（也就是epsilon取相乘開根，sigma相加除以二，gulp里應該可以自動算，還不清楚如何設(shè)置）。。。。但是結(jié)果老是算不對，求大家?guī)椭�，看看我哪里出錯了。。。多謝多謝��！ |
回復此樓» 猜你喜歡
308求調(diào)劑
已經(jīng)有3人回復
-
321求調(diào)劑
已經(jīng)有11人回復
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有8人回復
-
330求調(diào)劑
已經(jīng)有4人回復
-
304求調(diào)劑
已經(jīng)有7人回復
-
一志愿南京理工大學085701資源與環(huán)境302分求調(diào)劑
已經(jīng)有6人回復
-
南京大學化學376求調(diào)劑
已經(jīng)有3人回復
-
考研調(diào)劑求學校推薦
已經(jīng)有5人回復
-
323求調(diào)劑
已經(jīng)有3人回復
-
317求調(diào)劑
已經(jīng)有12人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 有什么軟件可以方便地計算出特定分子的Lennard-Jones勢能和碰撞直徑？
已經(jīng)有5人回復
- 如何利用GULP擬合分子間勢能曲線得出力場參數(shù)
已經(jīng)有10人回復
- 時間步長設(shè)置，截斷半徑設(shè)置
已經(jīng)有6人回復
- 對于MD，其中的勢函數(shù)參數(shù)選取，如何自己計算
已經(jīng)有5人回復
- ATP分子力場參數(shù)怎么做
已經(jīng)有4人回復
- 能不能用量子化學計算結(jié)果擬合Lennard-Jones勢，怎么擬合？
已經(jīng)有5人回復
- 關(guān)于勢函數(shù)的求助
已經(jīng)有23人回復
- 【求助】如何通過理論計算得到分子各態(tài)能級以及勢能曲線或曲面？
已經(jīng)有4人回復
oxox6085
專家顧問 (正式寫手)
-
專家經(jīng)驗: +129 - 模擬EPI: 1
- 應助: 69 (初中生)
- 金幣: 8061
- 散金: 222
- 紅花: 39
- 帖子: 473
- 在線: 636.8小時
- 蟲號: 612537
- 注冊: 2008-09-25
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
- 管轄: 分子模擬
|
沒人回答，好吧我自問自答一下，方便蟲友（好苦逼） 我后來發(fā)現(xiàn)我的input文件關(guān)于LJ勢那塊還是格式不對： 1，如果用A B的形式，需要： lennard Atom1 Atom2 A B Rmin Rmax flag1 flag2 2，如果是epsilon/sigma的格式： lennard epsilon Atom1 Atom2 epsilon sigma Rmin Rmax flag1 flag2 3，如果是epsilon zero的格式： lennard epsilon zero Atom1 Atom2 episilon sigma Rmin Rmax flag1 flag2 epsilon的時候sigma給的是勢能最低點，epsilon zero的時候sigma給的是勢能為0是的距離，mannual上有寫。 |
2樓2012-03-07 18:27:54
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 oxox6085 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
|
[考研] 一志愿北京化工大學0703化學318分，有科研經(jīng)歷，求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+7
|
.6lL 2026-03-18 | 7/350 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 0703化學調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|