 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 2263 | 回復(fù): 8 | ||
[求助]
lammps 離子液體，離子中相鄰隔的挺遠(yuǎn)的原子成鍵了，不知什么原因
|
|
我在進(jìn)行模擬的時(shí)候，首先用兩個(gè)分子（四個(gè)離子）測(cè)試了，結(jié)構(gòu)非常好，沒(méi)啥異常。但是增多原子后，發(fā)現(xiàn)其中一個(gè)陰離子（tfsi-）中F原子和相隔非常遠(yuǎn)的O原子（相隔五個(gè)原子）成鍵了，從而使整個(gè)離子變的亂起八糟。 我檢查了力場(chǎng)參數(shù)，采用的是opls—aa勢(shì)，和文獻(xiàn)中一致，電荷也守恒。 想問(wèn)問(wèn)各位蟲(chóng)子，這個(gè)是什么原因呢？該怎么解決呢？ 也不知怎么搞的，我遇到好多詭異的問(wèn)題。。。 |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
凌楓霜葉» 猜你喜歡
求材料調(diào)劑
已經(jīng)有6人回復(fù)
-
085600材料與化工
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有6人回復(fù)
-
286求調(diào)劑
已經(jīng)有9人回復(fù)
-
焦慮
已經(jīng)有12人回復(fù)
-
344求調(diào)劑
已經(jīng)有6人回復(fù)
-
266求調(diào)劑
已經(jīng)有9人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有18人回復(fù)
-
314求調(diào)劑
已經(jīng)有8人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- lammps進(jìn)行分子動(dòng)力學(xué)模擬之前是不是一定要進(jìn)行結(jié)構(gòu)優(yōu)化？
已經(jīng)有4人回復(fù)
- 誠(chéng)聘鋰電研發(fā)人才！
已經(jīng)有7人回復(fù)
- 求助？怎么從劍橋晶體庫(kù)得到物質(zhì)的結(jié)構(gòu)數(shù)據(jù)
已經(jīng)有5人回復(fù)
- 關(guān)于模擬離子液體力場(chǎng)的選擇（opls-aa)
已經(jīng)有9人回復(fù)
- 【討論】lammps data 文件
已經(jīng)有3人回復(fù)
ddll_slsl
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 101 (高中生)
- 金幣: 2939.2
- 紅花: 6
- 帖子: 670
- 在線: 255.3小時(shí)
- 蟲(chóng)號(hào): 2121784
- 注冊(cè): 2012-11-12
- 專(zhuān)業(yè): 動(dòng)力學(xué)與控制
2樓2013-02-04 10:27:01
3樓2013-02-04 11:15:01
|
我模擬的體系是emim陽(yáng)離子，和fsi陰離子。 我的in文件如下： LAMMPS input file units real boundary p p p atom_style full #uncharged molecules, with bond, angle, dihedral, improper dihedral bond_style harmonic angle_style harmonic dihedral_style opls improper_style cvff pair_style lj/cut/coul/long 12.0 pair_modify mix arithmetic tail yes kspace_style pppm 1.0e-4 #read_restart restart.mix.* read_data data.il #read_restart restart.mix.* #fix 1 all shake 0.0001 20 10 b 3 neighbor 2.0 bin neigh_modify delay 3 every 1 check yes #delay 3 velocity all create 300.0 4567712 rot yes dist gaussian compute 1 all temp #fix 2 all npt temp 300.0 300.0 100.0 iso 1.0 1.0 10.0 #fix 2 all langevin 300.0 300.0 200 48279 #fix 3 all nve fix 2 all nvt temp 300.0 300.0 50.0 minimize 1.0e-4 1.0e-6 100 1000 min_style cg min_modify dmax 0.2 timestep 0.05 dump 1 all custom 20 dump.atom.*.lammpstrj id x y z dump 2 all custom 200000 dump.force.* id mol type element fx fy fz #dump 3 all dcd 2000 dump.trj.lammpstrj dump_modify 1 flush yes dump_modify 2 flush yes #dump_modify 3 flush yes thermo_style custom step c_1 pe etotal vol lx ly lz thermo 1000 thermo_modify flush yes restart 100000 restart.mix.* run 100000 timestep 0.1 #thermo 1000 run 100000 我的log文件 24000 273.6706 272.92199 317.7888 125000 50 50 50 25000 325.54712 263.24033 316.61201 125000 50 50 50 26000 313.35585 255.15618 306.52916 125000 50 50 50 27000 3516.0726 1330.2863 1906.7272 125000 50 50 50 28000 1783.7174 315.87659 608.30734 125000 50 50 50 29000 1473.6364 568.98349 810.57815 125000 50 50 50 30000 1474.526 -32.645347 209.09516 125000 50 50 50 31000 671.23014 38.678556 148.72308 125000 50 50 50 32000 464.68752 -90.201486 -14.018499 125000 50 50 50 33000 11929.712 727.72594 2683.5372 125000 50 50 50 34000 1249.3459 -505.03535 -300.2119 125000 50 50 50 35000 6864.2858 180.49005 1305.8523 125000 50 50 50 36000 418.5098 -612.91378 -544.30138 125000 50 50 50 37000 174.50306 -677.42917 -648.82034 125000 50 50 50 38000 79.074924 -680.99101 -668.0271 125000 50 50 50 39000 149.36694 -684.89448 -660.40658 125000 50 50 50 40000 34276.857 762.54228 6382.0463 125000 50 50 50 41000 405.71061 -658.47541 -591.96137 125000 50 50 50 42000 52.579346 -689.37579 -680.7557 125000 50 50 50 43000 37.702359 -697.65724 -691.47615 125000 50 50 50 44000 77.636161 -698.36347 -685.63545 125000 50 50 50 45000 1068.4743 -541.76354 -366.593 125000 50 50 50 46000 769.90494 -593.81149 -467.58977 125000 50 50 50 47000 909.93159 -769.17108 -619.99276 125000 50 50 50 Rank 19 [Mon Feb 4 02:49:12 2013] [c13-3c2s3n1] application called MPI_Abort(MPI_COMM_WORLD, 1) - process 19 ERROR on proc 19: Bond atoms 26 27 missing on proc 19 at step 47662 (neigh_bond.cpp:49) 從第3000步開(kāi)始，能量突然驟減，我調(diào)出結(jié)構(gòu)圖來(lái)一看，是結(jié)構(gòu)出現(xiàn)了問(wèn)題 |
4樓2013-02-04 11:22:57
| 離子結(jié)構(gòu)圖和結(jié)果圖如下， fsi中的F原子和emim中的Hoverlap了，導(dǎo)致能量變?yōu)樨?fù)數(shù)了。 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 結(jié)構(gòu)1.jpg
- 附件 2 : result.jpg
2013-02-04 11:30:25, 11.79 K
2013-02-04 11:30:57, 21.26 K
5樓2013-02-04 11:32:15
ddll_slsl
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 101 (高中生)
- 金幣: 2939.2
- 紅花: 6
- 帖子: 670
- 在線: 255.3小時(shí)
- 蟲(chóng)號(hào): 2121784
- 注冊(cè): 2012-11-12
- 專(zhuān)業(yè): 動(dòng)力學(xué)與控制
6樓2013-02-05 21:59:50
7樓2013-02-08 09:06:00
xiuyeyangyun
鐵蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 24.5
- 帖子: 24
- 在線: 8.5小時(shí)
- 蟲(chóng)號(hào): 1259687
- 注冊(cè): 2011-04-09
- 專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
8樓2013-05-08 15:53:23
凌楓霜葉
木蟲(chóng) (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1952.4
- 散金: 89
- 帖子: 228
- 在線: 111.6小時(shí)
- 蟲(chóng)號(hào): 2013857
- 注冊(cè): 2012-09-19
- 性別: GG
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
9樓2019-04-03 18:34:39
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wup123 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +6 | Liwangman 2026-03-15 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
|
[考研] 考研調(diào)劑 +4 | 芬達(dá)46 2026-03-12 | 4/200 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|