| 3 | 1/1 | 返回列表 |
| 查看: 1063 | 回復: 2 | ||
[求助]
不知道為什么錯,求指導
|
|
LAMMPS (1 Feb 2014) ### REAX potential for GO ### This script is prepared on Feb. 9, 2011 ### QPHLL ## Initialization units real dimension 3 boundary p p p atom_style charge ## Atom definition read_data data.CHO orthogonal box = (2.1226 -0.257523 -0.855987) to (77.1226 74.7425 24.704) 4 by 4 by 1 MPI processor grid reading atoms ... 3078 atoms #Force Field, in ReaxFF, 1-C, 2-H, 3-O, 4-N, 5-S pair_style reax pair_coeff * * ffield.reax 1 2 3 # Settings # mass 1 12.0107 # mass 2 1.00794 # mass 3 15.9994 neighbor 2.5 bin neigh_modify delay 0 every 2 check yes # TimeUnit=fs timestep 0.25 compute reax all pair reax variable eb equal c_reax[1] variable ea equal c_reax[2] variable elp equal c_reax[3] variable emol equal c_reax[4] variable ev equal c_reax[5] variable epen equal c_reax[6] variable ecoa equal c_reax[7] variable ehb equal c_reax[8] variable et equal c_reax[9] variable eco equal c_reax[10] variable ew equal c_reax[11] variable ep equal c_reax[12] variable efi equal c_reax[13] variable eqeq equal c_reax[14] ## Output frequency thermo 10 thermo_style custom step cpu press vol temp etotal pe ke v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb v_et v_eco v_ew v_ep v_efi v_eqeq thermo_modify line one restart 10000 restart.* # AtomEye, scaled to box size, xs, ys, zs will be [0,1] dump 1 all cfg 1000 dump.*.cfg mass type xs ys zs vx vy vz fx fy fz #dump_modify 1 element C H O dump_modify 1 sort id # VMD dump 2 all dcd 100 movie.dcd dump_modify 2 sort id # XYZ dump 3 all xyz 10000 structure.*.xyz dump_modify 3 sort id # Similar to ReaxFF_xmolout, further modification needed! dump 4 all custom 10 xmolout type x y z dump_modify 4 sort id # Charge information dump 5 all custom 10 charge.q id type q dump_modify 5 sort id # Charge informaiton with coordinates dump 6 all custom 10 coordinate.q id type x y z q dump_modify 6 sort id minimize 1.0e-6 1.0e-6 100000 100000 WARNING: Resetting reneighboring criteria during minimization (../min.cpp:173) |

金蟲 (著名寫手)

銀蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 本人考085602 化學工程 專碩 +17 | 不知道叫什么! 2026-03-15 | 19/950 |
|
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 一志愿華中科技大學,080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
[考研] 材料與化工一志愿南昌大學327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 331求調(diào)劑(0703有機化學 +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|
|
[考博] 福州大學楊黃浩課題組招收2026年專業(yè)學位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|