| 4 | 1/1 | 返回列表 |
| 查看: 2146 | 回復(fù): 3 | |||
xiaowu787木蟲 (正式寫手)
|
[交流]
【求助】LAMMPS-REAX計算出錯 已有2人參與
|
|
剛學(xué)REAX,想自己做個輸入文件測試一下,結(jié)果出現(xiàn)如下錯誤: LAMMPS (24 Apr 2010) Reading data file ... orthogonal box = (0 0 0) to (3 8 2) 1 by 2 by 1 processor grid 4 atoms Error or end-of-file reading unit 4 on line: 55 rank 0 in job 268 pc07_59580 caused collective abort of all ranks exit status of rank 0: killed by signal 9 應(yīng)該是輸入文件沒有讀完就出現(xiàn)這個錯誤,請高手幫忙看看哪里存在錯誤? 由于只是為了測試,就做了一個1*1*1的MoO3的單胞,F(xiàn)在只是想測試程序的正常運算,力場里的數(shù)據(jù)可能不正確。 in文件 # this run is equivalent to GRASP testtatb units real atom_style charge read_data data.mo pair_style reax 10.0 1.0e-6 pair_coeff * * ffield.reax 1 2 neighbor 1.0 bin neigh_modify delay 0 every 1 check yes fix 1 all nve #Run one minimization to settle the lattice min_style cg minimize 1.0e-6 1.0e-6 100000 1000000 #print the entire lattice dump startdump all atom 2 autest.dump run 100 #but we only wanted printed out now undump startdump data 文件 #Position data for MoO3, 1 molecules 4 atoms 2 atom types 0.0 3.0 xlo xhi 0.0 8.0 ylo yhi 0.0 2.0 zlo zhi Masses 1 15.9990 2 95.9400 Atoms 1 1 0 1.72382 6.02831 0.92410 2 2 0 0.39584 1.40822 0.92410 3 1 0 2.06541 1.19984 0.92410 4 1 0 0.14781 3.06749 0.92410 ffield文件 Reactive MD-force field 39 ! Number of general parameters 50.0000 !Overcoordination parameter 9.5469 !Overcoordination parameter 127.8302 !Valency angle conjugation parameter 3.0000 !Triple bond stabilisation parameter 6.5000 !Triple bond stabilisation parameter 0.0000 !C2-correction 1.0496 !Undercoordination parameter 9.0000 !Triple bond stabilisation parameter 11.5054 !Undercoordination parameter 13.4059 !Undercoordination parameter 0.0000 !Triple bond stabilization energy 0.0000 !Lower Taper-radius 10.0000 !Upper Taper-radius 2.8793 !Not used 33.8667 !Valency undercoordination 7.0994 !Valency angle/lone pair parameter 1.0563 !Valency angle 2.0384 !Valency angle parameter 6.1431 !Not used 6.9290 !Double bond/angle parameter 0.3989 !Double bond/angle parameter: overcoord 3.9954 !Double bond/angle parameter: overcoord -2.4837 !Not used 5.7796 !Torsion/BO parameter 10.0000 !Torsion overcoordination 1.9487 !Torsion overcoordination -1.2327 !Conjugation 0 (not used) 2.1645 !Conjugation 1.5591 !vdWaals shielding 0.1000 !Cutoff for bond order (*100) 2.0038 !Valency angle conjugation parameter 0.6121 !Overcoordination parameter 1.2172 !Overcoordination parameter 1.8512 !Valency/lone pair parameter 0.5000 !Not used 20.0000 !Not used 5.0000 !Molecular energy (not used) 0.0000 !Molecular energy (not used) 3.6942 !Valency angle conjugation parameter 2 ! Nr of atoms; cov.r; valency;a.m;Rvdw;Evdw;gammaEEM;cov.r2;# alfa;gammavdW;valency;Eunder;Eover;chiEEM;etaEEM;n.u. cov r3;Elp;Heat inc.;n.u.;n.u.;n.u.;n.u. ov/un;val1;n.u.;val3,vval4 O 1.3825 4.0000 12.0000 1.9133 0.1853 0.9000 1.1395 4.0000 9.7602 2.1346 4.0000 33.2433 79.5548 5.8678 7.0000 0.0000 1.2104 0.0000 199.0303 8.6699 34.7289 13.3894 0.8563 0.0000 -2.8983 2.5000 1.0564 4.0000 2.9663 0.0000 0.0000 0.0000 Mo 2.4710 5.6504 95.9400 1.8000 0.3285 1.0000 0.1000 6.0000 13.0000 45.0000 4.0000 0.0000 0.0000 0.6062 6.1484 0.0000 0.0000 0.0000 152.6300 3.7659 0.0689 2.9902 0.8563 0.0000 -16.7660 3.1072 1.0338 8.0000 3.4590 0.0000 0.0000 0.0000 1 ! Nr of bonds; at1;at2;De(sigma);De(pi);De(pipi);p(be1);p(bo5);13corr;n.u.;p(bo6),p(ovun1) p(be2);p(bo3);p(bo4);n.u.;p(bo1);p(bo2) 1 2 113.2223 29.8045 86.9949 0.4438 -0.2465 1.0000 16.6234 1.0764 0.1199 -0.1696 7.0000 1.0000 -0.1374 5.4732 1.0000 0.0000 1 ! Nr of off-diagonal terms; Ediss;Ro;gamma;rsigma;rpi;rpi2 1 2 0.1453 1.8849 10.8505 1.6926 1.5319 1.6459 1 ! Nr of angles. at1;at2;at3;Thetao,o;p(val1);p(val2);p(coa1); 1 2 1 83.1375 10.0000 3.5278 0.0000 3.5962 0.0000 2.2128 0 ! Nr of torsions. at1;at2;at3;at4;;V1;V2;V3;p(tor1);p(cot1);n 0 ! Nr of hydrogen bonds;at1;at2;at3;Rhb;Dehb;vhb1 |
木蟲 (著名寫手)
木蟲 (正式寫手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿福大288有機化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 321求調(diào)劑 +7 | 何潤采123 2026-03-18 | 7/350 |
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
[考研] 【同濟軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[基金申請] 被我言中:新模板不強調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|