| 6 | 1/1 | 返回列表 |
| 查看: 3585 | 回復: 5 | ||
[求助]
用reaxff ,請大神幫忙看下我失敗的in和control文件 已有1人參與
|
|
in.file 如下 =========================================== boundary p p s units real atom_style charge #neighbor 2 bin #neigh_modify every 10 delay 0 check no read_data data.input pair_style reax/c lmp_control pair_modify mix arithmetic pair_coeff * * ffield Ti O C H fix 1 all qeq/reax 1 0.0 10.0 1.0e-6 reax/c fix 2 all nvt temp 10.0 10.0 100.0 min_style cg minimize 0.0 1.0e-8 1000 100000 #stopping tolerance for energy,force #max iterations,max number of force/energy evaluations #thermo 1000 # output thermodynamics every N timesteps #thermo_style multi write_restart restart.md clear boundary p p s units real atom_style charge neighbor 2 bin neigh_modify every 1 delay 0 check no read_restart restart.md pair_style reax/c lmp_control pair_modify mix arithmetic pair_coeff * * ffield Ti O C H fix 1 all qeq/reax 1 0.0 10.0 1.0e-6 reax/c fix 2 all nvt temp 10.0 100.0 100.0 #group Atom_fix id 1:312 #fix 3 Atom_fix setforce 0.0 0.0 0.0 group Ti type 1 fix 4 Ti setforce 0.0 0.0 0.0 timestep 1 thermo 10000 thermo_style multi fix walls all wall/reflect zlo EDGE zhi EDGE run 10000 unfix 2 fix 2 all nvt temp 100.0 100.0 100.0 dump 1 all xyz 1000 construction.xyz dump_modify 1 element Ti O C H dump_modify 1 sort id compute 3 all msd fix msd all ave/time 1 5 1000 c_3 file dump mode vector run 1000000 ================================================ control.file 如下 ========================================== simulation_name CHO_example ! output files will carry this name + their specific ext tabulate_long_range 1000 ! denotes the granularity of long range tabulation, 0 means no tabulation energy_update_freq 1 nbrhood_cutoff 1.5 ! near neighbors cutoff for bond calculations in A hbond_cutoff 6.0 ! cutoff distance for hydrogen bond interactions bond_graph_cutoff 0.3 ! bond strength cutoff for bond graphs thb_cutoff 0.001 ! cutoff value for three body interactions write_freq 1 ! write trajectory after so many steps traj_title CHO ! (no white spaces) atom_info 1 ! 0: no atom info, 1: print basic atom info in the trajectory file atom_forces 1 ! 0: basic atom format, 1: print force on each atom in the trajectory file atom_velocities 0 ! 0: basic atom format, 1: print the velocity of each atom in the trajectory file bond_info 1 ! 0: do not print bonds, 1: print bonds in the trajectory file angle_info 1 ! 0: do not print angles, 1: print angles in the trajectory file ================================================= 我計劃讓體系運行100W步,但是每次只能運行至80W步左右,這個問題一直困惑著我,按理說MD應該可以一直運行下去才對啊。。。 跪求大神解惑啊啊啊啊啊 ![]() ![]() ![]() ![]() |
金蟲 (正式寫手)
銅蟲 (小有名氣)

新蟲 (小有名氣)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調劑 +6 | 西南交通專材355 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 292求調劑 +9 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 10/500 |
|
|
[考研] 一志愿北京化工大學材料與化工(085600)296求調劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
[考研] 資源與環(huán)境 調劑申請(333分) +9 | holy J 2026-03-21 | 9/450 |
|
|
[考研]
|
孅華 2026-03-22 | 10/500 |
|
|
[考研] 274求調劑 +13 | 顧九笙要謙虛 2026-03-24 | 19/950 |
|
|
[考研] 環(huán)境專碩324分求調劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 一志愿哈工大,085400,320,求調劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 299求調劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
[考研] 求調劑 +3 | QiMing7 2026-03-25 | 3/150 |
|
|
[考研] 26考研-291分-廈門大學(085601)-柔性電子學院材料工程專業(yè)求調劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 招08考數(shù)學 +8 | laoshidan 2026-03-20 | 17/850 |
|
|
[考研] 【2026考研調劑】制藥工程 284分 求相關專業(yè)調劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 求調劑323材料與化工 +4 | 1124361 2026-03-24 | 4/200 |
|
|
[考研] 289求調劑 +9 | 懷瑾握瑜l 2026-03-20 | 9/450 |
|
|
[考研] 上海電力大學材料防護與新材料重點實驗室招收調劑研究生(材料、化學、電化學,環(huán)境) +4 | 我愛學電池 2026-03-23 | 4/200 |
|
|
[考研] 求調劑一志愿武漢理工大學材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考研] 344求調劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 環(huán)境學碩288求調劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調劑 +4 | 幸運的醬醬 2026-03-22 | 5/250 |
|