版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

laughingjudy

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 185.8
帖子: 33
在線: 31.5小時
蟲號: 2738781
注冊: 2013-10-20
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

[求助] 求幫忙看一個TIP4P水分子in文件，提示Invalid atom type in create_atoms mol command

IN文件如下：
# 3d TIP4P water model
units       real
variable    T equal 298.16
variable    P equal 1
atom_style full

dimension    3
boundary    p p p
molecule    H2O h2o.txt

lattice    fcc 8 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1
region       box block 0 4 0 4 0 4
create_box 1 box bond/types 1 angle/types 1
create_atoms  1 box mol H2O 4928458

mass       3 15.9994
mass       4 1.008
bond_style harmonic
bond_coeff 1 450 0.9572
angle_style harmonic
angle_coeff 1 55 104.52
pair_style lj/cut/tip4p/cut 1 2 1 1 0.15 9 8.5
pair_coeff 3 3 0.1550 3.1536

neighbor    2.0 bin
neigh_modify delay 5 every 1

velocity    all create $T 102486 dist gaussion
timestep    0.01
thermo       1000
fix 1 all shake 0.0001 100 1000 b 1 a 1
fix 2 all npt temp $T $T 100 iso $P $P 1000 drag 0.2

run       1000

h2o.txt如下：
#water

3 atoms
2 bonds
1 angles

Coords

1 0 0 0
2 0.9572 0 0
3 -0.9572*sin14.52 -0.9572*cos14.52 0

Types

1 1
2 2
3 2

Charges

1 -1.040
2 0.520
3 0.520

Bonds

1 1 1 2
2 1 1 3

Angles

1 1 2 1 3

Shake Flags

1 1
2 1
3 1

Shake Atoms

1 1 2 3
2 1 2 3
3 1 2 3

Shake Bond Types

1 1 1 1
2 1 1 1
3 1 1 1

log文件提示如下：
LAMMPS (10 May 2014-ICMS)
WARNING: OMP_NUM_THREADS environment is not set. (../comm.cpp:100)
  using 1 OpenMP thread(s) per MPI task
# 3d TIP4P water model
units       real
variable    T equal 298.16
variable    P equal 1
atom_style full

dimension    3
boundary    p p p
molecule    H2O h2o.txt

lattice    fcc 8 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1
Lattice spacing in x,y,z = 8 8 8
region       box block 0 4 0 4 0 4
create_box 1 box bond/types 1 angle/types 1
Created orthogonal box = (0 0 0) to (32 32 32)
  1 by 1 by 1 MPI processor grid
create_atoms  1 box mol H2O 4928458
ERROR: Invalid atom type in create_atoms mol command (../create_atoms.cpp:166)

不明白是少了語句呢還是哪里寫錯了。請各位大神幫忙解答啊

回復(fù)此樓

» 本帖@通知

@御劍江湖

» 猜你喜歡

085600 286分材料求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
275求調(diào)劑已經(jīng)有10人回復(fù)
286求調(diào)劑已經(jīng)有6人回復(fù)
083000學(xué)碩274求調(diào)劑已經(jīng)有5人回復(fù)
【求調(diào)劑】085601材料工程專碩 | 總分272 | 已經(jīng)有5人回復(fù)
275求調(diào)劑已經(jīng)有13人回復(fù)
0856，材料與化工321分求調(diào)劑已經(jīng)有8人回復(fù)
085701環(huán)境工程，267求調(diào)劑已經(jīng)有12人回復(fù)
320分，材料與化工專業(yè)，求調(diào)劑已經(jīng)有7人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2014-05-28 20:47:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

求助：運(yùn)行了gromacs，想查看體系的總原子數(shù)和水分子數(shù)！已經(jīng)有4人回復(fù)
求教給配位水分子加氫原子的方法！已經(jīng)有6人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 laughingjudy 的主題更新