24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3484  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料292调剂 +12 橘颂思美人 2026-03-23 12/600 2026-03-27 15:44 by caszguilin
[考研] 材料与化工(0856)304求B区调剂 +7 邱gl 2026-03-27 7/350 2026-03-27 15:27 by Joe率
[考研] 考研调剂 +10 呼呼?~+123456 2026-03-24 10/500 2026-03-27 11:46 by wangjy2002
[考研] 315调剂 +4 0860求调剂 2026-03-26 5/250 2026-03-27 11:23 by wangjy2002
[考研] 一志愿郑大085600,310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[考研] 324求调剂 +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-26 4/200 2026-03-27 08:04 by chemisry
[考研] 304材料求调剂 +4 钟llll 2026-03-26 4/200 2026-03-27 03:42 by wxiongid
[考研] 352求调剂 +4 大米饭! 2026-03-22 4/200 2026-03-26 16:40 by 不吃魚的貓
[考研] 289求调剂 +17 硕星赴 2026-03-23 17/850 2026-03-26 16:18 by 不吃魚的貓
[考研] 一志愿 南京邮电大学 288分 材料考研 求调剂 +3 jl0720 2026-03-26 3/150 2026-03-26 13:39 by zzll406
[考研] 化学调剂 +6 yzysaa 2026-03-21 6/300 2026-03-25 09:27 by aa331100
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
信息提示
请填处理意见