| 2 | 1/1 | 返回列表 |
| 查看: 510 | 回復(fù): 1 | ||
牛很牛新蟲(chóng) (小有名氣)
|
[求助]
水分子溫度趨近于零
|
| 目前我遇到的問(wèn)題就是在模擬過(guò)程中,盡管我設(shè)置了分子的溫度為一個(gè)定值,但是隨著計(jì)算的進(jìn)行,水分子的溫度逐漸趨近于0。這就是目前我遇到的最尷尬,最頭疼的問(wèn)題,目前我自己查找參數(shù)的設(shè)置是沒(méi)有問(wèn)題,但一定是哪里出了問(wèn)題,希望內(nèi)行的大佬指點(diǎn)一二。稍后貼出代碼 |
新蟲(chóng) (小有名氣)
|
dimension 3 boundary p p p units real #atom_style full neighbor 2.0 bin neigh_modify every 2 delay 2 lattice fcc 1 region 1 block -52.0 252.0 -52.0 252.0 -52.0 252.0 #all pair_style lj/cut/coul/long 15.0 kspace_style pppm 1.0e-5 atom_style full bond_style harmonic angle_style harmonic special_bonds lj 0.0 0.0 0.0 coul 0.0 0.0 0.0 read_data new.data 這個(gè)是前面的設(shè)置 fix 10 suan rigid molecule velocity all create 400 49284 #velocity Cl create 600 49284 #velocity suan create 600 49284 timestep 0.1 ############### compute ############# compute 1 NaCl temp #compute 2 Cl temp compute 2 suan temp compute 3 water temp compute 4 all pe compute ake all ke compute energy all group/group suan thermo_style custom step temp c_1 c_2 c_3 #press vol etotal c_3 c_1 c_4 c_2 c_ake c_gke c_energy thermo 1000 thermo_modify temp 1 fix 1 all nvt temp 450.0 450.0 100.0 #fix rec all recenter INIT INIT INIT dump 1 all custom 1000 ch0110.water id type x y z dump 2 NaCl custom 1000 ch0111.water id type x y z dump 3 suan custom 1000 ch0112.water id type x y z run 10000 #unfix 1 #fix 1 all nvt temp 250.0 250.0 5.0 #fix 2 all nve #fix 3 jia nve #compute cluster all cluster/atom 5.5 #dump 4 all custom 1000 ch0113.water id type x y z #run 500000 write_data restart.water 這是后面的設(shè)置 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 085600材料調(diào)劑 總分330 +3 | 池池丶 2026-03-03 | 3/150 |
|
|
[考研] 本科太原理工采礦工程,求調(diào)劑 +3 | onlx 2026-03-01 | 3/150 |
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專(zhuān)業(yè)的同學(xué)報(bào)考 +9 | pin8023 2026-02-28 | 12/600 |
|
|
[考研] 0856材料與化工,270求調(diào)劑 +14 | YXCT 2026-03-01 | 16/800 |
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 11/550 |
|
|
[考研] 0854總分272 +5 | 打小就是老實(shí)人 2026-03-02 | 6/300 |
|
|
[考研]
|
旅行中的紫葡萄 2026-03-03 | 4/200 |
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
[基金申請(qǐng)] 請(qǐng)問(wèn)大家,研究風(fēng)險(xiǎn)與應(yīng)對(duì)措施那里, 大家都怎么寫(xiě)呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
[考研] 266材料化工求調(diào)劑 +3 | 哇塞王帥 2026-03-03 | 3/150 |
|
|
[考研] 290分材料工程085601求調(diào)劑 數(shù)二英一 +8 | llx0610 2026-03-02 | 9/450 |
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
[考研] 材料工程274求調(diào)劑 +5 | Lilithan 2026-03-01 | 5/250 |
|
|
[考研] 材料學(xué)調(diào)劑 +10 | 提神豆沙包 2026-02-28 | 12/600 |
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|
|
[考研]
|
LYidhsjabdj 2026-02-28 | 4/200 |
|