| 6 | 1/1 | 返回列表 |
| 查看: 4900 | 回復(fù): 5 | ||
jidanaoyun新蟲 (初入文壇)
|
[求助]
lammps中LJ勢能下單位問題 已有1人參與
|
|
求助一個問題,抱著一點小希望能有前輩給予指點。金幣不多,海涵。。。 在lammps中,系統(tǒng)中有兩種原子A和B,原子間勢能都是L-J勢,系統(tǒng)單位為units lj 由于在lj單位制下,物理量都是無量綱的,其中 長度單位就是LJ勢中的參數(shù) sigma 此時比如我定義了盒子大小為[0 0 0 ]-[1 1 1],那盒子的邊長是 1*sigma A 還是 sigma B 呢 ? 還是說,我只能用real的單位制? 讀便手冊也沒找到答案 哎。。。。 |
木蟲 (小有名氣)
|
A和B是由你決定的,你愿意讓哪個做單位都行,關(guān)鍵是單位選定之后,做單位的那個就是1,另外一個相應(yīng)改變。比如sigma_A=2sigma_B,如果A作單位,B就是0.5;如果B作單位,A就是2。另外設(shè)置的時候注意各自的截斷半徑以及A與B發(fā)生作用時的mixing rule。 [ 發(fā)自手機版 http://www.gaoyang168.com/3g ] |

新蟲 (初入文壇)
木蟲 (小有名氣)

新蟲 (初入文壇)
新蟲 (初入文壇)
|
請問用lj勢可以不涉及集體的原子么?如ar,H,C等。不涉及兩種的具體原子該怎么寫?我不知道怎么寫,就用了具體的晶格常數(shù),我這樣寫的對么? # LJ 12-6 system with npt ensemble units lj atom_style atomic boundary p p p lattice fcc 1.073 lattice fcc 0.8 region box block 0 10 0 5 0 5 create_box 2 box create_atoms 1 box create_atoms 2 box mass 1 1.0 mass 2 1.0 velocity all create 0.01 872877 timestep 0.01 dump 1 all xyz 1000 lj.xyz pair_style lj/cut 2.5 pair_coeff 1 1 1.0 1.0 2.5 pair_coeff 2 2 1.6 1.1765 2.5 pair_coeff 1 2 0.2 1.0882 2.5 #pair_modify shift yes mix geometric neighbor 0.3 bin neigh_modify every 10 delay 0 check yes thermo 1000 fix 1 all npt temp 0.01 0.01 1.00 iso 0.0024 0.0024 100.0 run 50000 unfix 1 thermo 1000 運行結(jié)果總出現(xiàn)nan怎么回事?我想不用到集體的原子怎么寫?謝謝 |
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +3 | SYA! 2026-03-23 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研] 招08考數(shù)學(xué) +6 | laoshidan 2026-03-20 | 14/700 |
|
|
[考研] 333求調(diào)劑 +6 | 87639 2026-03-21 | 10/500 |
|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項目,四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|