| 5 | 1/1 | 返回列表 |
| 查看: 3588 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
用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 ================================================= 我計(jì)劃讓體系運(yùn)行100W步,但是每次只能運(yùn)行至80W步左右,這個(gè)問題一直困惑著我,按理說MD應(yīng)該可以一直運(yùn)行下去才對啊。。。 跪求大神解惑啊啊啊啊啊 ![]() ![]() ![]() ![]() |
金蟲 (正式寫手)
銅蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +3 | 阿姨復(fù)古不過 2026-03-27 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 13/650 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +7 | 大饞小子 2026-03-27 | 8/400 |
|
|
[考研] 291求調(diào)劑 +6 | 孅華 2026-03-22 | 6/300 |
|
|
[考研] 289求調(diào)劑 +7 | 新時(shí)代材料 2026-03-27 | 7/350 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +14 | 荷包蛋hyj 2026-03-26 | 14/700 |
|
|
[考研] 340求調(diào)劑 +4 | jhx777 2026-03-27 | 4/200 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 283求調(diào)劑(080500) +4 | A child 2026-03-27 | 4/200 |
|
|
[考研] 317求調(diào)劑 +5 | 十閑wx 2026-03-24 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國家線 +10 | 哈哈157349 2026-03-21 | 10/500 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +8 | zhn03 2026-03-25 | 9/450 |
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|