24小时热门版块排行榜    

查看: 3583  |  回复: 4

无事忙0930

金虫 (小有名气)

[求助] lammps我的data文件和in文件如下不知道哪里出错,不能运算,求教虫友细心阅读指教! 已有1人参与

data????????
LAMMPS data file. CGCMM style. atom_style full generated by VMD/TopoTools v1.5 on Tue Dec 06 08:29:29 +0800 2016
15140 atoms
10080 bonds
5080 angles          #?????0???5080#
0 dihedrals
0 impropers
4 atom types
1 bond types         #?????1???2#
1 angle types        #?????0???2#
0 dihedral types
0 improper types
0.019988 272.019988  xlo xhi
0.047490 272.047490  ylo yhi
0.108994 272.108994  zlo zhi

# Pair Coeffs
#
# 1  H
# 2  Na
# 3  O
# 4  S

Bond Coeffs

1 450 0.9572
2 364 0.1430  

Masses           #Masses????VMD????data????????г???????????????????#

1 1.0080000 # H
2 22.989771 # Na
3 15.999000  # O
4 32.060000 # S

Atoms

1 1 1 0.000000 104.268997 209.776993 98.169998 # H HOH
2 1 1 0.000000 104.505997 210.947006 99.283997 # H HOH
3 1 3 0.000000 103.817001 210.438995 98.767998 # O HOH
....................................................................................................

15140 40 2 0.000000 195.259003 201.108994 79.967003 # Na MOL

Bonds

1 1 1 3
2 1 2 3
3 1 4 6
..................................................

10080 1 15096 15097

Angles
1 1 1 3 2
2 1 4 6 5
..................................................

5000 1 14998 15000 14999
5001 1 15001 15003 15002
..................................................

5080 1 15238 15240 15239
end
in????????
echo screen
dimension     3
boundary        p p p
units real                                    
atom_style   full
bond_style   harmonic
bond_style   harmonic                      #?????#
angle_style  charmm
angle_style  charmm                        #?????#
pair_style    lj/charmm/coul/long 10 12      
pair_style    lj/charmm/coul/long 10 12    #?????#
pair_modify    mix arithmetic                 
kspace_style   pppm 1.0e-4               
read_data   waterliusuana.data
set type 1 charge 0.415                    #H?????2#
set type 3 charge -0.830                   #O?????4#
set type 4 charge  -1                      #Cl?????1#  #s#
set type 2 charge   1                      #Na?????3#
angle_coeff  1 55.00000 104.5200 0.0 0.0
angle_coeff  2 55.00000 104.5200 0.0 0.0           #?????????#
pair_coeff   1 1 0.102  3.188                      #???1 cl 2 h 3 na 4 o# #????1 h 2 na 3 o 4 s#
pair_coeff   4 4 0.4184 4.401                      #????????4 s#
pair_coeff   4 2 0.4184 3.492                      #????????4 s#
pair_coeff   4 3 0.5216 3.785                      #????????4 s#         #I<J#
pair_coeff   2 2 0.4184 2.583
pair_coeff   2 3 0.5216 2.876
pair_coeff   3 3 0.6502 3.169
group water type  1  3                             #???#
group liusuanna  type  2  3  4                     #???#
group waterliusuanna type  1 2 3 4                 #???#
velocity    all create 770.25 770099                             
neighbor   2.0 bin
neigh_modify     delay 0 every 1         
min_style    cg
minimize    1.0e-4 1.0e-6 500 1000 # etol ftol maxiter maxeval   stopping tolerance for energy (unitless)  for force (force units)  iterations  max number of force/energy evaluations#         
fix   1 all shake 1.0e-6 500 0 b 1 a 1   # tol=accuracy tolerance of SHAKE solution  iter = max # of iterations in each SHAKE solution #
fix   2 all npt temp 770.25 770.25 100.0 iso 27.78 27.78 1000.0         #gas yong npt #
thermo_style multi
thermo   100
timestep    1.0

dump mydump all custom 1000 dump.*.lammpstrj id type x y z
dump        1  all custom 1000 dump.lammpstrj id type x y z
dump                2 all image 10000 image.*.jpg type type &
                zoom 1.6 adiam 1.5
dump_modify        2 pad 4

dump                3 all movie 10000 movie.mpg type type &
                zoom 1.6 adiam 1.5
restart  100000 water.res
run 2000000

???д????????
F:LAMMPS 64-bit 20150904inguoshuailiusuanna>mpiexec -n 10 lmp_mpi -in
waternacl.in
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
LAMMPS (4 Sep 2015-ICMS)
WARNING: OMP_NUM_THREADS environment is not set. (../comm.cpp:90)
  using 1 OpenMP thread(s) per MPI task
dimension     3
boundary        p p p
units real
atom_style   full
bond_style   harmonic
bond_style   harmonic                      #需??改#
angle_style  charmm
angle_style  charmm                        #需??改#
pair_style    lj/charmm/coul/long 10 12
pair_style    lj/charmm/coul/long 10 12    #需??改#
pair_modify    mix arithmetic
kspace_style   pppm 1.0e-4
read_data   waterliusuana.data
Reading data file ...
ERROR on proc 0: Cannot open file waterliusuana.data (../read_data.cpp:1737)

job aborted:
rank: node: exit code[: error message]
0: kejilou709-PC: 1: process 0 exited without calling finalize
1: kejilou709-PC: 123
2: kejilou709-PC: 123
3: kejilou709-PC: 123
4: kejilou709-PC: 123
5: kejilou709-PC: 123
6: kejilou709-PC: 123
7: kejilou709-PC: 123
8: kejilou709-PC: 123
9: kejilou709-PC: 123@smutao@oxox6085
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yingxong2020

铜虫 (小有名气)

【答案】应助回帖

最后写的是读data文件错误,所以就是data格式不对吧?
2楼2017-01-20 16:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cly001

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by yingxong2020 at 2017-01-20 16:48:16
最后写的是读data文件错误,所以就是data格式不对吧?

read-data格式不对吧,后面应该是data.mubiao
3楼2017-05-15 19:46:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyfwdkkk

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by yingxong2020 at 2017-01-20 16:48:16
最后写的是读data文件错误,所以就是data格式不对吧?

大神你好,我想问一下,我想创建一个纳米团簇的模型,先region一个半径为100的大球,在region一个半径为10的小球,然后group这个小球,我想知道怎么用delete_atoms命令删除小球周围的那部分大球,使小球周围处于真空状态,在小球模型内进行模拟?求救啊大神
4楼2017-05-20 16:59:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HSUTB426

新虫 (初入文坛)

我也遇到这样的问题了 请问楼主解决了吗
5楼2017-07-12 17:20:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 无事忙0930 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +23 低垂的野花 2026-07-31 31/1550 2026-08-07 12:22 by speculiar
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode变化情况 +3 布布和一二 2026-08-07 9/450 2026-08-07 12:19 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +3 yang182083 2026-08-06 4/200 2026-08-07 11:07 by 虫友是什么虫
[基金申请] 基金中了 +11 laoda193707 2026-08-06 11/550 2026-08-07 10:54 by TT和NN
[基金申请] 听说今天filecode变了 +21 布布和一二 2026-08-06 44/2200 2026-08-07 10:11 by 布布和一二
[基金申请] 固定端突然变了,今天 +4 archvillain 2026-08-06 8/400 2026-08-07 08:06 by 带你看大海
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 6/300 2026-08-06 20:19 by 苏知砚
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 12/600 2026-08-06 15:11 by eddyin
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见