 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (2961)
- >導(dǎo)師招生 (269)
- >考研 (189)
- >考博 (81)
- >休閑灌水 (77)
- >蟲(chóng)友互識(shí) (66)
- >文獻(xiàn)求助 (54)
- >基金申請(qǐng) (41)
- >碩博家園 (23)
- >論文投稿 (14)
- >招聘信息布告欄 (13)
- >農(nóng)林 (8)
- >綠色求助(高懸賞) (8)
- >SciFinder/Reaxys (6)
- >找工作 (6)
- >博后之家 (5)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 3186 | 回復(fù): 1 | ||
月牙_2012新蟲(chóng) (初入文壇)
|
[求助]
用lammps作PMF不平衡
|
|
各位大神，我嘗試用LAMMPS中fix smd命令做PMF。體系是中間為石墨烯，是石墨烯左右兩邊是NaCl溶液。準(zhǔn)備分別做鈉離子、氯離子和水的PMF.但運(yùn)行的結(jié)果，卻是左右不對(duì)稱(chēng)。請(qǐng)各位大神幫我看看出現(xiàn)什么問(wè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ù)此樓» 猜你喜歡
招08考數(shù)學(xué)
已經(jīng)有11人回復(fù)
-
291求調(diào)劑
已經(jīng)有4人回復(fù)
-
276求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿070300浙大化學(xué)358分，求調(diào)劑！
已經(jīng)有4人回復(fù)
-
一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有5人回復(fù)
-
323求調(diào)劑
已經(jīng)有6人回復(fù)
-
352求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有8人回復(fù)
-
277材料科學(xué)與工程080500求調(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é)果不一致問(wèn)題！
已經(jīng)有11人回復(fù)
- MS生成Lammps的data文件，文件中原子的類(lèi)型和順序并不是自己想要的，請(qǐng)問(wèn)如何修改？
已經(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)問(wè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ù)
258568
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 33.5
- 帖子: 3
- 在線: 1小時(shí)
- 蟲(chóng)號(hào): 6492561
- 注冊(cè): 2017-05-10
- 專(zhuān)業(yè): 微/納機(jī)械系統(tǒng)
2樓2017-06-09 16:00:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 月牙_2012 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +6 | 生物工程調(diào)劑 2026-03-17 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 0856材料專(zhuān)碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +4 | 要好好無(wú)聊 2026-03-21 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂(lè)29 2026-03-18 | 5/250 |
|
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|