| 5 | 1/1 | 返回列表 |
| 查看: 1535 | 回復: 3 | ||
| 【懸賞金幣】回答本帖問題,作者callen123將贈送您 10 個金幣 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
callen123金蟲 (小有名氣)
|
[求助]
求助LAMMPS混合勢函數設置 已有1人參與
|
|
|
我的體系是水溶液含ZnOTF2(Zn,C,H,O,N,S)分子的混合體系。 想用reaxff勢函數描述ZnOTF2分子的CHONS元素,用文獻報告的參數描述Zn2+,再用SPC描述H2O分子。但是運行報如下錯誤,感覺是pair_style錯了: ------------- LAMMPS (23 Jun 2022 - Update 4) OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (src/comm.cpp:98) using 1 OpenMP thread(s) per MPI task Reading data file ... orthogonal box = (0 0 0) to (20 20 20) 4 by 4 by 4 MPI processor grid reading atoms ... 722 atoms read_data CPU = 0.074 seconds ERROR: Illegal pair_style command (src/REAXFF/pair_reaxff.cpp:187) Last command: pair_style hybrid/overlay reaxff lj/cut/coul/long 9.0 ------------- 輸入in文件主要部分如下: # Initialization units real dimension 3 boundary p p p atom_style charge # Atom definition read_data ZnOTF2_H2O.data # Force field settings pair_style hybrid/overlay reaxff lj/cut/coul/long 9.0 # Mass and atom types mass 1 32.065 # S mass 2 15.9994 # O mass 3 12.0107 # C mass 4 1.00794 # H mass 5 65.409 # Zn mass 6 1.00794 # H mass 7 15.9994 # O # Define other masses and atom types as needed set type 1 charge 0 set type 2 charge 0 set type 3 charge 0 set type 4 charge 0 set type 5 charge 0 set type 6 charge 0.4238 set type 7 charge -0.8476 # Force field parameters for molecule pair_coeff * * reaxff_JPCA S O C H fix reax all qeq/reaxff 1 0.0 10.0 1.0e-6 reaxff # Define other pair coefficients as needed # Force field parameters for Zn2+ pair_coeff 5 5 lj/cut/coul/long 0.046 1.568 pair_coeff 5 2 lj/cut/coul/long 0.046 2.367 pair_coeff 5 6 lj/cut/coul/long 0.046 0.0 # Oxygen-Zinc interactions pair_coeff 5 7 lj/cut/coul/long 0.046 2.367 # Hydrogen-Zinc interactions # SPC Water model pair_coeff 7 7 lj/cut/coul/long 0.1553 3.166 #8=Ow,7=Hw pair_coeff 7 6 lj/cut/coul/long 0.0 0.0 pair_coeff 6 6 lj/cut/coul/long 0.0 0.0 pair_coeff 6 1 lj/cut/coul/long 0.0 0.0 pair_coeff 6 2 lj/cut/coul/long 0.0 0.0 pair_coeff 6 3 lj/cut/coul/long 0.0 0.0 pair_coeff 6 4 lj/cut/coul/long 0.0 0.0 pair_coeff 7 1 lj/cut/coul/long 0.0 0.0 pair_coeff 7 2 lj/cut/coul/long 0.0 0.0 pair_coeff 7 3 lj/cut/coul/long 0.0 0.0 pair_coeff 7 4 lj/cut/coul/long 0.0 0.0 請大神指點,非常感謝! |
鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 復試調劑 +6 | 呼呼?~+123456 2026-03-08 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 歡迎211本科同學,過A區(qū)國家線,A區(qū)非偏遠一本,交叉學科課題組 +22 | lisimayy 2026-03-04 | 35/1750 |
|
|
[考研] 304求調劑(085602一志愿985) +4 | 化工人999 2026-03-09 | 4/200 |
|
|
[考研] 293求調劑 +4 | 上班不著吉 2026-03-09 | 4/200 |
|
|
[考研] 0817化學工程319求調劑 +5 | lv945 2026-03-08 | 6/300 |
|
|
[考研] 0703化學求調劑 +4 | 很老實人 2026-03-09 | 4/200 |
|
|
[考研] 337求調劑 +3 | 睡醒,。 2026-03-09 | 3/150 |
|
|
[考研] 求調劑,一志愿華中科大0702,數一英一,293 +4 | 小羅露一二 2026-03-07 | 4/200 |
|
|
[考研] 337一志愿華南理工材料求調劑 +4 | mysdl 2026-03-07 | 4/200 |
|
|
[考研] 081700學碩一志愿北京化工大學數二英一過六級有競賽求調劑 +5 | galaxary 2026-03-07 | 7/350 |
|
|
[考研] 一志愿211 化學305分求調劑 +3 | 0703楊悅305分 2026-03-05 | 3/150 |
|
|
[考研] 一志愿武漢理工085601,初試301分,請問能調劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研] 085600材料與化工 298 調劑 +11 | 小西笑嘻嘻 2026-03-03 | 11/550 |
|
|
[考研] 085600材料與化工(高分子)290分求調劑 +9 | wengyujian 2026-03-04 | 9/450 |
|
|
[考研] 268求調劑 +4 | 劉合華 2026-03-05 | 4/200 |
|
|
[考研] 材料與化工304求調劑 +7 | 邱gl 2026-03-05 | 10/500 |
|
|
[考研] 材料考研339求調劑 +3 | Karry*^_^* 2026-03-04 | 3/150 |
|
|
[考研] 334求調劑 +3 | Trying] 2026-03-03 | 4/200 |
|
|
[考研] 學碩材料275調劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
[考研] 295求調劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|