版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

月牙_2012

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 39.3
帖子: 9
在線: 12.1小時(shí)
蟲號(hào): 2180256
注冊(cè): 2012-12-11
專業(yè): 分離過(guò)程

[求助] 用lammps作PMF不平衡

各位大神，我嘗試用LAMMPS中fix smd命令做PMF。體系是中間為石墨烯，是石墨烯左右兩邊是NaCl溶液。準(zhǔn)備分別做鈉離子、氯離子和水的PMF.但運(yùn)行的結(jié)果，卻是左右不對(duì)稱。請(qǐng)各位大神幫我看看出現(xiàn)什么問題了。謝謝�。�！
units real
atom_style full
neighbor 2.0 bin
boundary p p p
read_data H-PMF.data

#lj potentials
pair_style  lj/cut/coul/long 10.0
pair_coeff 1 1 0.07 3.55
pair_coeff 2 2 0.07 3.55
pair_coeff 3 3 0.03 2.42
pair_coeff 4 4 0.1553 3.166
pair_coeff 5 5 0 0
pair_coeff 6 6 0.1 2.583
pair_coeff 7 7 0.1 4.401
pair_coeff 8 8 0.1553 3.166
pair_coeff 9 9 0 0
pair_coeff 10 10 0.1 2.583
pair_coeff 11 11 0.1 4.401
pair_coeff 12 12 0.1553 3.166
pair_coeff 13 13 0 0

pair_modify mix arithmetic
#special_bonds charmm

bond_style harmonic
bond_coeff 1 469 1.40
bond_coeff 2 367 1.080
bond_coeff 3 450 1.0
bond_coeff 4 450 1.0

angle_style harmonic
angle_coeff 1 63 120
angle_coeff 2 35 120
angle_coeff 3 55 109.47
angle_coeff 4 55 109.47

#special_bonds lj/coul 0.0 0.0 0.5
neigh_modify delay 0 every 1 check yes

kspace_style pppm 1.0e-4

#kspace_modify slab 3.0

special_bonds  lj 0.0 0.0 0.5 coul 0.0 0.0 0.5

#define groups
group GHE type 1 2 3
group shui1 type 4 5 12 13
group shui2 type 8 9
group LJ type 6 7 10 11
group shui union shui1 shui2
group and union shui LJ
group smd-w type 12 13

neigh_modify exclude group GHE GHE
#neigh_modify exclude type 9 9
#neigh_modify delay 0 every 1

#minimize 1.0e-4 1.0e-6 100 1000

#velocity and create 298.0 4928459 units box
velocity GHE set 0.0 0.0 0.0 units box

fix 1 GHE setforce 0 0 0
fix 2 and nvt temp 298.0 298.0 100.0
fix 3 GHE rigid/nvt group 1 GHE temp 298.0 298.0 100.0 force * off off off torque * off off off
fix 4 smd-w rigid group 1 smd-w force * off off on torque * off off off
fix 5 shui shake 0.0001 20 0 b 3 4 a 3 4
fix 6 smd-w smd cvel 20 -0.000005 tether NULL NULL 15 0.0

回復(fù)此樓

» 猜你喜歡

279分求調(diào)劑一志愿211 已經(jīng)有16人回復(fù)
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）已經(jīng)有5人回復(fù)
上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境）已經(jīng)有3人回復(fù)
急發(fā)核心期刊論文已經(jīng)有5人回復(fù)
291 求調(diào)劑已經(jīng)有5人回復(fù)
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有11人回復(fù)
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑已經(jīng)有7人回復(fù)
316求調(diào)劑已經(jīng)有7人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有4人回復(fù)
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑已經(jīng)有7人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

Lammps對(duì)單晶銅進(jìn)行弛豫，體積忽然增大，內(nèi)存異常，vmd查看發(fā)現(xiàn)體積基本不變已經(jīng)有13人回復(fù)
lammps 中做拉伸時(shí) sss和ssp邊界條件不同，結(jié)果為什么不一樣？？？？已經(jīng)有5人回復(fù)
求助LAMMPS兩次運(yùn)算結(jié)果不一致問題！已經(jīng)有11人回復(fù)
MS生成Lammps的data文件，文件中原子的類型和順序并不是自己想要的，請(qǐng)問如何修改？已經(jīng)有8人回復(fù)
MS里的forcite 模塊只能用自帶的那三個(gè)力場(chǎng)嗎？想用trappe力場(chǎng)是不是要用lammps? 已經(jīng)有6人回復(fù)
分子模擬中體系勢(shì)能和總能量為很大的正值已經(jīng)有8人回復(fù)
分子動(dòng)力學(xué)能量守恒的相關(guān)問題已經(jīng)有18人回復(fù)
想做粒子和固體（金屬）相互作用，不知道lammps能否實(shí)現(xiàn)？已經(jīng)有6人回復(fù)
lammps實(shí)驗(yàn)拉伸Al納米線的in文件的錯(cuò)誤提示好久找不出？已經(jīng)有8人回復(fù)
lammps 離子液體，離子中相鄰隔的挺遠(yuǎn)的原子成鍵了，不知什么原因已經(jīng)有8人回復(fù)
LAMMPS 的坐標(biāo)文件怎么構(gòu)建呢？不可能一個(gè)一個(gè)原子的輸入吧？已經(jīng)有8人回復(fù)
用MS建模導(dǎo)入LAMMPS計(jì)算，是不是不能用compass力場(chǎng)？已經(jīng)有9人回復(fù)
lammps求助：計(jì)算含有不同原子體系的熱導(dǎo)率的in文件該怎么修改已經(jīng)有7人回復(fù)
【求助】NVT系綜在MD模擬過(guò)程中的能量變化已經(jīng)有6人回復(fù)
【求助】LAMMPS中非周期性邊界條件和Born-Mayer-Huggins勢(shì)函數(shù)不會(huì)設(shè)置已經(jīng)有4人回復(fù)

1樓 2014-07-14 10:15:58

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

258568

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 33.5
帖子: 3
在線: 1小時(shí)
蟲號(hào): 6492561
注冊(cè): 2017-05-10
專業(yè): 微/納機(jī)械系統(tǒng)

請(qǐng)問你通過(guò)smd計(jì)算的pmf是用什么命令輸出的。

發(fā)自小木蟲Android客戶端

贊一下