| 6 | 1/1 | 返回列表 |
| 查看: 3586 | 回復(fù): 5 | ||
[求助]
用reaxff ,請(qǐng)大神幫忙看下我失敗的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 ================================================= 我計(jì)劃讓體系運(yùn)行100W步,但是每次只能運(yùn)行至80W步左右,這個(gè)問題一直困惑著我,按理說MD應(yīng)該可以一直運(yùn)行下去才對(duì)啊。。。 跪求大神解惑啊啊啊啊啊 ![]() ![]() ![]() ![]() |
金蟲 (正式寫手)
銅蟲 (小有名氣)

新蟲 (小有名氣)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +5 | 路癡小琪 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑推薦 材料 304 +10 | 荷包蛋hyj 2026-03-26 | 10/500 |
|
|
[考研] 334求調(diào)劑 +3 | 雨清天晴 2026-03-21 | 3/150 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 340求調(diào)劑 +3 | jhx777 2026-03-27 | 3/150 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 286求調(diào)劑 +13 | Faune 2026-03-21 | 13/650 |
|
|
[考研] 291求調(diào)劑 +9 | hhhhxn.. 2026-03-23 | 9/450 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 297求調(diào)劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
[考研] 315分求調(diào)劑 +5 | 26考研上岸版26 2026-03-26 | 5/250 |
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-25 | 7/350 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|